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A method for measuring half-lives with Geiger counters which is independent of counting 
losses is described. Foils which have been irradiated for various times are counted after a lapse 
of time adjusted to make the counting rates of the foils the same. Saturated activities of the 
foils will be different if an incorrect value of the decay constant has been assumed. Equations 
are given for obtaining a more correct value of the decay constant. The method is illustrated by 


data on the 54-minute period of indium. 





INTRODUCTION 


EASUREMENT of the half-life of an arti- 
ficial activity with a counter has usually 
required a correction for counting losses.! This 
correction is time-consuming and is not easy to 
make with precision.? In the following a method 
of making such measurements which doés not 
require this correction is reported. 


THEORY 


If the total number of counts N obtained by 
counting the emission of a sample for ¢, minutes 
is measured, the saturated activity A, (the initial 
counting rate of a sample irradiated for an 
infinite time) is given by the formula, 


A,=XN exp (Ats)/(1—exp (—)) 


* Now on leave from University of Texas. 

** Now at Cornell University. 
- 1See for example: H. Volz, Zeits. f. Physik 93, 539 
(1935); L. I. Schiff, Phys. Rev. 50, 88 (1936); W. B. 
Lewis, Proc. Camb. Phil. Soc. 33, 549 (1937); A. E. Ruark 
and F. E. Brammer, Phys. Rev. 52, 322 (1937); L. Alaoglu 
and N. Smith, Jr., Phys. Rev. 53, 832 (1938); J. Kurbatov 
and H. B. Mann, Phys. Rev. 68, 40 (1945). 

? H. Lifschutz and O. Duffendack, Phys. Rev. 54, 714 
(1938); J. Cork and J. Lawson, Phys. Rev. 56, 291 (1939); 
D. E. Hull and H. Seelig, Phys. Rev. 60, 553 (1941). 
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where ¢; is the time the sample was irradiated, 
ts; is the time between the irradiation and the 
start of the counting period, \=In 2/T, and T 
is the half-life in minutes. If the derivative of A, 
with respect to A is calculated and counting 
times selected to make this derivative zero, the 
value of A, will not depend strongly on a small 
error in \. Then if other counting schedules are 
selected with which the derivative is not zero, 
the calculated value of A, will differ from the 
above value if an incorrect value of \ has been 
assumed. Differentiation of Eq. (1) yields the 
formula, 








dA,/dk=KA,, (2) 
where 
ty exp (—At1) 
K= 1/A+4;— 
1 —exp (—At) 
—XMM 
* tz exp ( 2) - @) 
1—exp (—Af2) 


If one assumes that A, is approxiniately a linear 
function of \ one can write 


h\=)'+(A,—A,’)/K’'A,’, (4) 
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. Fic. 1. Schematic representation of calculated saturated 
activity vs. decay constant. 


where )’ is the assumed decay constant, A, is 
the correct saturated activity, A,’ is the saturated 
activity calculated from Eq. (1) with the assumed 
value of decay constant, and d is a more nearly 
correct value of the decay constant. 

In order for this method to be useful enough 
counts must be recorded to make statistical 
errors small, K must have predetermined values, 
and counting rates with the various schedules 
must be nearly the same. Since there are three 
conditions to be fulfilled and three parameters, 
t, and ¢t, and #¢s to adjust, appropriate schedules 
can always be obtained. 

Since all foils have the same activity when 
counted, the counting loss correction will be the 
same for all foils and will not affect the value 
of \ as given by Eq. (4). From Eq. (2) it can be 
seen that the sensitivity of the method (the 
percent difference in A, divided by the percent 
error in \) is equal to KX. Hence, for example, 
for indium where \=0.0128 a value of about 80 
for K is required to give a sensitivity of one. 
If one is attempting to detect an error in \ of 
one percent with a counting schedule adjusted 
to make K =80, one must measure activities to 
one percent. Figure 1 shows a schematic repre- 
sentation of A, as a function of for three 
different counting schedules. Since A, cannot 
vary with the counting schedule if the correct 
value of X is used in calculating saturated ac- 
tivities, all such curves must intersect at one 
point as shown. In practice, because of experi- 
mental errors, the three curves will not be con- 
current but will give three points of intersection 
and three values of \. If Eq. (4) is applied to any 
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resulting equations one obtains, 
N=N/+(A,”—A,')/(K’'A,'—K"A,"). (5) 


A,” and A,’ are saturated activities calculated 
with the assumed decay constant 0’ for counting 
schedules corresponding to K’ and K’. 


HALF-LIFE OF INDIUM 


A test of this method was made with indium.’ 
Data were taken on two Geiger counters with 
three different counting schedules. The three 
different schedules and corresponding values of K 
are given in Table I. K is calculated from 
Eq. (3) with \=0.012836. The times 4), te, and fs 
have the same significance as in Eq. (1). Table I 
gives the calculated activities corresponding to 
the three values of K. N; and A, are the number 
of foils measured with counter 1 and the calcu- 
lated activities. N2 and A,2 are the corresponding 
quantities measured on counter 2. A, gives the 
average activity. A,; and A, have been multi- 
plied by appropriate factors to make the sensi- 
tivities of the two counters the same. Using the 


TaB.eE I. Calculated values of K for three 
counting schedules. 











K hh ts ts 
min. min. min. min. 
—66 9 10 3 

0 20 10 63 
+73 52 10 123 





TABLE II. Activities calculated for three values of K. 








K A al A a A s 
min. M1 counts/min. counts/min. counts/min, 


N: 
—66 13 34882431 . 34867+ 97 34875+50 
7 





0 9 35002+75 349734136 34988477 
73 6 34904+90 349484136 34926+90 








TABLE III. Half-life of indium. 











K’; K” Half-life T 
min. min. min. 
—66; 0 0.012787 +0.000040 54.21+0.17 
0; 73 0.012860 +0.000047 53.90+0.20 
—66; 73 0.012825 +0.000021 54.05 +0.09 


mean 54.05+0.16 








*See also; E. Amaldi, O. d’Agostino, and E. Segré, 
Ricerca Scient. 5, 2 (1934); E. Amaldi, O. d’Agostino, 
E. Fermi, B. Pontecorvo, F. Rasetti, and E. Segré, Proc. 
Roy. Soc. 149, 522 (1935). 


two such curves and A, eliminated from the 
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three saturated activities in pairs one obtains 
from Eq. (5) values for the half-life of indium 
shown in Table III. In this table errors given 
are probable errors calculated from the root 
mean square deviations of the data. 


CONCLUSION 


These measurements were undertaken with the 
primary purpose of determining whether or not 
the half-life of indium was sufficiently different 
from 54 minutes to produce a significant sys- 
tematic error in experiments being performed 
with indium as a neutron detector. An error of 
a percent or more in T was considered sufficient 
to affect these experiments. Since the authors 
were not looking for errors of the order of a few 
tenths of a percent the measurement suffers 
from two uncertainties. First, while the indium 
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used was known to be quite pure, no check of the 
purity was made. Hence, it is possible that an 
impurity with a large cross section may have 
been present and affected the result. Second, no 
attempt was made to time operations with the 
precision necessary to obtain an accuracy of a 
few tenths of a percent. It is believed that the 
result should be correct to one percent or better 
and that the half-life of indium is 54.05+0.5. 
However, it is felt that the method is useful and 
it is hoped the measurements reported give an 
indication of the accuracy which may be obtained 
by this method of determining the half-life of 
an artificially induced activity. 

This work was carried out under contract 
between the University of California and the 
Manhattan District, Corps of Engineers, War 
Department. 
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(Received October 4, 1946) 


Measurements have been made of the yields of photo- 
fission in uranium and thorium together with a search for 
photo-fission in other heavy elements, using continuous 
x-rays from a 100-Mev betatron. Fission was detected in 
the presence of an intense background of x-rays by a dif- 
ferential ionization chamber and linear amplifier, the 
substance investigated being coated on an electrode of one 
chamber. A Victoreen r-thimble, surrounded by }-inch 
lead walls, was used to monitor the radiation. Curves were 
obtained of the number of fissions per roentgen unit for 
uranium and thorium. These are of similar shape, the 
uranium curve showing a rapid rise with increasing x-ray 
energy up to 18 Mev, followed by a gradual decrease as the 
maximum energy of the x-rays is further increased; the 
yield of fissions per roentgen at 100 Mev is about half that 
at 18 Mev. The ratio of uranium and’ thorium yields is 
very nearly two at all x-ray energies. No fissions were 


INTRODUCTION 


CCORDING to Bohr and Wheeler’s theory 
of the fission process,’ fission should be 


Fad pnesormmery! report of this work was presented in an 
invi r at the January, 1946 meeting of the American 


Physica iety. 


1N. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 


observed in intense 100-Mev irradiations of Bi, Pb, Tl 
Au, W, and Sm. Determination of cross sections from the 
yield curves is complicated by the continuous spectrum of 
the x-rays which has not been measured experimentally. A 
rough analysis of the data has been made in which a spec- 
trum is assumed for which the intensity is constant in each 
unit energy interval and the r-meter efficiency calculated 
roughly from a simplified picture of the generation of 
secondaries in the lead walls. The resulting analysis of the 
yield curves shows that the cross section for photo-fission 
as a function of quantum energy passes through a maxi- 
mum and then decreases and is extremely small above 30 
Mev. The maximum cross section is of the order of 
5X 10-** cm? for uranium and half that for thorium. In the 
other elements studied, the cross section must be below 
10-** cm?. 






possible for all heavy nuclei which lie well 
beyond the minimum of the packing fraction 
curve, provided sufficient excitation is provided 
to produce the necessary deformation of nuclear 
fluid which precedes division of the nucleus. 
Such excitation can be provided by particle 
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capture or by absorption of a y-ray into the 
nucleus. 

Shortly after the discovery of fission, Roberts, 
Meyer, and Hafstad* attempted, without success, 
to induce fission in uranium by bombardment 
with the 6.3-Mev y-rays from the F'%(p, vy) 
reaction. The negative result of this experiment 
was explained by Bohr and Wheeler,! who cal- 
culated that the cross section should be of the 
order of 10-?”? cm*. Using higher intensities, 
Haxby, Shoupp, Stephens, and Wells* later 
succeeded in observing photo-fission of uranium 
and thorium with 6.3-Mev y-rays from fluorine. 
The measured cross sections at 6.3 Mev were 
(3.51.0) 10-7 cm? for uranium, (1.7+0.5) 
X10-*7 cm? for thorium. The ratio of uranium 
and thorium cross sections was practically two. 
These results were confirmed by a group of 
Japanese investigators.‘ 

While natural y-rays have the advantage that 
only simple line spectra of radiation are in- 
volved, sources of -rays of sufficient energy and 
intensity for accurate photo-fission measure- 
ments are practically non-existent. Using con- 
tinuous x-rays with maximum energies up to 
20 Mev, Koch® has measured the yields of 
photo-fission in uranium and thorium and sought 
unsuccessfully for photo-fission in other ele- 
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Fic. 1. Differential ionization chamber for detection of 
fission in the presence of intense x-radiation. The material 
investigated for fission is coated on plate C. 


? R. B. Roberts, R. C. Meyer, and L. R. Hafstad, Phys. 
Rev. 55, 417 (1939). 
$7 Gna Shoupp, Stephens, and Wells, Phys. Rev. 59, 
1941). 
‘Arakatu, Vemura, Sonada, Shimizu, Kimura, and 
Kuraoka, Proc. Phys.-Math. Soc. Japan 23, 440 (1941). 


5H. W. Koch, Ph.D. Thesis, University of Illinois, 1944 — 


(unpublished). 
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ments. The failure to observe fission in lead and 
other heavy elements is somewhat surprising, 
since the critical energies for fission calculated 
according to the theory of Bohr and Wheeler! 
are well below the 20 Mev accessible to Koch. 
We have accordingly undertaken a search for 
photo-fission at still higher energies in a number 
of heavy elements, and have investigated the 
energy dependence of fission in uranium and 
thorium. 


APPARATUS 


The source of radiation employed in the 
present experiments was the 100-Mev betatron 
of the General Electric Research Laboratory.*® 
This generates a continuous spectrum of quanta 
of any desired maximum energy up to 100 Mev, 
the radiation appearing in pulses of short dura- 
tion with a 60-cycle repetition rate. 

Fission processes were detected by observing 
the ionization produced by fission fragments 
originating in a thin layer of material deposited 
on the inner wall of an ionization chamber. 
Because individual x-ray pulses from the beta- 
tron are very intense, it was necessary to use a 
balanced differential chamber. It was possible 
by balancing a double chamber to reduce the 
x-ray ionization pulses to the same level as the 
natural background noise of the amplifier; 
ionization produced by particles from material 
coated on the wall of one chamber could then 
readily be detected in the presence of intense 
x-radiation. 

Figure 1 shows details of the ionization 
chamber. The aluminum plates A, B, and C 
defined the two chambers, each 1 cm deep. The 
center plate, A, was connected to the grid of the 
first stage of a linear amplifier. Plates B and C 
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Fic. 2. Housing arrangement of ionization chamber and 
first stage of the linear amplifier. 


*W. F. Westendorp and E. E. Charlton, J. App. Phys. 
16, 581 (1945). 






~ —=— -_—- eee 














280 


260 


Fic. 3. X-ray fission yield 
for uranium. The ordinate is 
the number of fission Dj 


per roentgen unit recorded by 


the monitor. 


FISSIONS PER R 





PHOTO-FISSION IN HEAVY ELEMENTS 













1 1 ) 








were connected to individually adjustable 1000- 
volt stabilized sources of negative and positive 
potential, respectively ; balancing of the chambers 
was greatly facilitated by having these poten- 
tials variable. Plate C was readily interchange- 
able with other plates on which various sub- 
stances to be investigated were coated.’ The 
chamber was filled with air at atmospheric 
pressure. | 

Figure 2 shows the housing arrangement of the 
ionization chamber and the first stage of the 
amplifier. The entire structure was made very 
rugged to eliminate microphonic disturbances, as 
it was necessary to operate the chamber close to 
the betatron where the sound level was about 
100 decibels. The chamber housing was mounted 
on pads of sponge rubber two inches thick to 
eliminate mechanical vibration. 

Pulses from the ionization chamber were am- 
plified by a five-stage linear amplifier similar to 
that described by Waddell.* The output was fed 
to a vacuum-tube discriminating and thyratron- 
recording circuit. The output of the third stage 
of the amplifier was observed visually with an 
oscilloscope throughout each run to insure that 
pulses which were registered by the discrimi- 


7 Remarkably uniform coatings were deposited for us 
by Messrs. Dunbar, Huthsteiner, and Young of this 
laboratory under the supervision of Dr. Dushman. 
® Waddell, R. C., Rev. Sci. Inst. 10, 311 (1939). 
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nating circuit occurred only at the phase at which 
x-rays were produced; to facilitate this, the 
chamber was always operated slightly un- 
balanced. A slight 60-cycle pick-up from stray 
magnetic fields of the betatron did not interfere 
with operation of the amplifier. 

The radiation incident on the chamber was 
monitored by means’ of a standard Victoreen 
r-thimble and condenser r-meter. The thimble 
was placed in a lead cap of }” wall thickness 
directly behind the chamber. 


EXPERIMENTS 


Preliminary tests were made with a layer of 
thorium oxide on Plate C (Fig. 1) to determine 
the proper value of the discriminator setting in 
order that the recording circuit would register 
only pulses larger than those due to alpha- 
particles. The discriminator setting was carefully 
maintained throughout all subsequent measure- 
ments at such a value that no pulse less than 
four times the maximum height of the alpha- 
particle pulses was recorded as owing to a fission 
fragment. 

The ionization chamber was then placed in 
the center of the x-ray beam 11 feet from the 
target of the betatron; at this point, the diameter 
of the x-ray beam is comparable with the dimen- 
sions of the chamber. With uncoated plates in 
the chamber, an intense irradiation at 100 Mev 
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Fic. 4. X-ray fission yield 
for thorium. 
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was given. No pulses were observed; hence, any 
pulses observed with a coated plate must be 
attributed to the coating substance. 

Using a plate coated with a layer of 4 mg/cm? 
of uranium oxide, measurements were then 
made of the x-ray fission yield (i.e., the number 
of fission pulses observed per roentgen unit of 
radiation as recorded by the r-meter) for 
x-radiation of various maximum energies from 
10 to 100 Mev. For these measurements, the 
x-ray intensity was maintained at a sufficiently 
low value so that the chance of two fission pulses 
occurring in the same cycle was below one 
percent. The results are’ plotted in Fig. 3, in 
which the observed number of fissions per 
roentgen unit is given as a function of the 
maximum x-ray energy. The x-ray fission yield 
' rises rapidly with increasing x-ray energy up to 
about 18 Mev. As the maximum x-ray energy is 
increased above 18 Mev, however, the yield 
decreases and continues to fall until, at 100 Mev, 
the yield is less than half that obtained at 18 
Mev. The significance of this feature of the 
curve will be explained in the following section. 

After a thorough cleaning, the ionization 
chamber was checked by a run with an uncoated 
target for any possible uranium contamination, 
and a target plate was installed with a coating 
of 9 mg/cm? of thorium oxide. The x-ray fission 
yield for thorium was then determined; it is 
plotted in Fig. 4. This curve closely resembles the 
curve for uranium except that its maximum 
occurs at 20 Mev, slightly higher than for 
uranium. It will be noted that for x-ray energies 
above 25 Mev the yield from uranium is 2.0+0.1 
times the yield from thorium. Below 25 Mev, 


the ratio of the yields decreases slightly, the 
lowest value of 1.75 being reached at the maxi- 
mum of the thorium yield curve, because the two 
maxima are slightly displaced with respect to 
each other. Below 18 Mev, the steepness of the 
curves is such that it is difficult to determine the 
ratio of the yields with precision; however even 
near 10 Mev it appears to be approximately two. 

Upon completion of the work with thorium 
the chamber was again cleaned and checked for 
contamination. A target plate coated with 3 
mg/cm? of lead oxide was then installed and 
given an irradiation of 21 r with 100-Mev 
x-radiation. No fission pulses were observed. 
Another target plate coated with 7.3 mg/cm? of 
metallic lead deposited by evaporation was 
irradiated with 14 r of 100-Mev x-rays and with 
10r of 50-Mev x-rays; no fission pulses were 
observed in either irradiation. 

Similar experiments were carried out using 
other target materials. Substances investigated 
are listed in Table I; except for uranium and 
thorium, none gave any indication of fission. 

Upon the conclusion of these experiments, the 


TABLE I. Substances investigated for photo-fission. 
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thorium target was reinstalled and the excitation 
curve of Fig. 4 was rechecked, indicating that the 
negative results with other substances were not 
owing to failure of the apparatus. 


ANALYSIS OF RESULTS 


In attempting to interpret the excitation 
curves of Figs. 3 and 4, one must remember that 
these curves were obtained with continuous 
x-rays. The fission processes observed at some 
maximum x-ray energy E are caused by quanta 
of all energies between the fission threshold and 
E. It must also be remembered that the ordinates 
are the number of fissions per roentgen unit as 


‘ recorded by the r-meter under the conditions of 


the experiment. This ‘‘x-ray fission yield” is thus 
a ratio of the yields of two entirely different 
processes. Before the photo-fission cross section 
as a function of quantum energy can be deter- 
mined, it is necessary to know the spectral 
distribution of the radiation in the x-ray beam 
from the betatron and to calculate the reading 
in roentgen units of the monitor when it has been 
irradiated with a given intensity of radiation 
having this spectral distribution. Experimental 
knowledge of the spectrum is not yet available 
and the efficiency of the ionization monitor can 
be calculated only approximately. A provisional 
approximate analysis of the yield curves will be 
attempted, however. 

Let a target of area A containing N atoms per 
cm? be irradiated normally with x-rays of maxi- 
mum energy E of intensity such that n(e, E)de 
quanta per cm? of energy between ¢ and e+de 
are incident on the target. Let o(¢) be the cross 
section for photo-fission by a quantum of energy 
e. Then the number of fission processes which 
occur in the target will be 


F(E) -wa f a(e)n(e, E)de. (1) 


The number of fissions per atom of target 
material produced when the target is irradiated 
with continuous x-rays whose quanta have a 
maximum energy E is thus given by the area 
lying under a curve which is the product of the 
cross section for photo-fission by the spectral 
distribution of the quanta. The number of 
fissions observed is divided by the reading of the 
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ionization monitor in r-units and the quotient 
plotted as the x-ray fission yield corresponding 
to “x-rays of energy E.” 


X-Ray Spectrum 


Since the target of the betatron is thin, the 
x-ray spectrum can be expected to be practically 
the same as the spectrum calculated by Bethe 
and Heitler in their theoretical treatment of 
bremsstrahlung.® According to their result, the 


“ INTENSITY 
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Fic. 5. Bremsstrahlung spectrum as calculated by 
Bethe and Heitler. The dotted curve is the simplified 
spectrum assumed for analysis of the fission data. 


intensity per unit energy interval is approxi- 
mately constant except near the high and low 
energy ends of the spectrum. Since the decrease 
in fission yield for E>20 Mev indicates a rela- 
tively small fission cross section above this 
energy, little error will be introduced by “‘squar- 
ing off” the upper end of the spectrum by the 
assumption of a rectangular intensity distribu- 
tion. The spectral distortion introduced by 
absorption of the x-rays in their passage through 
the glass walls of the betatron doughnut and 
other intervening material will-also be neglected 
in the following analysis, since this distortion is 
slight except for quanta of very low energy, for 
which it tends to act as a filter. In Fig. 5, the 
theoretical intensity spectrum of Bethe and 
Heitler is given by the solid curve, and the 
approximate intensity spectrum employed in 
this treatment is shown as the dotted curve. 
The number of quanta per square centimeter 
with energy between « and e+de can thus be 


*H. A. Bethe and W. Heitler, Proc. Roy. Soc. 146A, 83 
Onto W. Heitler, The Quantum Theory of Radiation 
Oxford yop ay bt nen New York, 1944), second edition. 
B. Rossi and K. isen, Rev. Mod. Phys. 13, 240 (1941). 

















or 


os -——+ 
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Fic. 6. Relative response of the monitor (Victoreen 
r-thimble surrounded by }-inch lead walls) to continuous 
x-rays of the spectrum assumed with maximum ene 
given by the abscissa, the intensity being such that the 
constant KX in Eq. (2) is unity. 


written 
n(e, E)\de=Kde/e e<E, (2) 
=0 e>E. 


r-Meter Sensitivity 


The ionization recorded by the r-meter is the 
result of the passage through the thimble of 
negative and positive electron secondaries pro- 
duced by absorption or scattering of x-ray 
quanta in the walls of the 3-mm lead cap. To 
calculate this exactly would require taking 
account of the detailed energy and angular dis- 
tribution of the secondaries, of specific ionization 
as a function of electron energy, of scattering of 
the electrons, of tertiary processes, of the com- 
plicated geometry, and of the exact spectrum 
of the radiation. 

An approximate calculation has been made 
which it is believed describes the relative de- 
pendence of r-meter efficiency on the maximum 
x-ray energy, in which an extremely simplified 
set of assumptions has been used. 

Consider a beam of x-ray quanta with the 
spectral distribution assumed above incident 
upon an absorbing slab of thickness ¢. At a 
depth x in’ the slab, the number of quanta per 
square centimeter with energies between ¢ and 
e+de is given by ~ 


n(e, E, x)de=K exp(—rx)de/e €<E, (3) 
=0 - e>E&, 


where the absorption coefficient 7 is the sum of 
three partial coefficients 7,, 7-, and r, represent- 
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ing, respectively, pair production, Compton 
scattering, and photoelectric effect. Each absorp- 
tion act is assumed to produce two electrons of 
average energy ¢/2 if by pair production, one 
electron of average energy ¢/2 if by Compton 
effect, and one electron of average energy e if by 
photoelectric absorption. These electrons are 
assumed to proceed in the forward direction. 

The wall thickness, ¢, employed in these 
measurements is a fraction of a radiation length 
and is equal to the range of a 5-Mev electron. 
Tertiary processes in the wall can accordingly 
be neglected without introducing significant 
error. 

The number of electrons emerging into the 
space beyond the wall is then 


»(E) = f [(2refetrefetrefe)/elde. (4) 
0 


The factors f take account of absorption of 
quanta in the wall together with the fact that 
an electron of low energy can emerge from the 
wall only if its point of production lies within 
its range of the second surface: 


f= f exp(—rx)dx, 


or 


t 
f= f exp(—rx)dx, 
0 


according to whether or not the range R of a 
given secondary is less or greater than ?¢. In 
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Fic. 7. Analysis of the uranium yield curve of Fig. 4. 
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‘evaluating f, the range of an electron has been 
taken proportional to its energy. 

We now assume that all electrons from the 
wall which enter the thimble produce an equal 
average ionization, §, of 60 ion pairs per cen- 
timeter. The specific ionization of electrons is 
nearly independent of energy for electron ener- 
gies above 300 kev, amounting for air to 76 ion 
pairs per centimeter at 300 kev, decreasing to 
53 at 2 Mev, then rising slowly to 70 ion pairs 
per centimeter at 100 Mev. Low energy electrons 
which ionize appreciably higher than this are 
also highly scattered and are contributed from 
very thin layers of the wall; therefore, the error 
introduced by neglecting their higher ionizing 
power is not serious. The effect of multiple scat- 
tering and of the initial angular distribution of 
produced secondaries is somewhat reduced by 
having the thimble surrounded by the lead on 
all sides. 

The integration of (4) was carried out graphi- 
cally. For E=10 Mev, the number of secondaries 
contributed from all three absorption processes 
was calculated. To this were added the pairs and 
recoils resulting from absorption of quanta above 
10 Mev for a sequence of values of E up to 100 
Mev, the x-ray spectrum in each case having the 
form assumed above and with the constant K 
adjusted to unity. The result of this calculation 
is given in Fig. 6. Absolute accuracy cannot be 
claimed for it, but it is believed that its trend 
with energy closely parallels the actual variation 
of sensitivity of the r-thimble. 

The reading of the r-meter can then be written 


r(E) =e-8-v(E)-K 


=2.9X10-'-»(E£)-K e.s.u. (5) 


It follows that to produce a reading of 1 roentgen 
unit (i.e., one e.s.u. of charge per cubic cen- 
timeter) in the r-meter when the thimble is 
irradiated with 100-Mev x-rays and surrounded 
with }” of lead, there must be approximately 
2.5X10"/e quanta of energy e per unit energy 
interval per square centimeter. 


Analysis of the Excitation Curves 
The ‘‘x-ray fission yield” is the quotient 
®(E) = F(E)/r(£), (6) 
where F(Z) and r(£) are given by (1) and (5). 





PHOTO-FISSION IN HEAVY ELEMENTS 9 








From them it follows that the fission cross 
section is 
2.9X10-8E d(#yr) 


o(£) = ' 7 
) NA dE ”) 





Figure 7 shows the analysis of the data for 
uranium. It is evident that the cross section for 
photo-fission must be extremely small for 
quantum energies above 30 Mev. The resulting 
curve of o(£) vs. E is given in Fig. 8 for uranium. 
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Fic. 8. Cross section for photo-fission of uranium. 
Because of approximations made in the analysis and 
uncertainty of the experimental data for low x-ray energies, 
this curve must be considered of only qualitative sig- 
nificance. 


The determination of the cross section in- 
volves first drawing a smooth curve through the 
experimental points, multiplying the ordinates 
of this curve by a function known only roughly, 
and then differentiating the product. It is evident 
that the given data are not sufficiently precise 
to determine the cross section as a function of 
energy, even if the calculation of r(Z) were 
exact, since in this case the portion of the yield 
curve which is most significant to the deter- 
mination of the cross section is that below 30 
Mev. The experimental points in this region are 
subject to considerable fluctuation, and the true 
yield curve may be appreciably different from 
that shown in Fig. 3. It is also true that the 
effects of the approximations made in calculating 
r(£) are more serious at low energies. Therefore, 
it must be stressed that no significance can be 
attributed to the shape of the curve of o vs. E 
nor to the absolute values of the cross section 
given on the scale of ordinates. The curve of 
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Fic. 9. Four attempts to fit the experimental yield 
curve by assuming simple functions for the photo-fission 
cross section, identified by ‘‘e’’ in each figure. In each case, 
the curve denoted by “®” is the yield curve which should 
result from the assumed cross section with x-rays of the 
simple spectrum assumed in the text. The experimental 
curve is shown dotted. In ‘‘D’’ the cross section assumed is 
a resonance curve with 16-Mev resonance energy and 
3-Mev half-width at half-maximum. 


Fig. 8 is given for purposes of illustration and 
to emphasize the fact that the cross section must 
be very small above 30 Mev. 

The thorium data can be analyzed in the same 
manner. Owing to large fluctuation in the points 
below 40 Mev, it is not possible to compare the 
cross sections for uranium and thorium in detail. 
However, the data indicate similar variation 
with energy, the curve for thorium being dis- 
placed by 1 or 2 Mev owing to its higher thresh- 
old energy for fission. If it is permissible to 
assume that both targets were thick and that the 
average range of a fission fragment from thorium 
is equal to that of a fragment from uranium, 
then it can be concluded that the ratio of the 
photo-fission cross sections for uranium and for 
thorium at corresponding energies is 2.0+0.1 
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The problem of interpretation of the yield 
curves can also be approached in another way. 
One can assume several possible functions for 


@(E) and for each construct the corresponding 


yield curve which would result with an x-ray 
spectrum of the simple form assumed and with 
the r-meter sensitivity varying with maximum 
x-ray energy in the manner calculated above. 
One can then see how sensitive the result ob- 
tained by direct analysis is to the exact shape 
of the yield curve. 

Figure 9 shows several simple functions repre- 
senting possible approximate cross sections and 
the corresponding yield curves. In each figure, 
the experimentally obtained yield curve is 
shown dotted. It will be evident that only an 
extremely small photo-fission cross section can 
be possible if the decrease in yield above 30 Mev 
is to be accounted for on the basis of the assump- 
tions made above. 

According to the theory of Bohr and Wheeler,! 
photo-fission should be energetically possible and 
the fission fragments should have sufficient 
kinetic energy to be readily observable for all of 
the elements listed in Table I with the exception 
of samarium. However, of all the elements which 
were investigated, only uranium and thorium 
gave a positive result. The photo-fission cross 
sections of all other elements tested can hardly 
exceed 10-* times the cross section for uranium; 
ie., they are of the order of 10-** cm? or less. 

The decrease in the photo-fission cross section 
at high energies exhibited by uranium and 
thorium together with the absence of fission in 
the other heavy nuclei investigated is probably 
to be attributed to the competition of other 
types of photo-disintegration.'® 


10 G. C. Baldwin and G. S. Klaiber, Phys. Rev. 70, 259 
(1944). 
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The yield function for the disintegration of lithium bombarded by deuterons has been 


measured with detectors sensitive to neutrons, gamma-rays, and beta-rays. A gamma-ray of 
4.9-Mev energy was found in the course of this investigation. Resonances were found at 0.65, 
1.02, and 1.35 Mev corresponding to excited states of the intermediate nucleus (Be®) at 17.14, 
17.43, and 17.69 Mev. The resonance at 1.35 Mev was found only for the formation of radio- 
active Li* accompanied by the emission of a proton. Apparently proton emission is permitted 


and neutron emission forbidden from the 17.69-Mev state. 





HE yield of Li’, neutrons and alpha- 

particles from lithium bombarded by deu- 
terons of energies up to 1 Mev has been measured 
by Rumbaugh, Roberts, and Hafstad.! Their 
excitation curves showed a smooth increase in 
yield with energy until energies of about 700 kev 
were reached. From 700 to 1000 kev their yield 
curves continued to increase but less rapidly than 
at lower energies. Later experiments? which were 
carried out at Rice Institute on the excitation 
curves of carbon bombarded by deuterons of 
energies up to 2.0 Mev showed a number of reso- 
nances with half-widths of from 10 kev to 100 
kev. These experiments with carbon suggest that 
similar resonances might be found in the dis- 
integration of lithium by deuterons and so a 
similar study was undertaken for lithium bom- 
barded by deuterons; preliminary results of these 
experiments were reported at a meeting of the 
American Physical Society.**4 


APPARATUS 


The Rice Institute pressure Van de Graaff 
machine was used to produce bombarding par- 
ticles of homogeneous and variable energy. Thin 
targets of LiCl, LiOH, and metallic lithium were 
used. They were produced by evaporation from a 


-, Now at Illinois Institute of Technology, Chicago, 
inois. 

t Now at University of Wisconsin, Madison, Wisconsin. 

§ Now at Texas Company, Houston, Texas. 

1L. H. Rumbaugh, R. B. Roberts, and L. R. Hafstad, 
Phys. Rev. 54, 657 (1938). 

*W. E. Bennett, T. W. Bonner, E. Hudspeth, H. T. 
i es and B. E. Watt, Phys. Rev. 59, 781 (1941). 
an. . Bennett and B. E. Watt, Phys. Rev. 60, 167 


asa Bonner and H. T. Richards, Phys. Rev. 60, 167 





tungsten filament onto a polished disk of pure 
silver or copper. The LiCl targets blackened 


‘rapidly during bombardment and the yield of 


11 


neutrons increased with time as a result of the 
deposition of carbon. These targets were made 
thicker (50-100 kev) than the metallic lithium 
targets which were used later so that the carbon 
would have a relatively small effect. They were 
also moved frequently to bombard a new spot. 
The thickness of the targets was estimated from 
the amount of material evaporated, and also, 
somewhat more precisely, by comparing the yield 
from the thin target with that of a thick target of 
the same material. The latter calculation was 
possible after the yield function for the reaction 
was known. It also involved suitable correction 
for the range-energy relationship of the bom- 
barding particles in the thick target. The thick- 
ness of the target was expressed in kev, meaning 
the average energy lost by a bombarding particle 
in passing through it. 

The yields of beta- and gamma-rays were 
measured with Geiger counters filled with a self- 
quenching mixture of argon and alcohol. A scale 
of eight was used in the recording circuit. Coinci- 
dence counting was sometimes used, a bias being 
set to cut out pulses from one counter and record 
only the pulses of double size produced by 
simultaneous discharge of the two counters. To 
count gamma-rays, carbon blocks were used to 
absorb out the beta-radiation. 


THE GAMMA-RAYS FROM Li’+H? 


Only one gamma-ray (of 450-kev energy) has 
been reported® from the deuteron bombardment 


‘J. H. Williams, W. G. Shepherd, and R. O. Haxby, 
Phys. Rev. 52, 390 (1937). 
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of lithium. In the course of these experiments, a 
gamma-ray of higher energy‘ was found.® Its 
energy was determined by inserting aluminum 
sheets between two coincidence counters and 
recording the coincidence counting rate as a 
function of the thickness of aluminum. The curve 
is given in Fig. 1. The maximum range of 
Compton electrons as estimated from this curve 
is 8.70.4 mm of aluminum. This corresponds to 
an energy of 4.9+0.3 Mev. The relation between 


RATE 
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Fic. 1. The rate of coincidence counting as a function 
of the thickness of aluminum between counters. The end 
— at 8.7+0.4 mm ae Fo an energy of 4.9+0.3 Mev 
or the gamma-rays from Li’+H?. 


range and energy had previously been deter- 
mined by measurements on gamma-rays of 
known energy with the same apparatus.’ 

It is very probable that the high energy gamma- 
ray is associated with a neutron group as indi- 
cated by the following equation : 


sLi?-+,H?—,Be®+on!+10.1 Mev, 
*Be*—>,Be®+7+4.9 Mev. 


The spectrum of neutrons from this reaction 
shows evidence of the associated neutron group. 
Both Stephens* and Richards® found a wide 
group of neutrons at about 10.5 Mev. This group 


*W. E. Bennett, T. W. Bonner, H. T. Richards, and 
B. E. Watt, Phys. Rev. 59, 904 (1941). 

™W. E. Bennett, T. W. Bonner, and Bob E. Watt, 
Phys. Rev. 59, 793 (1941). 

*W. E. Stephens, Phys. Rev. 53, 223 (1938). 

°H. T. Ri s, Phys. Rev. 59, 796 (1941). 


is too wide and:of too low energy to be entirely 
accounted for by the 2.8-Mev level in Be® which 
breaks up into two alpha-particles, but it may be 
regarded as consisting of two unresolved groups, 
one associated with the 2.8-Mev level and the 
other with a level at 4 or 5 Mev which emits 
gamma-rays because of a selection rule which 
forbids the break-up into alpha-particles. A simi- 
lar level in Be*® at 17.5 Mev is responsible for the 
high energy gamma-rays from Li’ bombarded by 
protons. 

To detect 4.9-Mev gamma-rays in the presence 
of 450-kev gamma-rays, coincidence counters 
were used, the wall thickness being sufficient to 
prevent detection of gamma-rays of less than 
2-Mev energy. However, this procedure was 
found unnecessary and a single Geiger counter 
was used for the final data. The single counter 
was used to take a lead absorption curve under 
good geometrical conditions and the radiation 
seemed to be homogeneous. The 450-kev radia- 
tion was too feeble to give a noticeable effect 
under our conditions. The single Geiger counter 
did not respond appreciably to neutrons. A block 
of paraffin which reduced the neutron intensity 
to half, did not affect the counting rate in the 
Geiger counter. The neutrons produced radio- 
activity in the counter walls, but this was sub- 
tracted as part of the natural effect. 


THE YIELD CURVE FOR GAMMA-RAYS 
AND NEUTRONS 


The yield of gamma-rays and neutrons as a 
function of energy from 0.5 to 1.9 Mev is given in 
Fig. 2. The neutrons were observed in the 
forward direction. The scale of ordinates is arbi- 
trary and different for the two curves. The 
absolute yield of gamma-quanta was roughly 3 
percent of the yield of neutrons. This figure was 
arrived at by assuming that the Geiger counter 
was 2 percent efficient for counting gamma-rays 
and by using the data of Amaldi, Hafstad, and 
Tuve!® for the yield of neutrons from a thick 
target. 

The similarity of the gamma-ray and neutron 
curves shows that the same intermediate nucleus 
is involved. Since most of the neutrons are from 


' bombardment of the Li’ isotope, it is out of the 


10 FE. Amaldi, L. R. Hafstad, and M. A. Tuve, Phys. 
Rev. 51, 896 (1937). 
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Fic. 2. Yields of neutrons and 


rtially thin” : 


-rays from “ 
nage mbesded by 


targets of lithium 


deuterons. The thickness of the tar- 
gets was 55 kev for the gamma-ray 
curve and 90 kev for the neutron 
curve. The scale of ordinates is arbi- 
trary and different for the two 


curves. 
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question that the gamma-rays might have been 
produced by bombardment of the Li® isotope. 


THE YIELD CURVE FOR BETA-RAYS 
AND NEUTRONS 


Coincidence Geiger counters were used to de- 
tect the beta-rays. The wall thickness discrimi- 
nated against beta-rays of less than 2 Mev 
energy. The arrangement counted some gamma- 
rays. To see how many gamma-rays were 
counted, carbon absorbers were placed between 
the target and the counters to absorb out the 
beta-rays. Less than one percent of the counts 
remained after absorption of all the beta-rays 
from the target. 
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_ Fic. 3. Yields of neutrons and gamma-rays from “thin” 
lithium targets bombarded by deuterons. The yield is in 
arbitrary units, different for each curve. The neutron 
curves have been fitted at 0.55 Mev. The angular distribu- 
tion of neutrons changes with bombarding energy. 
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The yield curves for beta-rays and neutrons 
were taken simultaneously. The methane-filled 
electroscope used to detect neutrons was placed 
at 90° to the bombarding beam. The beta-rays 
were detected in the forward direction. The 
target was made by evaporating lithium metal in 
vacuum, and its thickness was estimated from 
the yield to be 14 kev. 

The excitation curves for beta-rays and for 
neutrons at 90° were run from 0.55 Mev to 1.4 
Mev. They are shown in Fig. 3. Another curve 
was run with the neutron detector in the forward 
direction. This was normalized to agree with the 
90° curve at 0.55 Mev and plotted in Fig. 3. 

The neutron curve shows resonances at 0.65 
and 1.02 Mev. The neutrons show a departure 
from a spherical symmetry which is different for 
the two resonances. If there is spherical sym- 
metry for the first resonance, then the second 
resonance shows a marked maximum of neutron 
intensity in the forward direction. Further 
studies of the angular distribution of the neutrons 
from the reaction would be desirable. 

For beta-rays, the first resonance appears to be 
at 0.75 Mev, about 0.1 Mev higher than the 
corresponding neutron resonance. It may, how- 
ever, correspond to the same state of the inter- 
mediate nucleus. The beta-rays are produced in 
the reaction: 


3Li’+,H?—;Li*+,H!—0.26 Mev, 
3Li8—,Be®+8. 





The yield of beta-rays is equal to the yield of low 
energy protons. The protons must penetrate a 
potential barrier to escape and, therefore, the 
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yield of protons would be expected to differ from 
the yield of neutrons by a factor which increases 
with energy in the low energy region. This factor 
would shift a maximum in the yield curve toward 
higher energies. This shift does not occur for the 
second resonance at 1.02 Mev; but this second 
resonance is narrower and the proton is more 
energetic, so it is less affected by the potential 
barrier. This argument seems to be contradicted 


E. P. NEY AND F. C. ARMISTEAD 


by the fact that the first resonance is stronger 
than the second for proton emission; however, 
other differences in the states of the intermediate 
nucleus may affect the strength of a resonance, A 
third resonance at 1.35 Mev in the beta-ray curve 
corresponds to a state of Be® for which neutron 
emission is forbidden or at least is much weaker 
than for the other resonances at 0.65 and 
1.02 Mev. 
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The Self-Diffusion Coefficient of Uranium Hexafluoride* 


Epwarp P. Ney AND FontTAINE C. ARMISTEAD** 
Rouss Physical Laboratory, University of Virginia, Charlottesville, Virginia 
(Received October 9, 1946) 


A method is reported for measuring the self-diffusion coefficient of a vapor by observing 
with a mass spectrometer the rate at which material enriched in one isotope diffuses into 
normal material in an apparatus of known geometry. A value of pD (where p is the density 
and D the diffusion constant) was obtained for uranium hexafluoride at 30°C. Knowledge of 
the viscosity together with the constants reported here allow calculation of the molecular force 


law. 


hoe the kinetic theory of transport phenomena 
pD is a prominent factor whose value must 
be known. Heretofore its value has been deter- 
mined from the coefficient of viscosity, n, on the 
basis of the equation pD=en where « is at best 
an uncertain quantity. The following is a report 
of the measurement of pD in which D is under- 
stood to-be the diffusion coefficient for one gas 
diffusing into another where the only difference 
between the gases is their isotopic ratios. 

This measurement represents an excellent 
approximation to the self-diffusion coefficient, 
which strictly speaking cannot be measured.! 
If one may speak of the self-diffusion of a gas 
whose molecules are composed of more than one 
isotope, one should probably define it as the 


*This paper is based on work performed under Con- 
tract No. OEMsr-398 with the Office of Research and 
Development at the University of Virginia, in 1942 and 
1943. “The information presented in this paper will appear 
in Division X of the Manhattan Project Technical Series 
as of the contribution of the Physics Department, 
University of Virginia.” 

** Now at Clinton Laboratories, Oak Ridge, Tennessee. 

1 migeomaly the self-diffusion coefficient can be defined 
only for a gas composed of identical molecules, e.g., all 
UF, or all U*F,. 


diffusion of one gas in which the isotopes are in 
their normal ratio into another with normal 
ratio. The ideal experiment to measure this 
would consist of tagging one group of molecules 
differently from the other group in such a way 
that the tagging would not affect the phe- 
nomenon of diffusion. If uranium hexafluoride 
slightly enriched in U?** and uranium hexafluoride 
slightly enriched in U** diffuse into each other 
exactly as normal uranium hexafluoride diffuses 
into normal then isotopic concentration is just 
such an ideal tag. As can easily be seen the 
measurement of the changing concentration of 
one isotope as the diffusion takes place gives 
directly the changing concentration of one group 
in the mixture of both groups of molecules. 


I. DESCRIPTION OF APPARATUS 
AND PROCEDURE 


The apparatus is shown schematically in Fig. 1. 


The diffusion bulbs F and B were two coaxial 
copper cylinders of volumes 925 and 259 cubic 
centimeters, respectively. They were connected 
by the copper pipes C (length 6.80 cm; area 
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Fic. 1. Schematic drawing 
of the diffusion apparatus 
(not to scale). 


0.3215 cm*) and E (length 6.78 cm; area 0.3238 
cm?) and by the Neoprene tubing, D, of area 
0.2714 cm’. In order to obtain different relaxation 
times, the length of the Neoprene was made 2.08 
cm in six of the runs and 6.10 in two of them. 
Connected to bulb F is a null type sylphon 
pressure gauge G. The pressure above the sylphon 
could be made equal to the pressure in the bulbs 
by operating valves R and Q to obtain a null 
reading of the indicator needle which was ob- 
served in a microscope field. The neutralizing 
pressure was read on the mercury manometer NV 
and the oil manometer O, the latter being capable 
of showing a pressure change of 0.06 mm Hg. 

Material from bulb B was drawn through the 
adjustable capillary leak A into the ion source 
of the mass spectrometer. The leak A and the 
valve I add negligible volumes to the diffusion 
bulbs. 

The system used for purifying the samples and 
filling the bulbs was connected to them through 
the valve J. This system consisted of a pressure 
gauge H and a manifold to which samples could 
be connected at L and M. The diffusion system 
was submerged in a 12-gallon water bath which 
was equipped with a vigorous stirrer and with 
a cover to minimize evaporation. A Beckman 
thermometer was placed beside each bulb. 

Measurements of the dimensions of the ap- 
paratus were straightforward with the exception 
of the cross section of the Neoprene tube (D). 
This cross section was decreased by the pressure 
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on the outside and by the thin uranium hexa- 
fluoride reaction coating on the inside. Both of 
these effects were measured. The Neoprene used 
did not undergo any continuous reaction with 
uranium hexafluoride. 

In order to establish that the samples were 
pure after admission to the diffusion bulbs, and 
that they remained pure throughout the course 
of an experiment, the impurity peaks were 
examined in the mass spectrometer. Samples 
with impurities less than one percent by volume 
were obtained for the experiment. 

A certain amount of material was drawn by 
the mass spectrometer. This withdrawal causes 
some mass flow from F to B. An upper limit to 
this effect was determined by the following test. 
A sample of uranium hexafluoride was admitted 
to the diffusion bulbs. The pressure in the bulbs 
and the intensity of the ion peaks in the spec- 
trometer were made the same as those used in 
the diffusion runs. After a careful reading on the 
the gauge (G) had been obtained, the Neoprene 
(D) was pinched and left closed for 15 hours. If 
there had been a rate of withdrawal of such mag- 
nitude as to affect the measurement of the dif- 
fusion constant by as much as one percent there 
would have been an instantaneous and observ- 
able change at the gauge (G) upon unpinching D. 
No change was noted. 

Before making any diffusion runs, the entire 
apparatus was pretreated. When it was necessary 
to replace some part of the apparatus, the pre- 
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treatment was repeated before more runs were 
made. A vacuum was maintained in the bulbs 
when runs were not being made. 

In making the runs, samples were attached at 
L and M and purified. One sample had been 
enriched in U***F., the other in U?*5F,. The first 
sample was then admitted to the diffusion 
system (through valve J) to a pressure of 1 cm 
of mercury. The Neoprene (D) was then pinched 
bulb F was pumped out and filled to 1 cm with 
material from the second sample. The character 


of the first sample was then the character of the 


material in bulb B at the start of the diffusion. 
Only moderate care was taken to remove traces 
of the first sample from F and from the walls of 
the purifying system before introducing the 
second sample. The presence of a slight amount 
of the first sample in F would not affect the 
calculation of pD, but would merely reduce the 
over-all concentration change measured. Before 
starting the diffusion the Neoprene tube was 
unpinched for 10 seconds in order to assure equal 
pressures in F and B. (A) was opened and peaks 
were brought into the mass spectrometer; the 
Neoprene was unpinched, and the diffusion 
started. Continuous measurement of the 235/238 
ratio was made in the mass spectrometer by the 
null method.” This consists essentially of balanc- 
ing a fraction of the 238 current against the 235 
current, and reading this fraction on a voltage 
divider. Readings were taken every ten minutes 
for a period of two to three hours, the relaxation 
time of the system. The resolution of the mass 
spectrometer remained good throughout the 
runs. The Beckmann thermometers showed tem- 
perature agreement to better than 0.01°C rise 
of one bulb over the other per hour. A variation 
of the order of ten times this could have been 
tolerated without affecting the accuracy of the 
measurement. 


At the end of the relaxation time the leak A 


was closed and the peaks in the mass spec- 
trometer dropped to their residual value. A 
standard sample was measured before and 
after the diffusion run. This permitted correction 
of the diffusion points for apparatus drift. This 


2A. O. Nier, E. P. Ney, and M. G. Inghram, “Null 
method for the comparison of two ion currents in a mass 
spectrometer,”” Chicago Meeting of the American Physical 
ety, June 22, 1946. 


correction was, on the average, about 4 percens 
of the total change measured in the diffusion ryp, 
Careful measurements of the pressure in the 
bulbs were made before and after the runs. 

Several methods were used to obtain the 
value of the concentration after complete mixing 
(C*). : 

In five of the runs C* was measured 20 hours 
after the start of the experiment by making an 
analysis of the material in the diffusion bulbs 
(i.e., comparing the 235/238 ratio with that of 
standard material); in two of the runs C® was 
obtained by pumping the material out of the 
bulbs into a copper trap and making a routine 
analysis of this; and in one of the runs C* was 
obtained by calculation from the known initial 
values of the concentration in bulbs F and B, 
Of course, in all cases C* could be thus calculated, 
and this was done as an added check. The agree- 
ment was found to be within the accuracy of 
measurement. 

In six of the runs the diaphragm in the sylphon 
gauge (G) was pressed against the inside plug of 
the gauge by admitting air above it. The dia- 
phragm thus acted as a valve closing off the main 
volume of the sylphon gauge from bulb F. In 
these instances the calculations were made on 
the basis of the two-bulb problem. In two of the 
runs the diaphragm of the sylphon gauge was 
forced against the upper stop of the gauge by 
pumping a vacuum through the valve Q. The 
volume of the sylphon gauge then communicated 
with the bulb F throughout the experiment. The 
calculations for this case was made on the basis 
of the three-bulb problem (described later). 


Il. THEORY 


The following is a list of symbols to be used in 
developing the theory: 


C;'=concentration in grams/cc of U**5F, in bulb B at 
time ‘‘?”’ 
C,'=concentration in grams/cc of U?*F, in bulb F at 
time ‘‘t” 
D=coefficient of diffusion in cm*/sec. 
V2=volume of bulb B in cc 
Vi=volume of bulb F in cc 
l=effective length of connecting tube 
A =effective area of connecting tube 
C” = equilibrium value of the concentration after complete 
mixing 
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The simplest case of the diffusion problem 





_ will be worked out first and the necessary cor- 


rections will be calculated later. 

The assumptions made in the simple theory 
are: 

(A) The concentration gradient in the con- 
necting tube is linear. 

(B) The concentration gradient is all in the 
connecting tube. 

(C) The volume of the connecting tube is neg- 
ligible compared with the volume of the bulbs. 

Under these conditions we can write the equa- 
tion for Q, the rate of transfer of U?*5F, by dif- 
fusion : 


Q=-—DA gradC. (1a) 
Because of assumptions A and B above: 
Q= —(DA/I)(C2'— Ci"). (1b) 


The rate of transfer in terms of the time is 
given by: 

Q= V26C,'/6t. (2) 
Equating (1) and (2), we obtain the differential 
equation 


5C2'/5t = —(DA/1V2)(C2'— Ci"). (3) 
In view of the relation (4) between C,' and 


C,', we can write Eq. (5) using 8 and y as pre- 
viously defined. 


V2C2'+ Vii! = VoC®, (4) 
5C,'/5t= —yC2'+B8. (5) 

We also have the boundary condition : 
C2*=C, at t=0. (Sa) 


The solution of (5) subject to the boundary 
condition is: 


C‘= C*+(C°— C”)e-7': (6) 


Equation (6) is the solution of the diffusion 
equation for the simple case. In order to apply 
it to our experiment it is necessary to modify it 
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by making proper corrections for assumptions 
A, B, and C. 

The assumption of a linear concentration 
gradient enters in the above as a substitution of 
A(C:—C,)/l for A-gradC. Since the concentra- 
tion gradient is linear in each element of length 
l; with cross section A;;1/A must be taken to 
mean }>°;1;/A; which is the equivalent single //A 
for the connecting tube. 

Since in our experiment, the //A for the bulbs 
was one-half of one percent of the //A for the 
connecting tube the assumption that the gradient 
is all in the connecting tube is justified providing 
a small end correction is made. For the purpose 
of making this correction, we assume that the 
conditions at the ends of the connecting tube 
are the same as they would be if the bulbs had 
infinite volumes. The magnitude of the correction 
for the analogous electrical case has been shown 
by Rayleigh*® to be an effective increase in the 
length of the tube at each end by about 0.82 
times the radius. This additional length is con- 
sidered only in >> ;/;/A,. It can also be shown by 
a method similar to that of Loschmidt‘ that the 
concentration gradient in the connecting tube is 
set up in a time of the order of 4 minute, with 
our geometry. 

The connecting tube has a volume that can 
not be neglected. Because of the symmetry of 
the connecting tube about the point of pinching 
the Neoprene connector, the average concen- 
tration in the connecting tube (Cz) is at any 
time the mean of the concentrations at the ends 
of the connecting tube. 


Cr =3(Ci'+Cr'). (7) 

Equation (4) must then be replaced by: 
V2Co'+ ViCi'+C ,Al=(Vo+Al)C*.. (8a) 
And using (7) and. grouping terms, (8a) becomes: 


(V2+Al/2)C2'+(Vi+Al/2)Ci' 
=(Vo+Al)C*. (8b) 


Equation (6) can be made consistent with this 
condition if one replaces in y, Vi by (Vi+Al/2), 
V2 by (V2+Al/2), and Vo by (Vo+Al). 


3 See Maxwell, Electricity and pee (Oxford, 1892), 
p. 434; Lord Rayleigh, Theory of Sound (The MacMillan 
Company, 1878), Vol. II, pp. 291-295. 

*Loeb, Kinetic Theory of Gases (McGraw-Hill Book 
Company, Ince, New York, 1934), pp. 268-272. 
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The quantity which embodies all of the above 
corrections, and which will replace y, will be 
called a. 


ts D | (Vot+Al) () 
Lily/ADL(Vi+-Al/2)(V2+Al/2)) 


The equation which describes the diffusion is 
then: 








Cc*—C;' 
=e~*, (10) 
Cc*—C,° 
The constant 1/a is the relaxation time of the 


system. 

As was mentioned earlier, in some cases the 
diffusion was allowed to proceed with the sylphon 
gauge open to the adjacent bulb. The solution 
for the three-bulb problem will be indicated here. 
We use the same notation which we used 
previously with the addition of: 


l,/A2=cffective 1/A between sylphon gauge and 
bulb F. 
V;=volume of sylphon gauge (G). 
C;'=concentration of U"*5F, in (G) at time “‘t.” 


The equations which must be satisfied are: 





















































Vi—=—(Ci'—C,'), (11a) 
bt l 
6C;' DA; DA 
r= (Cs*— Cy‘) +—(C2'—Ci'), (1 1b) 
bt l, l 
(Vit Vot V3)C*® — ViCit— V2C3' 
ye . (11) 
Vs 
° | 
ba? FIGURE 2. 
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Fic. 2. Graph of logio(C®— C')’as a function of the time. 
The slope of the least squares solution, together with the 
geometrical constants gives the diffusion coefficient. 
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With the boundary conditions: 
C.=C,°; (12a) 
Ci:=C;=C,°=C;°; when t=0. (12b) 


The solution of Eqs. (11) and (12) has for 
one of our numerical cases the form: 


Ci4—C* 


= 0.966 exp (—9.49DtX 10-5) 
C,°—C* 





+0.034 exp(—16.45DiX 10-5), 


If in Eqs. (11) C3'=C,', the equations reduce 
to the two-bulb case, for which the corresponding 
numerical solution is: 


C,'—C* 

C°—C* 
The value of D as calculated assuming a two- 
bulb problem, when the diffusion actually 


follows the three-bulb equations can be shown 
to be too small by approximately 0.8 percent. 





=exp(—9.72X10-*D2). 


Ill. CALCULATIONS AND RESULTS 


The time derivative of the logarithm of 
|C*—C,"| gives (from (10)): 


6 
— log.| C*—C,'| = —a. 
ét 


The slope of the plot of log.| C.*—C*| against 
time thus gives —a. A graph of logio(C*— C;') 
against time is shown in Fig. 2. The line indicated 
was obtained by the method of least squares. 
The C,' of the theory was obtained from the 
235/238 ratio measured in the mass spectrograph. 

The quantity D/a was calculated from the 
measured geometric constants (9). The diffusion 
constant was then obtained by multiplying this 
D/a by the experimental a obtained from the 
least-squares solution. 

In order to know pD (where p is the density), 
p was calculated assuming the ideal gas law to 
hold. This assumption has been shown to be 
justified at one centimeter pressure and room 
temperature. 

The results obtained are presented in Table I. 
As is indicated, half of the runs were made with 
C,' decreasing, and half with C;' increasing. 
Calling A the concentration difference between 
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TABLE I, Observational results. 

















Concen- OD ien* 
Cm Hg Two-bulb or C+ increasing tration 

Run Date made pressure Temp three-bulb or decreasing difference Geometry —— 
5 Nov. 10, 1942 1.08 cm 27.6°C III bulb decreasing A A 235 
6 Nov. 16, 1942 0.94 23.9° III increasing 4 A 215 
7 Nov. 25, 1942 1.08 27.4° II increasing A A 254 
8 Feb. 12, 1943 1.09 26.4° II decreasing 4A A 238 
9 Feb. 15, 1943 1.10 26.2° II increasing 4A A 234 
10 - Feb. 18, 1943 1.11 26.2° II decreasing 4A A 237 
11 Feb. 25, 1943 1,07 27.2° II decreasing 44 B 244 
12 March 1, 1943 1.08 26.9° II increasing 4A B 215 


Average (pD)s0°= 234 micropoises 








the bulbs at the start of runs 5, 6, and 7, the con- 
centration difference between the bulbs at the 
start of runs 8, 9, 10, 11, and 12 was about 4A. 
The average pD was the same for both of these 
conditions. This indicates that there is no 
memory effect in the mass spectrograph which 
would keep it from following the changes in 
concentration at the proper rate. Geometry A 
indicates that the Neoprene (D) was 2.08 cen- 
timeters long, geometry B means that the 
Neoprene was 6.10 centimeters. The values of 
pD corrected to 30°C are shown in the last 
column of Table I. The average value for (pD) 30° 
was 23443 micropoises. The cumulative error in 
the measurement of the geometrical constants is 
about 3 percent. The error in measuring the 
pressure is less than one percent. In view of this: 


(pD) s0° = 


Knowledge of pD makes it possible to deter- 
mine ¢ in the theoretical relation pD =en, where 
n is the viscosity coefficient. This ¢ is of peculiar 
interest in that a value for the force field ex- 


234+9 micropoises, 





* The value of (pD)m* was calculated from pD at the temperature measured (for example from (pD)m.s) by utilizing the known temperatur « 
dependence of viscosity. This was justified because of the smallness of the correction. 





ponent of the molecule may be deduced from it. 
According to classical theory « should range from 
1.200 for hard sphere molecules to 1.54 for mole- 
cules interacting as the inverse fifth power of 
their separation. 

The value of the ratio of pD to 7 at 30°C is 
1.31. This gives for the molecular force law 
F=const./r™. 
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The scattering cross sections of the proton, deuteron, and the carbon nucleus have been de- 
termined for 4.1-, 12.5-, and 13.5-Mev neutrons by detecting the corresponding recoil protons. 
A coincidence set of three proportional counters was used. Similar measurements were made 
with 14-Mev neutrons from the reactions D+Li. These were detected by the 9.9-min. activity 
arising from the reaction 29Cu®--on'—>29Cu™-+-2on'. The m—2n disintegration of the deuteron 
has been investigated for 14-Mev neutrons and cross section found to be of the order of 9-4 
X10-* cm*. The scattering cross sections at 4.1, 12.5, and 13.5 Mev are, respectively, for H 
1.73, 0.69, 0.694 10-* cm*; for D 1.79, 0.78, 0.864 10-* cm?; for C 1.99, 1.40, 1.23 10-* 
cm?, The experimental results are compared with those of other experimenters and with theo- 


retical predictions. 





1. INTRODUCTION 


NFORMATION about the nuclear forces can 
be deduced by comparing the experimentally 
determined energy-dependence of the neutron- 
proton scattering cross section with the corre- 
sponding theoretical predictions. The more valu- 
able experiments for such a purpose are those 
carried on with mono-energetic neutrons. Such 
neutrons are produced according to reaction 


1D?+ ,D?—,He?+ on'!, (1) 


whose spectrum energy is variable with both the 
observation angle and the voltage used to ac- 
celerate the impinging deuterons. 

Of the numerous experiments performed with 
reaction (1) we shall mention that of Aoki! who, 
by changing the observation angle, was able to 
determine the proton-neutron cross section for 
five different values of the energy of the im- 
pinging neutrons ranging between 2.14 and 2.76 
Mev. Other experiments performed with reac- 
tion (1) refer to neutron energies varying, from 
case to case, between 2.4 and 3 Mev.? Besides 
these measurements and those performed with 
slow neutrons,’ there have been others carried 

* Instituto di Fisica dell’ Universita, Centro di Fisica 
Nucleare del C.N.R. 

t Instituto Superiore di Sanita. 

1H. Aoki, Phys. Rev. 55, 795 (1939). 

2E. T. Booth and C. Hurst, Proc. Roy. Soc. Al61, 248 
(1937); R. Ladenburg and M. H. Kanner, Phys. Rev. 52, 
911 (1937); W. H. Zinn, S. Seely, and V. W. Cohen, Phys. 
Rev. 56, 260 (1939). 

3E. Amaldi and E. Fermi, Ricerca Scient. 7-1, 223, 
393 (1936); Phys. Rev. 50, 899 (1936); L. Simons, Phys. 
Rev. 55, 792 (1939); Kgl. Danske Vid. Sels. Math.-Fys. 
Medd. 16 (1940); V. W. Cohen, H. H. Goldsmith, and J. 


Schwinger, Phys. Rev. 55, 106 (1938); H. B. Hanstein, Phys. 
Rev. 57, 1045 (1940). 
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on with neutrons of about 0.2-0.3 Mev produced 
by nuclear photo-effect* or by the reaction® 


1D?+ 6C®—,N®+ on’. (2) 


For energies larger than 3 Mev there is a lack 
of published experimental data. As far as we 
know only Salant, Roberts, and Wang* and 
Salant and Ramsay’ have performed measure- 


ments with the neutrons from the reaction 
1D? + ;Li’—,Be®+ on! (3) 


detected by means of the activity of 9.9-min. 
period, produced in copper according to the 
reaction 


29 Cu® +- ot'—>99Cu™ +29, 


29 Cu®’—ogNi® + er. 


(4) 
(5) 


It is to be noticed that reaction (4) has an appre- 
ciable yield only for energies of the impinging 
neutron larger than 11 Mev,’ so that the results 
of these authors refer to neutrons of about 14 
Mev. : 

Finally, considering all the measurements per- 
formed up to the present, it is to be noticed 
that the experimental methods are different for 
different energy, excepting the measurements of 


4A. Leipunski, L. Rosenkewitsch, and D. Timoschuk, 
Phys. Sov. 10, 625 (1936); W. E. Good and G. Scharff- 
Goldhaber, Phys. Rev. 58, 89 (1940). 

5E. Amaldi, D. Bocciarelli, and G. C. Trabacchi, 
Ricerca Scient. 11, 121 (1940). 

* E. O. Salant, R. B. Roberts, and P. Wang, Phys. Rev. 
55, 984 (1939). 
( oun O. Salant and N. F. Ramsay, Phys. Rev. 57, 1075 
1 ’ 

8 R. Sagane, Phys. Rev. 53, 492 (1938). 
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Aoki, which cover a rather narrow energy 
interval. 

For a strict comparison with theory it would 
be very convenient to have measurements per- 
formed with mono-energetic neutrons whose en- 
ergy could be varied within wide limits. 

In Section 2 we shall explain a method which 
satisfies rather well such a condition ; it has been 
applied to the determination of the scattering 
cross sections of hydrogen, deuterium, and carbon 
for 4.1-, 12.5-, and 13.5-Mev neutrons. 

In other sections we attempt to clarify some 
interesting points about the neutron-proton and 
neutron-deuteron collisions.°® 


2. THE METHOD 


Our experimental arrangement is shown in 
Fig. 1. Neutrons emitted by the source S impinge 
on the paraffin layer P whose thickness (2 mm) 
is negligible with respect to the absorption of 
neutrons and large with respect to the range of 
recoil protons. Protons projected in the forward 
direction and fast enough to go through the 
Al-absorber A, are registered by a coincidence- 
set of three proportional counters, ¢1C2¢3. 

The triple coincidences are caused only in part 
by protons recoiling from P; in part they are 
caused by random coincidences whose number, 
depending on the type of neutron source and on 
the geometry of the arrangement can be de- 
termined by inserting in position O’ an absorber 
thick enough to stop all recoil protons. There- 
fore, all measurements were performed by count- 
ing the triple coincidences with this last absorber 
in position O (obturator opened) and then in 
position O’ (obturator closed) and taking the 
difference between the two numbers so obtained. 

The absorber A determines the minimum en- 
ergy Ey of the registered recoil protons. All the 
measurements were performed with protons scat- 
tered in the forward direction ; hence, the energy 
of recoil protons is equal to the energy of the im- 
pinging neutrons. If Emax is the maximum energy 
of the neutrons emitted by the source S, the 
difference between the numbers of triple coin- 
cidences counted with obturator open and closed 
is caused only by recoil protons from neutrons 
whose energy lies in the interval Ey— Emax. That 

*M. Ageno, E. Amaldi, D. Bocciarelli, and G. C. Tra- 


ROG Naturwiss. 31, 231 (1943); Nuovo Cimento 1, 253 
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is, our apparatus is a detector sensitive only to 
neutrons lying in a narrow energy interval close 
to the end of the spectrum of the neutrons from 
the source. This last limitation can be eliminated 
by replacing the thick obturator with another 
one chosen in such a way that only protons of 
energy equal or larger than E’ can be registered 
when both the absorber A and the new obturator 
are in the path of the particles recoiling from P. 

Under these last conditions the difference be- 
tween the numbers of triple coincidences counted 
with this new obturator open and closed is caused 
only by recoil protons from neutrons whose en- 
ergy lies between the two arbitrarily fixed values 
E, and E’. 

As we shall see in Section 4, our measurements 
of the cross sections have been performed, by 
means of the first of these methods, for neutrons 
of energy equal to 4.1, 12.5, and 13.5 Mev. 
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Fic. 1. Arrangement of source and scattering chamber. 
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With the second method a measurement at 
8 Mev was attempted; however, the result was 
rather uncertain on account of the large statisti- 
cal error affecting a rather small difference be- 
tween large numbers. 

In order to.determine the total (transmission) 
cross section of a given element for neutrons of 
the chosen energy interval, it is sufficient to 
insert, at about half-way between S and P, a 
scatterer D containing that element. The ratio of 
the numbers of recoil protons observed with and 
without such a scatterer in the position shown in 
Fig. 1, gives the transmitted fraction of imping- 
ing neutrons which is connected by a simple 
well known relation to the desired total cross 
section. 

In conclusion for every value of both the 
energy Ey (and Emax) and the thickness of the 
scatterer, we must perform four measurements 
under the following conditions: with obturator 
open and without scatterer: OO; with obturator 
open and with scatterer: OS; with obturator 
closed and without scatterer: CO; and with ob- 
turator closed and with scatterer: CS. 


3. THE EXPERIMENTAL ARRANGEMENT 


The experimental arrangement was about the 
same as that used in an earlier research; there- 
fore, details of the construction and the voltage 
(2140 v) applied to the counters and on the regis- 
tering set will not be given here. 

The scattering chamber was changed as shown 
in Fig. 1. The paraffin arc used in our preceding 
research was taken away from its support G and 
replaced by a 2-mm thick paraffin layer P of 
25.X40 mm? area, supported by a brass plate 1 
mm thick. 

The distance between the paraffin surface and 
the first wall of the third counter was 110 mm. 
The walls of the counters were of Al 1.93 mg/cm? 
thick. The scattering chamber was filled with 
commercial argon at atmospheric pressure. As 
shown in Fig. 1 the distances between the source 
S and scatterer D(d,;) and between scatterer 
D and paraffin P(d.) were both 200 mm. 

The neutron source S was a target of Be, B, 
or .Li bombarded with deuterons of 0.8 Mev pro- 


10E. Amaldi, D. Bocciarelli, B. Ferretti, and G. C. 
Trabacchi, Naturwiss. 30, 582 (1942); Ricerca Scient. 13, 
502 (1942). 
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TaBLE I. Data on scattering as determined by 
proton. Each transmission determination is the resy 
200 measurements except as indicated. 


y rool 
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Thickness + ee* 4 recoil 
(cm) (g/em*) OO © rotors 
D+Be; Eo =3.2 Mev; Emax =4.6 Mev; E=4.1 Mev 
6601 4802 732 659 0.706 
6055 3381 727 538 0.534 rey +r 
6403 2235 791 454 0.317 +0.011 


5836 3329 698 475 0.555+0,015 
6305 2602 743 «545 «0.370 +.0.012 


Cc . . 6010 3336 702 526 529 +0.014 
le 6345 2808 720 541 0.403+0.012 
8.0 J 6208 2191 753 S515 0.30740.011 
D+B; Eo=9.4 Mev; Emax =13.2 Mev; E =12.5 Mev 
6.0 5.2 2815 1683 923 655 0.543+0.022 
485 4.74 2647 1795 889 692 0.627+0.024 


4.0 6.35 2843 2036 853 753 0.64440.024 


Scat- 
terer 





CuHes 


CuDn 








D+Li; Eo=11.5 Mev; Emax =14.5 Mev; E =13.5 Mev 

13706 10382 
2 13094 

8 9462 
94 14063 


5343 
5369 
3664 


5096 
5462 


4248 
5136 
3018 


4128 
3318 
1773 


3970 
3512 


3680 
3469 


0.748 +0.019 
0.518 +0.0161 
0.445 +0.0162 


0.702 +0.017 
0.600 +0.018 


0.817 +0,0203 
0.581 +0.016 
0.385 +0.015¢ 


295 2.6 
; 7322 

4354 

10269 
74 13027 8054 
.3.17 12211 
5 13560 

9 8461 


10183 
8366 
3350 








1 Result of 160 measurements. 
2 Result of 150 measurements. 
3 Result of 210 measurements. 
4 Result of 130 measurements. 


duced by the one million-volt equipment of the 
Istituto Superiore di Sanita. The target was 
reduced to a disk of 20-mm diameter by means 
of an iron screen. The neutrons were observed 
at 90° with respect to the direction of the im- 
pinging deuterons. 

The scatterers had an area equal to 30X30 
mm? and thicknesses varying between 20 and 90 
mm. 

In order to count the triple coincidences al- 
ternatively in the four conditions OO, OS, CO, 
CS, we have used the same automatic device 
used in our preceding research, conveniently 
modified ®: the scattering chamber containing 
the paraffin P and the counters remained in a 
fixed position while the scatterer D, supported 
by a long 2 mm thick aluminum arm, was peri- 


' odically inserted and taken away from the path 


of the neutrons coming from S and impinging 
on P, 

The triple coincidences were registered for 28 
sec. in each of the four above-mentioned arrange- 
ments; then the arrangement was changed auto- 
matically in 2 sec. By alternating so frequently 
the four measurements, the effect of slow varia- 














tions of the neutron emission or detection sensi- 
tivity was made as small as possible. By means of 
a convenient caliper the distances and the align- 
ment were frequently checked during the 


measurements. 
4. THE MEASUREMENTS 


The measurements were performed with neu- 
trons from the three reactions Be(D, n), B(D, n), 
Li(D, m). The more critical point in all measure- 
ments of the transmitted fraction of neutrons is 
to ascertain that the number of neutrons scat- 
tered by the surrounding materials and room- 
walls is negligible with respect to the number of 
neutrons impinging on the detector directly from 
the neutron source. 

Although the experiments described in sub- 
section (6a3) and (6a4) of our earlier paper!® had 
given very satisfactory results on this point, we 
performed two more checks with D+Li neutrons 
at Ep=11.15 Mev. 


(a) We have verified that the number of registered pro- 
tons varies proportionally to the square of the distance 
r=d,+d; of the paraffin P from the source S. By increasing 
this distance from r’ = 16.6 cm to r’”’ = 29.6 cm, the number 
of registered protons decreased by a factor 3.12+0.15, in 
satisfactory agreement with (r’’/r’)? = 3.18. 

(b) In our apparatus the target was cooled by means of 
a water layer about 2 mm thick which could give rise to 
a few scattered neutrons. In order to check this point we 
substituted for the scattered D a 6-cm thick paraffin ring 
whose inner and outer diameter were, respectively, 4.5 
and 9.0 cm. We found that the paraffin ring increased the 
number of registered protons by a factor 1.109+0.031. 
Considering the thickness and the geometrical disposition 
of the water-cooling, we can conclude that it gives rise to 
a negligible number of neutrons scattered to the detector 
with energy high enough to be registered. 


We measured the fraction of neutrons trans- 
mitted by scatterers of light paraffin, heavy 
paraffin, and graphite. The melting point of the 
light paraffin was between 68° and 70°C so that 
we can assume the formula Cs2H¢s. The melting 
point of heavy paraffin was higher than 77°C so 
that we can assume the formula CssD7; the 
graphite used was ‘‘Acheson.” 

In the measurements with neutrons from Be, 
B, and Li the Al absorber A was, respectively, 
0 mg/cm?, 124 mg/cm?, and 175.5 mg/cm? thick. 
Taking into account the thickness of the five 
counter-walls (5X 1.93 =9.6 mg/cm? Al), the dis- 
tance (11 cm of argon) between the paraffin P 
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and the first wall of the third counter, and the 
stopping-power of aluminum," we found that the 
minimum ranges of recoil protons, in the three . 
above-mentioned conditions, were, respectively, 
17.5, 103, 140 air cm, corresponding to Eo =3.2, 
9.4, 11.15 Mev. 

The results of our measurements are given in 
Table I; in addition to the minimum energy Eo, 
the maximum energy Emax, and the mean en- 
ergy E are also given. For the neutrons from B 
and Li Ensx was deduced from the absorption 
curves of recoil protons given in our earlier 
paper.’ A similar curve was measured in the 
case of the neutrons from Be which turned out 
to agree very well with the results of other 
authors.” i 

The mean energy E was calculated by formyla 
(11) of our earlier paper from the above-mentionyd 
absorption curves for B and Li neutrons, and the 
spectrum of Bonner and Brubaker for Be ngu- 
trons. In Table I, besides the thickness of ;he 
scatterers, we give the numbers of recoil protijns 
registered in the 4 dispositions OO OS CO CS, 4ad 
the transmitted fraction 5 


T =(OS—CS)/(00—CO). 46) 


The errors are calculated from the square yoot 
of the number of counted pulses, by means of the 
usual propagation rule. In some cases we have 
calculated the standard deviation which turped 
out to be about 1.3 times the error calculated 
from the square root of the number of pulses. 

In the case of the neutrons produced by means 


of reaction 
,D?+,B* ¥—>,C™ 124 om) (7) 


the measurements were performed for a single 
value of the thickness of the scatterer because 
the mean energy of these neutrons does not differ 
very much from that of neutrons produced in 
reaction (3) which has the advantage of a much 
higher intensity. 

We have also attempted to measure the trans- 
mitted fraction of the neutrons of the second 
group of reaction (7) which extends between 7.3 
and 9.1 Mev. For such a purpose the obturator 
was replaced by an Al absorber 60.6 mg/cm? 


1M. S. Livingstone and H. A. Bethe, Rev. Mod. Phys. 
9, 276 (1937), 

%T. W. Bonner and W. M. Brubaker, Phys. Rev. 50, 
308 (1936). 











24 AGENO, AMALDI, BOCCIARELLI, AND TRABACCHI 


thick; the Al absorber A was 54.8 mg/cm? thick; 
and the paraffin P 43 mg/cm’, corresponding to 
43 cm air. Under these conditions with obturator 
open and closed we counted only recoil protons 
whose range was equal or larger than, respect- 
ively, 66 cm of air (7.3 Mev) and 97 cm of air 
(9.1 Mev); the calculated mean energy E was 
8 Mev. 

From these measurements we found, for each 
scatterer studied, a result intermediate to that 
obtained at F =4.1 Mev (D+Be) and at E =13.5 
Mev (D+Li). However, the statistical errors 
were appreciably higher than that reported in 
Table I, as we must expect, taking into account 
that in this case we deal with relatively small 
differences between rather large numbers. There- 
fore we do not report in detail these results which, 
although in agreement with all other measure- 
ments, are not usable for a strict comparison 
with theory. 

In order to check our results by a different 
method we repeated the measurements of the 
transmitted fraction of the D+Li neutrons, us- 
ing, as detector, the activity induced in copper 
according to reaction (4). This method *’ is 
very simple; it is, however, to be noticed that 
the value of the threshold of reaction (4) is 
rather uncertain and the dependence of its cross 
section on the energy of the impinging neutrons 
is practically unknown. For this reason the agree- 
ment between the results obtained by detecting 
fast neutrons by means of recoil protons and by 
means of the copper activity, must be considered 
as a reciprocal check of the two methods. 

In order to have geometrical conditions simi- 


TABLE II. Scattering of neutrons from the D+Li re- 
action detected by use of the Cu® 9.92-min. activity. 











Num- Num- 
ber of ber of 
counts counts 
without with 
scat- Zero scat- Zero Transmission 
cm g/cm? terer effect terer effect al T 
Scatterer CsaHes 


2.95 2.6 28016 4820 21696 4970  0.721+0.009 0.748+0.019 
6.0 5.2 25829 4900 16848 4857  0.573+0.009 0.518+0.016 
8.95 7.8 29768 5170 15648 5060  0.430+0.008 0.445+0.016 


Scatterer CasDrn 


3.0 2.94 29502 5530 23190 5530  0.730+0.009 0.702+0.017 
485 4.74 27597 5150 18640 5160  0.592+0.008 0.600+0.018 


Scatterer C 


2.00 3.17 25044 5080 22576 4980  0.843+0.011 0.817+0.020 
6.0 9.5 27746 5000 17712 4870  0.564+0.008 0.581+0.016 








10.0 159 28202 4980 13840 4910  0.384+0.007 0.385+0.015 











lar to those used in the preceding measurements, 
the copper detectors were plates, 1 mm thick, 
bent in the form of cylinders 25 mm in diameter 
and 40 mm long so that their section (25x49 
mm”) was equal to the area of the paraffin sheet 
of Fig. 1. The distance between source and de- 
tector was again 40 cm and the scatterer was 
placed at half this distance. 

To reduce, as much as possible, the number of 
scattered neutrons, all surrounding materials 
were removed; the scatterer and the detector 
were supported by an aluminum plate, 1 mm 
thick and 3 cm wide, marked with convenient 
lines, which permitted the distances and the 
alignment of scatterer and detector in the suc- 
cessive experiments to be reproduced accurately, 
During irradiation the detector was screened 
with cadmium 0.78 g/cm? thick to eliminate the 
copper activity of 5-min. period produced by 
slow neutrons. 

The influence of slow intensity-variation in the 
neutron emission was eliminated by irradiating 
at the same time two copper detectors placed at 
the same distance from the source, in two direc- 
tions forming an angle of about 30°: the scatterer 
was placed in the path of the neutrons impinging 
on one of these two copper cylinders. 

After ten minutes irradiation, the two detec- 
tors were placed around two thin-walled glass 
counters and the number of pulses was registered 
by means of a “‘scale-of-16” and a “‘scale-of-8.” 
The counting continued for 19 minutes, starting 
one minute after the end of irradiation. The 
measurement was repeated with the same scat- 
terer placed in the path of the neutrons impinging 
on the other copper cylinder and the results of 
the two measurements were added together. The 
zero effect was measured before and after each 
measurement. 

The results are summarized in Table II; the 
values given are the sum of the pulses counted 
in four measurements. In column 8 we give the 
transmitted fraction J’ deduced from the pre- 
ceding data, while in column 9 we reproduce, for 
comparison, the transmitted fraction of neutrons 
T deduced by means of (6) from the recoil- 
proton measurements (Table I, D+Li). 

The agreement is very satisfactory so that it 
seems to us reasonable to assume that the meas- 
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TABLE III. Scattering cross sections from the data of Tables I and II. 














Molecular 






























« X10" « X10" 
Scatterer a x’ kK” T-K o X10" (average) (average) 
D+Be £=4.1 Mev 
Cc 6.34 0.003 0.003 0.523 2.04+ 0.08 
A=6.5 cm 9.51 0.004 0.007 0.392 1.994 0.07 o(C)=1.99 +0.04 
protons 12.68 0.005 0.012 0.290 1.97+ 0.06 
CD 2.94 0.003 0.003 0.549 192+ 8.7 
ane om 4.74 0.004 0.008 0.358 204+ 7.0 1995.4 o(D)=1.79 20.08 
protons 
Ca2Hes 1.4 0.003 0.000 0.703 188+ 12 
A=4.9 cm 2.6 0.006 0.001 0.527 185+ 7.6 178+4 o(H)=1.73 +0.06 
protons 5.2 0.008 0.006 0.303 173% 5.2 
D+B E£=12.5 Mev 
Cc 6.35 0.003 0.001 0.640 1.40+ 0.12 o(C)=1.40 +0.12 
protons 
CasDre 4.74 0.003 0.002 0.622 95+ 7.7 o(D)=0.78 +0.12 
protons 
CszHes 5.2 0.006 0.002, 0.535 90.54 6 o(H)=0.69 +0.11 
protons 




















#=13.5 Mev 3.17 0.001 0.000 0.816 
os 9.5 0.003 0.002 0.576 
A=10.6 cm 15.9 0.004 0.008 0.373 
protons 
f=14 Mev 3.17 0.001 0.000 0.842 
S 9.5 0.003 0.002 0.559 
Cu 15.9 0.004 0.008 0.372 
B=13.5 Mev 
CssDr2 2.94 0.002 0.001 0.699 
A=9.2 cm 4.74 0.003 0.002 0.595 
protons 
B=14 Mev 2.94 0.002 0.001 0.727 
CssD72 4.74 0.003 0.002 0.587 
u 
£=13.5 Mev 2.6 0.003 0.000 0.745 
CuHes 5.2 0.006 0.002 0.510 
A=9.4 cm 7.8 0.007 0.004 0.434 
protons 
E=14 Mev 2.6 0.003 0.000 0.718 
CaHee 5.2 0.006 0.002 0.566 
Cu 7.8 0.007 0.004 0.419 








D+Li #=13.5 Mev 


1.29% 0.16 
1.17+ 0.06 
1.254 0.0 


1.22+0.037 


on 


o(C)=1.23 +0.015 


1.09+ 0.08 
1.234 0.03 
1.244 0.02 


1.23+0.016 


115+10 


o(D) =0.864+0.028 
105+2.2 


86.5+2.6 


+ 

+ 3.5 

+ 2.3 5.0+1.35 
+ 1.9 


SY OO 


(H) =0.694+0.019 












urements, performed with both methods, cannot 
be affected by serious systematical errors. 

As a further check of the measurements per- 
formed with a copper detector, we have plotted, 
on a semilogarithmic scale, the decay curve of 
Cu obtained by adding the pulses counted in 
each minute, in all the measurements used to 
deduce Table II. Such a decay curve turned out 
to be a straight line corresponding to a half pe- 
riod of 8.30.4 min. As this value is appreciably 
lower than that obtained by Crittenden™ (9.92 
+0.05), one could suspect a slight influence, in 










8 E. C. Crittenden, Phys. Rev. 56, 709 (1939). 








our measurements, of the 5-min. period. Although 
our decay curves can be analyzed only with 
difficulty into two exponentials of 5- and 9.9- 
min. periods, we have calculated under this 
assumption the percentage of the total number of 
pulses which could be attributed to slow neu- 
trons. We found less than five percent so that 
we considered it unnecessary to introduce any 
correction for this reason. 


5. RESULTS 


If the transverse dimensions of neutron-source 
S, detector P and scatterer D are negligible with 
respect to the distances d,, dz, the percentage 
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of the transmitted neutrons varies exponentially 
with the thickness z of the scatterer 


en alA, (8) 


A=1/Nce is the mean free path, ¢ the total cross 
section, and N the number of nuclei per cm’. 

If 6 and M are the density and the atomic 
weight of the scatterer we have 


2 zéo 6 o 0) 
A M1.67X10-% M1.67X10-%4 


where A=2z6 is the thickness of the scatterer in 
g/cm’. 

In our experimental disposition the transverse 
dimensions of the neutron source and detector 
were small but not negligible with respect to 
d; and d;. Therefore we must introduce two cor- 
rections caused by: 





(1) “scattering in,” that is, neutrons which are scattered 
through such a small angle that they still impinge on the 
detector. 

(2) “double scattering,” that is, neutrons which impinge 
on the detector on account of a second collision in the 
scatterer. 


Both these effects increase the percentage T of 
the transmitted neutrons with respect to the 
value (8) one would expect under ideal geo- 
metrical conditions. Therefore we can write 


T(z) =e-/4+-K'+K”", (10) 


where K’ and K” correspond to the above-men- 
tioned effects. 

By simple geometrical considerations we found 
for our experimental conditions 


a(0) 
K’ =0.007——-(1 —e-*/4) (11) 
é 


for scatterer of atomic weight larger than 1, and 


a(0) 
K’=0.014——(1—e-*) (12) 


for hydrogen. o(0) is the differential cross section 
in the forward direction and ¢ is the differential 
cross section averaged over all directions. If the 
scattering is isotropic in the center of gravity 
system we have o(0)/¢=1. 

In practice, we have used for D and C, 


K'=0.007(1-T); (11’) 
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Fic. 2. Absorption curves of recoiling particles. 


for H 
K’=0.014(1—T) (12’) 
in the case of neutrons of 4 Mev, and 
K’=0.018(1—T) (12”) 


in the case of neutrons of 12.5-13.5 Mev as 


suggested by our preceding results." 

K” was calculated numerically in the case of a 
scatterer of 30X30 mm? section and z=A=90 
mm thick. We found, for nuclei of atomic weight 
larger than 1, 

K"’ =0.008; (13) 
for hydrogen 
K"” =0.002. (14) 


The low value of K” in this case is caused by 
the energy-loss that a neutron undergoes in colli- 
sions against hydrogen. As we must expect that 
K" varies roughly proportionally to (z/A)? we used 


K"” =0.008(z/A)? (13’) 
in the case of D and C, and 
K"” =0.002(z/A)? (14’) 


in the case of hydrogen. 
Of course in the case of paraffin we used the 
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average of the corrections due to carbon and 
hydrogen. 

Table III contains the cross sections calculated 
by means of (9) and (10); we give also the 
values of K’, K’, and T—(K’+XK”) in order to 
show that the corrections represent always only 
a few percent of the transmitted fraction T of 
neutrons. The mean free path given in column 
1 is obtained from a semilogarithmic graph of the 
uncorrected transmitted percentage of neutrons; 
therefore it is only an approximate value. In 
the first column we have also indicated the de- 
tector: recoil-protons or Cu-activity. In this last 
case, according to Salant and Ramsay,’ we as- 
sume for the mean energy of the detected neu- 
trons E=14 Mev. Such a high value is justified 
by the presumably rapid increase of the cross 
section of process (4) with the energy of the im- 
pinging neutrons. 

Column 6 contains the cross section per mole- 
cule whose weighted averages are given in column 
7. The atomic cross sections are given in the fol- 
lowing column. Finally for D+Li neutrons we 
have averaged the values obtained with two 
different detection methods although the cor- 
responding mean energies are slightly different. 

Such a procedure appears justified considering 
that the two averaged values are equal to one 
another within the experimental errors, and that, 
under all known reasonable assumptions about 
the proton-neutron interaction, the scattering 
cross section does not change appreciably in a 
narrow energy interval. 


6. COMPARISON OF THE NEUTRON-PROTON AND 
NEUTRON-DEUTERON COLLISION 

Although the results contained in the last 
column of Table III are to be discussed in Sec- 
tion 7, we will now compare the neutron-proton 
and neutron-deuteron cross sections. From Table 
III we see that they are equal, within the ex- 
perimental errors, at 4.1 Mev, while at higher 
energies the deuteron cross section is slightly 
larger than the proton cross section. From the 
data of Table III we have, at 14 Mev, 


op/on =1.22+0.05. (15) 


For such a high energy the neutron-proton cross 
section is caused only by the scattering, while 
in the case of deuterons we must consider also 
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the (m, 2m) process 
on! + ,D?—>,H!+2on'+Q, (16) 


where Q is equal to the mass-defect of the deu- 
teron, i.e., Q@=—2.17 Mev. Therefore we can 
write 

op =0.+¢4, (17) 


where ¢, is the cross section for elastic scattering 
and oq the cross section of process (16). 

To investigate reaction (16) and to evaluate 
the corresponding cross section, we have per- 
formed the following experiments: 

(1) We replaced the light paraffin P with a 
heavy paraffin layer and measured the absorp- 
tion curve of the recoiling particles. This curve 
is reproduced, as a function of the range in cm 
of air, in Fig. 2 (curve D); in the same figure we 
give also the absorption curve of protons re- 
coiling from light-paraffin (curve H). 

As we have shown in an earlier paper,” the ordi- 
nate N(Eo) of the two curves of Fig. 2, are given 
by the equation 


Nog ts) C) 
N(E)=— J f f(E)o(E)dE, (18) 


o “E 


where F=dE/dR is the energy loss of the re- 
coiling particles in light or heavy paraffin, No 
the number of protons or deuterons per cm’, 
f(£)dE the number of neutrons emitted by the 
source with energy between E and E+dE, o(E) 
the corresponding cross section, and g a pure 
geometrical factor. Under the reasonable assump- 
tion that o(Z) varies only slowly in the energy 
interval involved, and considering that No was 
about the same in the light (7.4410") and 
heavy paraffin (7.35 X10”), we can write 


1 
N(Eo) eee), (19) 


where 
@ 


#(E,) = f f(E)dE (20) 
E 


Ee 


is the twice integrated spectrum of the emitted 
neutrons and é@ the cross section averaged over 
the energy interval extending from E,» to the 
maximum energy Emax of the neutron spectrum. 

(2) The energy of a deuteron recoiling in a 
head-on collision with a 15-Mev neutron is 
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8/9 of 15 or 13.3 Mev corresponding to a range 
in air Rp=112 cm. Considering that Ens: for 
the D+Li neutrons is certainly lower than 
15 Mev, curve D of Fig. 2 gives evidence that 
the paraffin also emits a few protons. In order to 
settle this point, we took the paraffin layer away 
from its position P in Fig. 1 and arranged on the 
two opposite faces of a brass plate 2 mm thick 
two layers of light and heavy paraffin each 1 
mm thick. The brass plate was secured to the 
pivot G so that we could change very rapidly 
the light and heavy paraffin and compare the 
numbers of recoiling particles. Using an Al ab- 
sorber A so thick that the registered particles 
had a range equal or larger to 116 cm of air, we 
found 
Ni(116)/N(116) = 14+2.2, (21) 


i.e., the heavy paraffin bombarded with D+Li 
neutrons emits protons of energy larger than 10 
Mev whose number is equal to about 7 percent 
of the number of protons of the same energy 
recoiling from light paraffin. In Fig. 2 the two 
circles plotted at R= 116 cm represent the result 
of this measurement. 

(3) The number of protons emitted from the 
heavy paraffin we used is so small that it could 
be caused completely or in part by light hy- 
drogen impurity. Assuming again that the maxi- 
mum energy of the D+Li neutrons is 15 Mev, 
the maximum energy of protons emitted in the 
forward direction in reaction (16) is equal to 
12.5 Mev corresponding to 171 cm of air. There- 
fore we repeated a measurement similar to that 
described in (2) for this higher value of the range. 
We found 

N,(171)/N.(171) =11+6. (22) 


The two circles marked in Fig. 2 at R=171 cm 
represent this result. 

Although the experimental error is very high, 
we can conclude that the observed protons are 
due at least in part to hydrogen impurity. 

The heavy paraffin is supposed to contain no 
more than five percent of light paraffin,* which 
is consistent with our preceding results. 

In any case, from the above results it is im- 
possible to form any conclusion about the exist- 


*The Norsk Hydro-Elektrisk Kvoelstofaktieselskab, 
from which we bought the heavy paraffin, kindly com- 
municated these data to us. ' 
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ence of process (16) except that its cross section 
must be very small for neutrons of 13-14 Mey. 

(4) Finally instead of comparing the numbers 
of particles, emitted from light and heavy paraf. 
fin, having the same range, we have compared 
the numbers of recoiling particles due to ney- 
trons of the same energy. The experiment was 
performed for neutrons of energy equal or larger 
to 13 Mev which give rise to recoiling protons 
and deuterons of, respectively, 13 and 11.15 Mey. 

The corresponding ranges in air (R,=183, 
Rp =88 cm) were obtained by using two conveni- 
ent Al absorbers in the two cases (223 and 87.6 
mg/cm?). We remember that in this experiment 
the distance from paraffin to the first wall of the 
third counter was 20 cm (of argon). We found, in 
72 measurements, 


Nw =57654+134 N,=Ny=2724+122 
Nw /N,=2.12+0.11. 


These data, represented in Fig. 2 by two crosses, 
were used to adjust the two absorption curves 
on the same scale. 

Now we must consider that the particles 
emitted from the heavy paraffin are of three 
types: recoil deuterons, disintegration protons, 
and recoil-protons due to light paraffin impurity. 
The number of these last can be deduced by 
curve H of Fig. 2 once we know the percentage 
(five or’ seven percent) of the impurity: for 
instance the lower curve of Fig. 2 was obtained 
by multiplying curve H by 0.05. 

In the following we shall show some simple 
calculations in both cases of five and seven 
percent impurity. 

We have 


5 
5% Np=Ne—-——Nx=3665, (23) 
~ Ms 100 

7 
1% Np=Ne—-—Nu=2865.  (23’) 

100 


Neglecting a constant factor and applying (19) 
to our case, we have 


1 1 
Np =—<o. (13) +—o¢8(10.9), 
Fp Fy 


(24) 





1 
Na =—son®P(13). 
Fu 





Intr 
(26) 


in t 
valu 
of tl 
Rec: 

















The energy of the neutrons which gives rise to 
disintegration protons in the forward direction 
fast enough to be registered (8.6 Mev corre- 
sponding to R,=88 cm) is 10.9 Mev. From (24) 


we deduce 
Np Fu o.(13) -, ga(10.9) (10.9) 


ae oe i . (25) 
Na Fp cn(13) on(13) (13) 





From the range-energy curve" we have: for 
protons of energy between 13 and 14 Mev, 
Fa=3.84 10-*? Mev/cm of air; for deuterons of 
energy between 11.55 and 12.5 Mev, Fp=7.5 
10-2 Mev/cm of air: 


Fy/Fp=0.5. (26) 


The factor ©(10.9)/#(13) can be obtained from 
curve H of Fig. 2 


®(10.9 Mev) _@(134cm) 
6(13 Mev) (183cm) _ 





(27) 


Introducing the numerical values (23), (23), 
(26), and (27) in (25) we find 


a(13) _oa(10.9) 
6 








5 1.35=0.5 . 28 
70 on(13) on(13) (28) 
o(13)  oa(10.9) 
7 1.05 =0.5 .  (28' 
% on(13) on(13) (8) 


We remember that the cross sections appearing 
in these relations are averaged from the given 
values of the energy up to the maximum energy 
of the neutrons of D+Li reaction (14.5 Mev). 
Recalling (15) and (16) we can write 


o./on =1.22x, op/oH=1.22(1—x) (29) 


where 
x=0./(¢6-+¢a) (30) 


is the percentage of neutron-deuteron collisions 
which gives rise to recoil deuterons. 

Under the reasonable assumption that o_ does 
not vary very much when the energy of the im- 
pinging neutrons changes from 11 to 14.5 Mev, 
we can introduce the relations (29) in (28) and 
(28’) and calculate the value of x. We find 


5% x=(6—1.1)/(6—0.5)=0.89, (31) 
71% x=(6—0.9)/(6—0.5)=0.93. (31’) 
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The result is not very sensitive to errors of 
the ratio ©(10.9)/#(13) or to the percentage of 
light paraffin. In order to have x =1 it would be 
necessary to have Np/Ny=1.22X0.5=0.6, ie., 
a percentage of light paraffin equal to 13-14 
percent, which seems definitively too high. 


7. DISCUSSION 


Excepting for the results of the measurements 
at 12.5 Mev (D+B), which are not very ac- 
curate for the reasons explained in Section 4, all 
other values of the cross sections are affected by 
errors of the order of a few percent. As we have 
already pointed out in Section 5, the good agree- 
ment between the results obtained at 14 Mev 
using two detectors of quite different type, can be 
considered as a satisfactory reciprocal check of 
the applied methods. By averaging these two 
groups of results, the experimental errors were 
correspondingly lowered. We can now compare 
our results with that of other authors and with 
the theoretical predictions. 

As far as we know the neutron-proton and the 
neutron-deuteron cross sections were not yet 
measured for neutrons of about 4 Mev. Going to 
the theoretical predictions, we can use, for the 
neutron-proton collision, the formula deduced by 
Kittel and Breit,“ considering the scattering of 
the S wave by a square potential-well of the 
Heisenberg-Majorana type. It differs from pre- 
ceding simpler formula only because the finite 
value of the radius of nuclear forces was con- 
sidered and assumed equal to the electronic 
radius, as suggested by the proton-proton scat- 
tering experiments. 

The scattering of the P wave seems to be 
negligible at such a low neutron energy under 
all usual assumptions about the nuclear forces. 
From that formula we find, at 4.1 Mev, 


o=1.95X10- cm’, 


i.e., a value appreciably higher than our experi- 
mental result (1.73+0.0610~-* cm’); we can- 
not, however, exclude the possibility that such a 
difference is mostly caused by experimental 
errors. 

The neutron-deuteron collision was investi- 
gated theoretically by Hécker™ who has calcu- 


“4 C. Kittel and G. Breit, Phys. Rev. 56, 744 (1939). 
16 H. Hicker, Physik. Zeits. 43, 236 (1942). 
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lated the cross sections for elastic scattering, 
disintegration according to reaction (16), and 
radiative capture. 

He used, for the interaction between the three 
particles, the Voltz-potential, i.e., a linear com- 
bination of Majorana, Heisenberg, Bartlett, and 
Wigner forces decreasing with the distance ac- 
cording to a gaussian curve. The relative values 
of these four terms were adjusted by Voltz to 
obtain the well-known saturation effect of the 
forces and the experimental values of the deu- 
teron constants. 

At 4.1 Mev we need to consider only the elastic 
scattering whose cross section was, unfortun- 
ately, calculated by Hécker only up to about 3 
Mev where it is equal to 1.5<10-* cm?. How- 
ever, considering the regularly decreasing shape 
of the Hécker curve, we can conclude that our 
experimental value is appreciably higher than 
the theoretical one. 

The same conclusion can be drawn from the 
experimental results of other authors: at 2.4 
Mev Kikuchi and Aoki” found 1.98 10-** cm? 
where the Hécker value is 1.60 X 10-* cm?; Zinn, 
Seely, and Cohen" found 2.17+0.08 x 10-** cm? 
at 2.88 Mev; the value given by Hécker is 
1.55 X 10-* cm?. 

Also for thermal energy the experimental 
values obtained by various authors, and par- 
ticularly by Dunning and Carroll,'!® are higher 
than the theoretical one. 

Considering that in neutron transmission ex- 
periments all systematical errors tend to give 
cross sections lower than the right values, it 
seems to us reasonable to attribute the above 
stated discrepancy to defects of the theory. It 
is however to be noticed that the above men- 
tioned experimental data differ even more from 
preceding calculations of other authors.”° 

The results obtained at 12.5 Mev (D+B) al- 
though in agreement with our other measure- 
ments, are not accurate enough for a strict com- 
parison with theoretical predictions. 


16H. Voltz, Zeits. f. Physik 105, 537 (1937). 

17S. Kikuchi and H. Aoki, Proc. Phys. Mat. Soc. Japan 
21, 75 (1939). 

18 W. H. Zinn, S. Seely, and V. W. Cohen, Phys. Rev. 
56, 260 (1939). 
( 19See L. Motz and J. Schwinger, Phys. Rev. 58, 26 
1940). 
21. I. Schiff, Phys. Rev. 52, 242 (1937); S. Fligge, 
Zeits. f. Physik 108, 545 (1938). 








For neutron energies of about 14 Mey the 
neutron-carbon and neutron-proton cross sec. 
tions have been measured by Salant and Ram. 
say’ by detecting the D+Li neutrons by means 
of the copper activity produced according to 
reaction (4). These authors performed measure- 
ments both at 90° (Z=14 Mev) and 0° (£=15 
Mev) with respect to the direction of the im. 
pinging deuterons on lithium. The results ob- 
tained in the first group of measurements can 
be compared with our results: they found 


o-=(1.27+0.04) x 10-4, 
on = (0.70+0.06) X 10-4 cm? 


in very good agreement with our results 


o-= (1.230.015) X10, 
on = (0.694+0.019) X 10-* cm, 


According to the Kittel-Briet formula the proton 
cross section at 14 Mev is 0.66 X 10- cm’. 

The contribution to the total cross section due 
to the scattering of the P wave wascalculated 
by Kittel and Breit'* and by Ferretti* for neu- 
trons of 16 Mev, under the assumption that the 
nuclear forces are of the neutral type studied by 
Bethe.” According to these authors the P wave 
contribution represents nine percent of the cross 
section deduced by the scattering of the S wave. 
With the neutral meson forces studied by Rarita, 
Schwinger, and Nye* the P-wave contribution 
would represent 46 percent of the cross section 
deduced by the S-wave scattering. 

As we have already pointed out in an earlier 
paper,!® the small difference between the experi- 
mental value and the scattering cross section 
deduced: by the Kittel-Breit formula is in better 
agreement with the neutral meson theory of 
Bethe than with the neutral meson theory of 
Rarita and others. 

Finally we consider the neutron-deuteron colli- 
sion experiments at 14 Mev. As we have shown 
in section 6, for neutron energy of 13-15 Mev, the 
cross section of the (m, 2m) deuteron disintegra- 
tion process is equal to about 10 percent of the 
total collision cross section. 

In this connection it is to be noticed that, 


21 B. Ferretti, Ricerca Scient. 12, 843, 993 (1941). 

2H. A. Bethe, Phys. Rev. 57, 260, 390 (1940). 

% W. Rarita, J. Schwinger, and H. A. Nye, Phys. Rev. 
49, 209 (1941). 
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according to the numerical values of the!co- 
efficients appearing in (26) and (29), our method 
of observation is 12 times more sensitive to dis- 
integration-protons than to recoil-deuterons. The 
value we found for the disintegration cross 
section has to be considered an upper limit be- 
cause it would decrease by increasing the as- 
sumed percentage of light-paraffin impurity. For 
this reason we consider it impossible to state 
with certainty, from our experimental results, 
the existence of process (16) for neutrons of 14 
Mev. If it takes place as one would conclude by 
assuming a reasonable value for the light paraf- 
fin impurity, the corresponding cross section can 
be calculated from the total cross section given 
in Table III and the values (31), (31’) of x: 


og= (0.10.86 —0.05 X 0.86) X 10-4 
=(0.09—0.04)10-“ cm?. (32) 


Process (16) was investigated by Bagge* who 
observed in a Wilson-chamber filled with heavy 
hydrogen the spectrum of recoil particles of the 
Rn+Be neutrons. From an irregularity observed 
in the spectrum at about 11.0 mm of hydrogen, 
this author draws the conclusion that process 
(16) takes place for neutrons of energy higher 
than Ey,=3.45 Mev (threshold of the process) 
with a cross section of about 


oa=0.3X10-*% cm’. 


Such a value is between three and seven times 
larger than that we found at 13-15 Mev. 

If we do not assume an irregular dependence 
of the cross section of process (16) from the 
energy of the impinging neutrons, we must con- 
clude that the irregularity observed by Bagge is 
not due to the (m—2n) disintegration of deuteron. 

Barkas and White* have investigated the dis- 


% E. Bagge, Physik. Zeits. 43, 226 (1942). 
% W. H. Barkas and M. G. White, Phys. Rev. 56, 288 
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integration process of deuterons under proton 
bombardment. With protons of 5.1 Mev they 
found a cross section equal to 


0.014 X 10-* cm’. 


The influence of the electrostatic repulsion on 
the value of the cross section has been calculated 
by Hécker who deduced from the Barkas and 
White result the following value of the cross 
section of process (16) for neutrons of 5.1 Mev 


o4=0.02 X10- cm’. 


This value seems to be in reasonable agreement 
with our result as the cross section would in- 
crease by a factor 2-5 on increasing the energy 
of the impinging neutrons by a factor 3. 

The cross section of process (16) has been cal- 
culated by Hicker who found an increase with 
the energy E of the impinging neutron, roughly 
proportional to the square of the difference 
E—Eu(Eu=3.45 Mev threshold of the process) : 
at 14 Mev the absolute value would be 


og=2.5X10- cm’, (33) 


that is, at least 30 times larger than our experi- 
mental value (32). 

Such a high difference between the experi- 
mental and the theoretical value can be recog- 
nized also independently from the experiments 
of section 6, by simple comparison of the total 
cross section at 14 Mev given in Table III 
(0.864+0.028)10-", with the theoretical value 
that can be deduced by adding to (33) the theo- 
retical elastic scattering cross section: this can 
be roughly evaluated to about 0.5-1.010-" 
cm? by extrapolation of the Hécker correspond- 
ing curve. 

We can conclude that, in contradiction with 
theoretical predictions at 14 Mev, and probably 
also at lower energy, the disintegration process 
(16) has a cross section very small with respect 
to the elastic scattering cross section. 
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In this paper there are presented a discussion of the calculation of range-energy values for 
protons in air and aluminum and a tabulation of the results. The calculations have been done 
with considerable accuracy and at sufficiently small energy intervals to allow good graphical 
interpolation. These figures, as well as those for the rate of energy loss in both cases, have 
been compiled over a wide range of energies—up to 10” ev 





N view of the efforts currently being given to 

the production of very high energy protons 
and the general interest in mesons of various 
energies, it will be convenient to have readily 
available a fairly extensive tabulation of range- 
energy relations for such particles. Protons up 
to 10! ev in air and aluminum have been chosen 
in this instance; the range of mesons up to 
10° ev can be deduced easily from Tables I 
and II. Calculations have been made previously 
to this to obtain similar relations"? ; it is intended 
herein to provide data compiled either with 
more careful numerical methods, or with ex- 
tension of certain energy ranges for the sake of 
better accuracy, or both. 

The ranges proposed are those a proton would 
have under the circumstance that it lose energy 
along its path solely through ionization and 
excitation of the atoms of the stopping material. 
These processes certainly constitute the most 
important mechanism for slowing-down over a 
rather wide range of energies; for example, 
meson production will not take place until a 
proton has at least 100 Mev and probably does 
not become important until a much higher 
energy. In addition to total ranges, the rate of 
energy loss by ionization and excitation will be 
of interest; these figures are included in the 
tabulation. 

The figures in this paper represent mean 
ranges since they are calculated from an ex- 
pression which was derived to give the average 
energy loss per unit thickness of stopping ma- 
terial. For a brief review of the theoretical 
derivation, see references 1 and 2. The complete 


1B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 249 


(1941). 
a les and R. Ladenburg, Phys. Rev. 60, 761 
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expression is* 
dE 4a(ez)*e?N 


dx mv? 


2mv? 
x {2 n ’ —In(i—?) -6" — Cet, 





where 


ez =charge of the incident particle, 

v=velocity of the incident particle, 

N=number of atoms per cm? of stopping material, 

Z=atomic number of stopping material, 

I =average ionization potential of stopping material, 
m=electron mass, 

B=v/c, c the velocity of light, 

Cx is a correction term which must be applied i in case v 
is comparable to the velocity of the K-electron of the 
stopping material (but large compared to that of all others). 
Reference 3 gives a discussion of its calculation and figures. 


There are several restrictions placed on the 
validity of this form for dE/dx, the most im- 
portant of which are (1) that the incident particle 
be much more massive than an electron, (2) that 
the incident energy be much less than M*c?/m 
where M is the mass of the incident particle and 


m that of an electron (less than 10° Mev for . 


protons), (3) that the proton have a velocity 
considerably greater than that of L electrons in 
the stopping atoms, and (4) that the proton 
energy be large enough so that electron capture 
and loss are of no consequence. A minimum 
figure suggested for the last of these is 0.1 Mev. 
Thus, all conditions are seen to be well satisfied 
for a proton with an energy between 15 Mev 
and 10‘ Mev in air and Al. 

Let us consider the ranges of protons in air. 
An accurate curve of E vs. R(Z) up to 15 Mev 
for this case is presented in Fig. 28 of the 


*M. S. Livingston and H. Bethe, Rev. Mod. Phys. 9, 
263 (1937). 
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Livingston-Bethe article. To find ranges beyond 
this point, one merely evaluates 

dE’ 
—dE' /dx 
Straightforward numerical integration was found 
to be the only satisfactory procedure. Even for 


E 
R(E) = R(15) +f 
15 


TABLE I. Rate of energy loss and ranges for protons in air. 











E —dE/dz R E —dE/dz R 
(Mev) (Mev/cm) (cm) (Mev) (Mev/cm) (cm) 
15 3.574x<10~? 2.385 X10? 400 0.3308<107? 754.4 X10? 
17 = 3.230 2.975 500 0.2904 1073 
19 2.951 3.623 600 0.2789 1420 
210 «2.721 4.329 700 0.2649 1788 
3 2.529 5.092 800 0.2547 2174 
2 2.364 5.910 850 0.2507 2372 
30 ©: 2.040 8.194 900 0.2472 
35 «1.801 10.81 950 0.2442 2776 
40 = «1.617 13.74 1000 0.2416 2982 
45 1.472 16.99 1250 0.2325 4038 
50 = s:1.352 20.53 1500 0.2277 5126 
60 = «1.173 28.49 1750 0.2252 6231 
70 =: 1.040 37.56 2000 0.2240 7344 
80 0.9393 47.69 2250 0.2236 8462 
90 ©=—-:0.8594 58.83 2500 0.2237 9580 
100 0.7943 70.95 2750 0.2242 10,700 
120 0.6950 97.95 3000 0.2250 11,810 
140 0.6227 128.5 4000 0.2290 16,220 
160 0.5678 162.1 5000 0.2335 20,540 
180 0.5246 198.8 6000 0.2378 24,790 
200 «=—:0.4899 238.3 7000 0.2419 28,950 
250 0.4265 348.2 8000 0.2457 33, 
300 0.3840 472.1 9000 0.2492 37,100 
350 0.3526 608.0 10,000 0.2524 41,080 











TABLE II. Rate of energy loss and ranges for 
protons in aluminum. 











—dE/dz —dE/dz 

E ev R E Mev/ R 
(Mev) (mg-cm~*) (mg-em~*) | (Mev) (mg-cm~*) (mg-cm ~*) 

1 3.45 60 0.8458x10-2 39.83 x 102 

1.5 6.69 70 0.7516 52.40 

2 11.5x107 10.8 80 0.6794 66.42 

2.5 9.85 * 15.6 90 0.6222 81.82 

3 8.62 21.0 100 0.5757 98.54 

3.5 7.69 27.3 120 0.5047 135.8 

4 6.96 34.5 140 0.4530 177.7 

45 6.37 42.1 160 0.4136 224.0 

5 5.88 50.3 180 0 274.3 

5.5 5.47 59.0 200 0.3576 328.4 

6 5.12 69.1 250 0.3120 478.7 

6.5 4.82 79.2 300 0.2813 648.0 

7 4.55 90.0 350 0 833.4 

75 431 101.3 400 0.2428 1033 

8 4.10 113.2 500 2.201107 1.467108 

8.5 3.92 125.6 600 2.054 1.938 

u 3.75 138.8 700 1.952 2.438 

95 3.59 152.4 800 1.879 2.961 
10 3.45 166.7 850 1.851 3.229 
10.5 3.32 181.4 900 1.826 3.501 
ll 3.21 196.6 950 1.802 3.777 
11.5 3.10 212.5 1000 1.785 4.055 
12 2.99 229.0 1250 1.721 5.484 
12.5 2.90 246.1 1500 =: 1.688 6.952 
13 2.8168 263.7 1750—s:11.671 8.441 
13.5 2.734 281.8 2000 1.664 9.941 
14 2.658 .6 2250 1.663 11.44 
15 2.518 3.393 X10? 2500 1.665 12.95 
17 2.281 .228 2750 1.670 14.45 
19 2.089 5.146 3000 1.677 15.94 
21 1.930 6.143 4000 =«11.710 21.85 
23 1.796 7.218 5000 =—s«1.747 27.63 
25 1.682 8.369 6000 1.782 33.30 
30 1.456 11.57 7000 1.815 38.86 
35 1.289 15.23 8000 1.845 44.32 
40 1.160 19.33 9000 1.873 49.70 
45 1.058 23.85 10,000 1.898 55.01 
50 0.9743 28.78 
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the low energies considered it proved to be not 
sufficiently accurate and also too laborious to 
perform the calculation by reducing the integral 
to the difference of standard exponential integrals. 

Values of the constants c, e, m, and N were 
taken from the tables of Birge,* the last being 
made to correspond to conditions of 760 mm Hg 
and 15°C. The proton mass, Mpc*, was computed 
from these tables to be 937.6 Mev. The average 
ionization potential in air, 80.5 ev, and Z=7.22 
are the values as found in reference 3. dE/dx is 
rather insensitive to changes in J; for small 
values of E (~15 Mev), a 3 percent change in J 
is reflected as about 0.5 percent in dE/dx, the 
sensitivity decreasing with increasing energy. 

A similar procedure was followed in calcu- 
lating the range of protons in Al. For energies 
up to 13 Mev, they are obtained from the work 
of Livingston and Bethe; see Figs. 30 and 34 of 
that article and the accompanying text. Beyond 
this point, numerical integration was again em- 
ployed. The same data for Cx as used before 
are here also satisfactory. This correction due to 
K ionization is more important initially for Al 
than for air since the binding energy and velocity 
of the innermost electrons are greater for Al by 
a considerable amount. Experiments by Wilson® 
on the stopping power of Al relative to air 
provide a value for the average ionization poten- 
tial of 150 ev. 

The values for the rate of energy loss and 
range as given in the following tables are believed 
to be correct (consistent with the constants cited 
in this paper) to a few parts in the fourth figure 
where this is given. There is a serious question 
as to the validity of the values of dE/dx in Al for 
very low energies. The condition 


E>(M p/m) X (ionization potential of L electrons) 


is not very well satisfied in the neighborhood of 
E=2 Mev. This indicates that a correction “C,”’ 
should be applied, but it is not available. How- 
ever, in the range above 13 Mev, where ranges are 
obtained using the theoretically computed energy 
loss, Cy will have become completely negligible. 

The author wishes to express acknowledgment 
to Professor H. A. Bethe for various helpful 
suggestions made by him in connection with 
these calculations. 








* For remarks on the validity of dE /dz for E<13 Mev, see discussion at the end 
of the paper. 


*R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
*R. R. Wilson, Phys. Rev. 60, 749 (1941). 
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In a non-stressed cubic crystal the axes of pairs of 
adjacent solute atoms are randomly orientated along the 
several permissible crystallographic directions. It is shown 
that in the presence of a stress the equilibrium orientation 
is no longer random, but that certain directions become 
preferred, and it is further shown that the continual 
striving of the crystal to maintain equilibrium causes it to 
manifest all the characteristic anelastic effects, stress 
relaxation, internal friction, etc. In particular, its internal 
friction will be anomalously high in that temperature 
range where the time of relaxation for the establishment 


of equilibrium is comparable with the period of oscillation, 
The theory is developed for the dependence of anelastic 
effects upon temperature and upon crystallographic orien- 
tation, and it is found that in b.c.c. crystals these effects 
vanish for a tensile stress applied along one of the (100) 
axes. The magnitude of the anelastic effects associated with 
this stress induced preferential orientation may be very 
large. Thus in alpha-brass past experiments indicate that 
under optimum conditions the associated internal friction 
is at least 20 times as large as that due to all other causes, 





O real metal is perfectly elastic in the sense 
that stress and strain are single-valued 
functions of one another. The application of a 
constant stress is always followed by one or more 
types of relaxation, such as the relaxation of 
temperature fluctuations, of shear stress across 
viscous regions, etc., and consequently by a 
gradual creep, which in most cases is recoverable. 
These relaxations are also responsible for stress 
relaxation when a constant strain is applied, for 
a variation of the elastic moduli with frequency, 
and for internal friction. That property of a 
metal in virtue of which stress and strain are not 
unique functions of one another in the pre-plastic 
range has been called anelasticity. The various 
sources of anelasticity have recently been re- 
viewed by the author.' In the present paper an 
analysis is given of one source of anelasticity 
which has heretofore been overlooked, the re- 
laxation of the preferential orientation of pairs 
of solute atoms. 

The solute atoms in a solid solution usually 
have a different size from that of the solvent 
atoms, and therefore give rise to an elastic 
distortion of the lattice. When the solution is of 
the substitutional type, and when the lattice has 
cubic symmetry, the distortion produced by each 
isolated solute atom will also have cubic sym- 
metry. When two solute atoms are adjacent to 


* This research was supported by ORI (Contract No. 
N6ori-20-IV). 
1C. Zener, Metals Tech. (Aug. 1946). 
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one another, the distortion of the lattice by the 
pair will no longer have cubic symmetry. Thus 
if the solute atoms are larger than the solvent 
atoms, the lattice will be extended more along 
the axis of the pair than along any other axis. 
When the average stress throughout the crystal 
is zero, the pair axes will be randomly dis- 
tributed over the several permissible crystal- 
lographic directions. When, however, a stress is 
applied to a crystal, the equilibrium distribution 
of pair axes is no longer random. Thus if the 
solute atoms are larger than the solvent atoms, 
and a tensile stress is applied, the solute pairs 
will tend to be oriented preferentially along 
those permissible axes which make the least 
angle with the tensile axis. This preferential dis- 
tribution of pair axes will, however, require some 
time to be established, and consequently strain 
lags behind stress, i.e., the crystal is anelastic. 
The internal friction in an alpha-brass (solid 
solution of zinc in copper) crystal was observed 
by the author? in 1942 to be anomalously high 
with a narrow temperature band in the vicinity 
of 400°C. The observations are reproduced in 
Fig. 1. At that time no explanation was given 
of this anomaly. Recently an interpretation was 
advanced based on the assumption that the 
equilibrium distribution of zinc atoms was non- 
uniform,! an assumption which has no cor- 
roborative support. In view of the above con- 
siderations, it looks as though the anomaly shown 


2C. Zener, Trans. A.I.M.E. 152, 122 (1943). 
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in Fig. 1 is a manifestation of the stress induced 
preferential orientation of pair axes, and that 
all solid solution crystals will show similar 
anomalies. The temperature variation of the 
internal friction is explained as follows. At low 
temperatures the time of relaxation + for the 
“establishment of equilibrium distribution of pair 
axes is very long compared with the period of 
vibration P. The distribution of pair axes there- 
fore remains essentially unaltered during vibra- 
tion, strain is essentially in phase with stress, 
and hence the internal friction is extremely low. 
At very high temperatures the time of relaxation 
t is very short compared with the period P, and 
consequently the distribution of pair axes is at 
every instant essentially the equilibrium dis- 
tribution. Strain is again essentially in phase with 
stress, and so the internal friction is again very 
small. Only in the intermediate temperature 





range, where 
r=P, (1) 


and hence where the pair axis distribution is 
neither constant nor equivalent to the equi- 
librium distribution, does strain lag behind 
stress appreciably, and hence is the internal 
friction large. 

If the rate at which the distribution of pair 
axes tends to approach the equilibrium distribu- 
tion is proportional to the difference between the 
instantaneous and equilibrium distribution, it 
may be shown! that the internal friction (Q-) 
varies with the angular frequency of measure- 
ment W as 


TW 


— (2) 
1+7°w? 


Q'=4 


Here Ag is defined in terms of the unrelaxed and 
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relaxed Young’s moduli, Ey and Ep, respectively, 
by 
Ag=(Ev—Enz)/Ez, (3) 


and is called the relaxation strength. The ob- 
servations on the brass crystal were found to 
obey this formula provided + was assumed to 
vary with temperature through a heat of activa- 
tion type of relation, i.e., 









(4) 


In Section I an analysis is given of the manner 
in which the relaxation strength varies with 
temperature and with the material constants of 
the metal. In Section II an analysis is given of 
the manner in which the relaxation strength 
varies with the crystallographic orientation, and 
a comparison is given of the relaxation strengths 
in the case of a pure tensile test, Az, and in a 
pure torsion test, Ag. 

The theory of the anelasticity arising from the 
preferential orientation of pair axes is formally 
very similar to that arising from the stress 
induced macroscopic diffusion of solute atoms 
studied by Gorsky,’? and to the stress induced 
preferential distribution of the interstitial carbon 
and nitrogen atoms in alpha-iron studied by 
Snoek‘ and Polder.’ 


r= reeH/RT, 


\ 






















I. TEMPERATURE DEPENDENCE OF 
RELAXATION STRENGTH 





In this section we shall compute the relaxation 
strength for a particular case where the crystal 
axes are so oriented as to simplify the com- 
putations. The influence of crystal orientation 
is reserved for the next section. 

An f.c.c. lattice will be considered in which a 
tensile stress is applied along the [111 ] direction. 
The pair axes of the solute atoms may lie along 
any of the six (110) directions. Of these six 
directions three [011], [101], [110] have a 
direction cosine of (2/3)! with respect to the 
[111] axis, while three, [011], [101], [110], are 
normal to the [111] direction. If the solute 
atoms are larger than the solvent atoms, a tensile 
stress along the axis of the specimen will therefore 
tend to orient the pair axes from the second 

















oe ae, rome ik. Zeits. Sowjetunion 8, 457 (1935). 
~ 6, 591 (1939); 8, 711 (1941); 9, 862 
(19d), my, we Weekblad 39, 454 (1942). 
5D. Polder, Phillips Res. Rep. 1, 1 (1945). 
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group, the group, into the first, the » group, 


of the (110) directions. 

The tensile strain ¢ will be regarded as a 
function of the tensile stress ¢, the temperature 
T, and the number of pair axes in, the p group 
per unit volume m,. The tensile strain could also 
be regarded as a more complex function of 
atomic distribution, such as the relative orien- 
tation of next nearest neighbors. However, it is 
anticipated that only a small error will be intro- 
duced by neglecting the correlation in positions 
of pairs of solute atoms which are not nearest 
neighbors. According to the present assumption, 
a variation in the tensile strain is specified by 


de=Ey'do+adT+ (de/dn,)dny. (5) 


Here Ey is the isothermal unrelaxed tensile 
elastic modulus, i.e., the tensile modulus when 
mM, remains unchanged, and a is the linear 
thermal expansion coefficient when m, remains 
unchanged. We shall now compute the change 
in the elastic modulus occasioned by , changing 
in such a way as to maintain equilibrium. Under 
equilibrium conditions 


y= ne /(1 eV), 


where m is the number of pairs per unit volume, 
and U is the work necessary to transfer a pair 
axis from an m to a p direction, the transfer 
being performed at constant strain and at con- 
stant temperature. If, as is normally the case, 
U/kT is small compared with unity, the above 
equation reduces to 


N,p=n/2—nU/4RT. (6) 


It remains to obtain an estimate of U. It is 
anticipated that U will not depend appreciably 


upon temperature. It will therefore be estimated — 


by considerations at the absolute zero tem- 
perature where thermal effects are absent. The 
change in energy per unit volume is then given by 


= ode+ Udny. 
The differential of E—ge is therefore 
, —eda+ Udny. 
Since this is a perfect differential, 
aU /do = —de/dny. 


We now introduce the assumption that U is a 
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function only of , and is independent of n,. We 
thereby exclude any effects arising from the 
mutual interaction of neighboring pairs of atoms, 
an interaction which may, in certain cases, be of 
importance, perhaps even leading to an ordered 
arrangement of the pairs at zero stress. This 
assumption leads to 


U=—(de/dn,)o. (7) 
When we now substitute Eqs. (6) and (7) into 
Eq. (5), we obtain 
Ep — Ep =n(de/dn,)?/4kT, 
and therefore, from Eq. (3), 


Ag=nEy(de/dn,)*/4kT. (8) 


An estimate of the relaxation strength Ag 
will now be given for the case of alpha-brass with 
30 percent zinc. In this solution® 


de/8C =0.071, 


where C is the zinc concentration. Upon using 
this value, and upon taking m as identical with 
the number of zinc atoms per unit volume, and 
upon replacing (de/dn,) by r(de/dn), we obtain 
for the relaxation strength at 400°C 


Az=0.075r’. 


It is anticipated that the numerical coefficient r 
will be less than, but of the order of magnitude 
of, unity. In the crystal examined Ag was? 0.025, 
a value whose order of magnitude is in agreement 
with the above theoretical computation. 


II. ORIENTATION DEPENDENCE OF 
RELAXATION STRENGTH 


In an experimental study of the anelasticity 
arising from stress induced preferential orien- 
tation of pair axes it is highly desirable that 
single crystals be used in order to avoid the 
anelasticity associated with the viscous’ grain 
boundaries. Since it is difficult to grow a set of 
single crystals all of the same orientation, it is 
desirable that a theoretical formula be at hand 
which gives the dependence of the relaxation 
strength upon crystal orientation. 

In the type of relaxation we are here con- 
sidering no relaxation will be induced in cubic 
_ €E. A. Owen and L. Pickup, Proc. Roy. Soc. 140, 179 


(1933). - 
7 See the coming paper by T. S. Ké. 
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crystals by pure hydrostatic pressures. The 
relaxation under an arbitrary stress system may 
therefore be represented in terms of the relaxa- 
tion in the two shear moduli cq, and (¢1;—¢12) /2. 
The relaxation strengths of these two moduli 
will be denoted by 6 and 4’, respectively. Thus 


5= { (Cas)u — (Cas) e} /(Caa)r- (9) 


In order to compute Ag and Ag in terms of 6 
and 6’, use will be made of the following formulae 
for the tensile and shear moduli: 


E =s,,—[2(s11—S12) — Sas], (10) 
l= $444+2[2(511—Si2) —Saa 1, (11) 

where 
O=7PyvP +22 +781, (12) 


with 71, Y2, Y3 being the direction cosines of the 
specimen axis with respect to the principal axes 
of the specimen. We now subtract Ey from 
Er, and divide by Ey. The _ right-hand 
member will now be transformed from the s's 
to the c’s. Towards this end we introduce the 
notation 

. S=Sq4, C=C, . 
s’=2(sir—Si2),  ¢’ = (C11 — C12) /2, 
and observe that 

s=1f/e; s’=1/e’; sii +2sy2.=1/(C114+2¢12); 

5 = (cu —Ccr)/cr=(Sr—Ssu)/su; 

0 =(c'u—c'r)/c'r=(s'r—S'v)/S'v. 
Upon dividing the right-hand member by ss’, 
and upon using these relations, we obtain 


co’ /3—(cd’—c’d)o 
(13) 





gE = 
Caa(Cir +12) 


———(c—c’) 
2(C11 + 2¢12) 


In a similar manner one obtains 
c'6+2(cd’ —c’d)@ 
Ag= , (14) 
c’+2(c—c’)o 


It remains to determine which, if any, of the 
two unknown quantities 6 and 6’ are zero in 
cubic lattices. In b.c.c. lattices the axes passing 
through nearest neighbors lie along one of the 


8E. Schmid and W.[Boas, Kristallplastizitét (Berlin, 
1935), p. 23. 























four (111) directions. Pairs of solute atoms 
lying along these four directions will be affected 
equally by a shear stress o2:—¢y, across the (110) 
planes, but unequally by a shear stress o,, across 
a (001) plane. The first type of shear stress will 
therefore not cause a preferred orientation of 
pair axes, i.e., 5’ is identically zero, and therefore 
the Ag due to interaction of nearest neighbors is 
zero when a tensile stress is applied along one of 
the principal axes. On the other hand, the shear 
stress oz, will tend to cause a preferred orien- 
tation, and hence 4 is not zero. These conclusions 
may be arrived at in a more elegant manner. 
Let 111, “i11° ++ be the number of pairs of solute 
atoms with axes lying along the directions [111], 
[111], ---. The potential energy of the lattice 
can then contain three interaction terms which 
are linear in the m’s and in the strains. One inter- 
action term will contain the product of the sum 
of the m’s and of the sum ¢€22+€yy+€22. This term 
will cause no relaxation, and hence will not be 
further considered. The remaining two terms 
represent interaction of the orientation variables 
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with shear strains. These terms must have the 
same symmetry as the lattice. One term is 


ox { 10111(€ye t+ €s2 t+ €xy) + i11(€ys — €xz — Exy) 
+n1i1( — €ys t+ €sz— Ezy) +1111( —€Eysz —€z2 +€,y) } f 


No interaction term can be formed of the ¢,,, 
€y, ss Strains, other than that representing a 
uniform dilation, which has the cubic symmetry 
of a cubic lattice. 

In f.c.c. lattices the axes passing through 
nearest neighbors lie along one of the six (110) 
directions. These pairs are affected unequally by 
both types of shearing stress, and hence both 3 
and #6 are different from zero. This conclusion 
is vindicated by the existence of two shear 
interaction terms which satisfy the symmetry 
relations. These are 


B{ (M011 +031) (2€cz — Ey — Exe) 
+ (1101 +01) (Zéyy —€Ess— €zz) 
+ (2110+ i10) (2€s2— Exe — Ey) } 
and 


¥{ (011 — Moi1) €ye t+ (101 — M01) €s2 + (%110 — M10) Exy}. 
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HE problem of the interaction of matter 
and fields has not been satisfactorily solved 
to this date. The root of the trouble in present 
field theories seems to lie in the assumption of 
point interactions between matter and fields. 
On the other hand, no relativistically invariant 
Hamiltonian theory is known for any form of 
interaction other than point interactions. 

Even for the case of point interactions the 
relativistic invariance is of a formal nature only, 
as the equations for quantized interacting fields 
have no solutions. The uses of source functions, 
or of a cut-off in momentum space to replace 
infinity by a finite number are distasteful arbi- 
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It is usually assumed that space-time is a continuum. This assumption is not required by 
Lorentz invariance. In this paper we give an example of a Lorentz invariant discrete space-time. 








trary procedures, and neither process has yet 
been formulated in a relativistically invariant 
manner. It may not be possible to do this. 

It is possible that the usual four-dimensional 
continuous space-time does not provide a suitable 
framework within which interacting matter and 
fields can be described. I have chosen the idea 
that a modification of the ordinary concept of 
space-time may be necessary because the “ele- 


~ mentary” particles have fixed masses and 


associated Compton wave-lengths. 
The special theory of relativity may be based 
on the invariance of the indefinite quadratic form 


S=Cf—x?—y¥—32, (1) 
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for transformations from one inertial frame to 
another. It is usually assumed that the variables 
x, y, 3, and ¢ take on a continuum of values and 
that they may: take on these values simul- 
taneously. This last assumption we change to the 
following: 

x, y, 2, and é are Hermitian operators for the space-time 


coordinates of a particular Lorentz frame; the spectrum 
of each of the operators x, y, 2, and ¢ is composed of the 
possible results of a measurement of the corresponding 
quantity; the operators x, ¥, 2, and ¢ shall be such that the 
spectra of the operators x’, y’, 2’, formed by taking linear 
combinations of x, y, z, and ¢, which leave the quadratic 
form (1) invariant, shall be the same as the spectra of 


x, ¥, , and é. 


In other words, we assume that the spectra of 
the space-time coordinate operators are in- 
variant under Lorentz transformations. It is 
evident that the usual space-time continuum 
satisfies the above definition; it is, however, not 
the only solution. The principal result in this 
paper is that there exists a Lorentz invariant 
space-time in which there is a natural unit of 
length. We hope that the introduction of such a 
unit of length will remove many of the divergence 
troubles of present field theory. 

It can be shown easily that the introduction 
of a smallest unit of length in space-time will 
force one to drop the usual assumptions of com- 
mutativity of x, y, z, and ¢, otherwise the assump- 
tion of Lorentz invariance of the spectra of the 
operators x, y, 2, and ¢, if they commute, implies 
continuous spectra. 

To find operators x, y, z, and ¢ possessing 
Lorentz invariant spectra, we consider the 
homogeneous quadratic form! 


—9? =o? — 9-92? — 93° — 10, (2) 


in which the n’s are assumed to be real variables. 
The variables 49 to », may be regarded as the 
homogeneous (projective) coordinates of a real 
four-dimensional space of constant curvature? 
(a De Sitter space). We now define x, y, 2, and ¢ 
by means of the infinitesimal elements of the 
group under which quadratic form (2) is in- 


1 We could also have taken the quadratic form 
P=ne— n° — 12° — nF +n2. 
This leads to another Lorentz invariant discrete space-time. 
*I am indebted to W. Pauli in a private communication 
for this interpretation of the variables 0, m -- -. 
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variant. We take 
x =1a(40/dn1—110/9n4), 
y =1a(940/dn2—20/An4), 
z= 1a(940/0n3—30/On), 
t = (ia/c)(949/Ono+100/On4), 


in which a is the natural unit of length, c is the 
velocity of light. The operators x, y, z, and ¢ 
are assumed to be Hermitian and may be re- 
garded as operating on single valued functions 
of mo, 1:-:. From the assumption that x, y, and 
z are Hermitian operators of the form (3), one 
can show that each of them has a spectrum con- 
sisting of the characteristic values, ma where m 
is an integer, positive, negative, or zero. The 
operator, ¢, has a continuous spectrum extending 
from minus infinity to plus infinity. The spec- 
trum of each of the operators x, y, z, and ¢ is 
infinitely degenerate. 

Transformations which leave quadratic form 
(2) and 4 invariant are covariant Lorentz trans- 
formations on the variables 7, m2, 93, and mo. 
When the transformed variables 1’, 91’, 72’, 13’, 
n« are substituted into Eqs. (3), it is found that 
x, y, s, and ¢ undergo contravariant Lorentz 
transformation. It is evident that the new 
operators x’, y’, 2’, and ¢#/ which are formed by 
replacing 10, 71, 72, 3, 7 in Eqs. (3) by m0’, m’, 
no’, ns’, m4 are linear expressions with real con- 
stant coefficients in x, y, z, and ¢, and are con- 
sequently Hermitian operators if x, y, z, and ¢ 
are Hermitian. The functions on which these 
operators operate are left invariant except for 
change in argument. Thus, we see that x’, y’, 
z’, and ¢’ which are Hermitian operators of the 
same form in the variables no’, 91’, 2’, 93’, 74 as 
x, y, 2, and-t were in the variables mo, 91, 72, 73, 4 
have the same spectra as do x, y, 2, t. 

We define six additional operators by the 
equations 


(3) 


L,=th(30/dn2—120/9n3), 
Ly =th(10/8ns—139/8m:), 
L,=th(928/8n1—110/8n3), 
Mz =th(08/01+110/8n0), 
M, = th(00/8n2+720/Ano), 
M,=ih(00/dn3+1730/dn»), 


(4) 
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in which h is Plank’s constant divided by 27. 
Now, Lz, Ly, L:, Mz, M,, M, are, with the ex- 
ception of the factor th, the infinitesimal elements 
of the four-dimensional Lorentz group. It may 
be shown that each of the above operators 
commutes with quadratic form (1) when the 
values of x, y, z, ¢ given by (3) are used. This is 
another way of stating the Lorentz invariance of 
quadratic form (1). It will be observed that L., 
Ly, L., Mz, M,, M, do not involve m4, and that 
as a consequence the Lorentz transformation 
which leaves form (2) and 7, invariant leaves 
quadratic form (1) invariant, and as has already 
been shown, leaves the spectra of the operators 
x, y, , and ¢ invariant. We see from these facts 
that the usual assumptions concerning the con- 
tinuous nature of space-time are not necessary 
for Lorentz invariance. This result is the mini- 
mum objective of this work. 

The ten operators defined in (3) and (4) have 
a total of forty-five commutators. Only six of 
these commutators differ from the ordinary ones 
and these six are 


[t, x ]=(ta*/he) Mz, 
[t, y]=(#a*/hc)M,, (5) 
[t, 2] = (ta*/hc) M,. 


[x, y] ag (ia?/h)L., 
Ly, 2] = (ia*/h) Le, 
[z, x ]=(ia*/h)L,, 


We see from these commutators that if we take 
the limit a—0 keeping / and c fixed, our quan- 
tized space-time changes to the ordinary con- 
tinuous space-time. 

The commutators for the quantities L., L,, L., 


M:., M,, M, are those of the infinitesimal ele-_ 


ments of the Lorentz group, and for this reason 
we introduced the factor of <h in their definition. 
Thus, L,, Ly, L, have the usual properties of 
quantum angular momentum. The commutators 
which involve one of the operators x, y, 2, or t, 
with one of the operators L,, Ly, L,, Mz, My, M, 
are independent of a, the natural unit of length. 
Thus all of the commutators pass in the limit 
a—0 to their usual expressions. 

In addition to the ten quantities defined in 
(3) and (4) we define four additional operators, 
having the transformation properties of space 
or time displacement operators, or the equiva- 


lent, of energy and momentum operators.* 
be=(h/a)(m/n), Pe=(h/a) (s/n), 
py=(h/a)(n2/ns), Pe= (he/a)(no/n). 


By means of algebraic manipulations one can 
show that 


1 
L.=yp.s—2py; Mz=-xpitctp.; etc. (7) 
c 


The angular momenta have their usual ex- 
pression in terms of coordinate and momentum. 
One must exercise care as to the order in which 
the factors in (7) are written since they do not 
commute. The four quantities p., py, ps, py; 
commute with one another, and have the same 
commutators with L,, L,, L., Mz, M,, M, as do 
the usual expressions for the space or time dis- 
placement operators. In addition, each has a 
continuous spectrum running from minus infinity 
to plus infinity. Their commutators with the co- 
ordinates and time are not the same as usual and 
are given by 


(x, ps ]=thl1+(a/h)*px*); 

[t, pe] =shl 1 — (a/he)*p,?); 

Lx, by J=Ly, b2]=th(a/h)*pzpy; 

[x, pr]=c*[pz., t]=ih(a/h)*p.p:; etc. 


Here we note that if all the components of the 
momentum are small compared to h/a and the 
energy is small compared to hc/a then these 
commutators approach those which are given in 
orinary quantum mechanics. Further, as we take 
the limit a—0 these commutators change to their 
standard values. 

The fact that these new commutators between 
coordinate arid momentum differ from the old 
ones appreciably only for large values of the 
momentum, and that they differ by large 
amounts when |p|>/a implies that a field 
theory based on quantized space-time will give 
substantially the same results as the usual field 


(8) 


*The most general form for the energy momentum 
operators with the correct transformation properties is 


=2 1 (™) = -"n(*) 
er Td ee ed OF 
in which f(n4/n) is a dimensionless function of its argu- 


ment. A choice of f(n4/n) other than the choice which gives 
Eqs. (6) may be more useful for physical problems. 








6) 


8) 








theory for all processes which do not involve 
large components of the momenta, but will 
produce large effects for processes which do 
involve large components of the momenta. We 
might expect that the usual atomic and molecular 
formulas will suffer no appreciable change, while 
expressions for self-energy, polarization of the 
vacuum, and possibly nuclear forces will be con- 
siderably altered. Alterations in the last men- 
tioned quantities are certainly necessary. 

We note here that the commutator relations 
(7) have a solution for x, y, 2, ¢. 














0 re] a 
| + (0/8)*P6( b+ Py — 
Op: Op: Opy 
0 0 
+); +e—)], 
Op, Op: 
(9) 
0 0 re 
1=in{ —— (a/hc) *po( p+ y— 
Op: Op: Opy 
0 0 
+p5—+)— | etc. 
Op, Op: 


This solution of the commutator relations (8) is 
completely analogous to the solution of the con- 
ventional commutator relations between coor- 
dinate and momentum for the coordinates in 
terms of derivatives with respect to the mo- 
mentum. It would not be surprising to find in 
our case that no corresponding solution for the 
momenta in terms of derivatives with respect to 
the coordinates exists since our coordinates and 
time do not have continuous spectra in general 
and do not commute. 

It is not difficult to verify that if the ex- 
pressions for x, y, z, and ¢ given by (9) are sub- 
stituted in the commutators (5) and the cor- 
responding L,, Ly, L., Mz, My, M, computed, 
then these last mentioned quantities satisfy 
their usual commutation relations. 

It might appear that x, y, z, and ¢ as given by 
Eqs. (9) are not Hermitian operators. It may be 
shown, however, that they are Hermitian opera- 
tors when the correct volume elements of group 
space are used. The group which concerns this is 
the transformation group which leaves quadratic 
form (2) invariant. An infinitesimal volume 
element, dr, of group space is given in terms of 
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Dz, Py, P, and p; by the formula 


pa hdpdpAp Ap: (10) 
ac(p.?+-p,y?+p.?+(h/a)*—(p./c)*) 


To this point the invariance properties of 
quantized space-time have been considered only 
with respect to rotations and Lorentz trans- 
formations which leave the origin fixed. A con- 
tinuum of translations is not admissable in.this 
space, indeed one can prove that if x, y, 2, and / 
are such that quadratic form (1) is invariant 
under infinitesimal displacements, the space x, 
y, z, and ¢ must be a continuum. Translations of 
the origin may be introduced as follows. Let 
S(no, 1, 72, 73, 74) be an arbitrary unimodular 
single-valued homogeneous function of the degree 
zero in the variables, mo, 71, 72, 73, ns. Let A be 
any one of the operators defined by Eqs. (3), 
(4), or (6). Form the new operators A’ =SAS in 
which § means the complex conjugate of S. It 
may then be shown that the primed operators 
satisfy exactly the same commutation relations 
as did the original operators and that they are 
Hermitian. The primed operators thus have 
exactly the same transformation properties with 
respect to Lorentz transformations as did the 
original operators. In particular, the displace- 
ment operators pz, py, Pz, Ps are not changed by 
this process. This is a result which we must 
demand for the space and time displacement 
operators. If we choose S=exp[im arctan(1:/ 14) ], 
we find that this results in a translation of ma 
in the x direction. It should not be expected, nor 
is it possible in this space, to find a function S 
which corresponds to sharply defined transla- 
tional values of x, y, z, and ¢ simultaneously. The 
relation between two different frames of reference 
cannot be set up more precisely than the com- 
mutation relations between the coordinates 
permit. 

We will not discuss in any detail in this paper 
the limitations placed upon the simultaneous 
measurability of x, y, z, and ¢ due to the non- 
commutativity of these quantities. Some pre- 
liminary calculations which I have made indicate 
that these limitations are not serious enough to 
interfere with the ordinary description of atomic 
phenomena in terms of a continuous space-time 
nor with our usual macroscopic theory. 
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The problem of the normal reflection of a shock wave is solved by an approximate analytical 
integration of the hydrodynamical equations. The solution given here leads to nearly the same 
numerical results as an exact method, based on a graphical integration of the hydrodynamical 
equations, which has been described by Chandrasekhar. A method of computing the complete 
pressure-time curve at the reflector is given and applied to reflection in a class of fluids obeying 
the Tait adiabatic equation of state. It is found that in compressible fluids (gases) the pressure 
on the reflector is prolonged and that the impulse delivered to it exceeds the value predicted 
by the acoustic theory. In slightly compressible media (liquids and solids), on the other hand, 
the blow is shorter and the impulse delivered to the reflector is Jess than one would expect from 
the acoustic approximation. The method given here is also applicable to the reflection of 























gravity waves on the surface of a liquid. 














INTRODUCTION 





HEN a shock wave is reflected from a solid, 

the pressure on the reflector at first rises 

to a high peak and then decays to the value 

previously existing. For weak shocks, or sound 

waves, which strike a rigid reflector, this initial 

peak is of such a height that the pressure ratio 

({-) across the reflected shock front and the 

pressure ratio ({;) across the incident shock front 
are equal, 
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so that the overpressure at the reflector is twice 
the overpressure in the incident shock. In the 

- limiting case of very strong and normally incident 
shocks, on the other hand, 


§=[(3y-1)/(y—-1) Ik. 


For a fluid whose adiabatic constant, y, is 1.4, 
,/{; is 8. The rise in pressure at the first moment 
of reflection and the flow patterns simultaneously 
present have been studied at all angles of inci- 
dence and the theory is able to describe in an 
approximate way both the regular and the 
irregular, or Mach, reflection. One is, however, 
rarely interested in the initial, or peak, pressure 
alone unless the period of the reflector is very 
short compared to the duration of the incident 
wave. Similarly in the collision of equal and 
opposite shocks, which is hydrodynamically 
equivalent to reflection from a rigid wall, it is not 





















*Now at Argonne National Laboratory, Chicago, 
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just the first instant of the collision but rather the 
subsequent history of the interaction which is 
interesting. It is the complete course of the 
reflection of one shock, or of the interaction of 
two shocks, which we wish to discuss here. 


QUALITATIVE DESCRIPTION OF RESULTS 


For the case of normal incidence alone has the 
problem been solved. The possibility of using a 
very simple numerical method based on the 
Riemann form of the equations of motion was 
indicated by Chandrasekhar! who considered re- 
flection in air. Here we replace the numerical 
method by an analytical one and extend the work 
to other media besides air. All media covered by 
our calculations are subject to the Tait equation 
of state for isentropic changes. 


(p+m)v7 =k. (1) 


Here is overpressure ; v is specific volume; 7, 7, 
and k, which are characteristic of the medium, 
are constants for isentropic processes but are 
functions of the entropy. When Eq. (1) is applied 
to ideal gases, y is the adiabatic constant lying 
between 1 and 1.67. The same equation is also 
applicable to matter in the condensed state, if + 
is now interpreted as the internal pressure and y 
is an empirical constant, not simply related to the 
specific heats: e.g., for water, r=3000 atmos., 
y=7.15; for freshly detonated explosives which 


1S. Chandrasekhar, a report of limited circulation issued 
by Ballistics Research Laboratory, Aberdeen, Mary- 
nd. 
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REFLECTION OF SHOCK WAVES 


have not yet expanded, y=3. It has also been 
pointed out that an analogy may be set up be- 
tween the equations of gas dynamics and of 
gravity waves on the surface of a fluid and this 
analogy corresponds to the case of y=2; our 
results are in fact applicable to gravity waves if 
pressure is interpreted as h?, where h is the height 
of the gravity wave.” 

Although the numerical method used by 
Chandrasekhar is rigorous, the analytical method 
employed here contains certain approximations 
and the two methods were, therefore, compared 
for an incident linear shock (i.e., the pressure falls 
linearly with time behind the peak) of 1.5 atmos. 
in air (the particular example computed by him) 
and an incident linear shock of 1800 atmos. in 
water. The pressure time curves at the wall are 
shown in Fig. 1. The two methods agree well and 
differ considerably from the acoustic approxi- 
mation. 

The results of our calculations can be most 
conveniently expressed in terms of J/g pdt, the 
time integral of the pressure at the wall, where T 
is the total time required for reflection. This 
integral is the momentum delivered to the wall or 
the difference in momenta between the incident 
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Fic. 1. Pressure-time curves at reflector. The curve at 
the left corresponds to reflection in air of an incident pulse 
whose peak overpressure is 0.5 atmos. The other curve 
describes reflection in water of an incident pulse whose 
peak overpressure is 1800 atmos. The ordinate is over- 
pressure in atmos. and the unit of time is the duration of 
the incident pulse. © Computed by Chandrasekhar’s 
method ; — analytic approximation; —— acoustic approxi- 
mation. 


?T. von Karman, “Flow in Compressible Fluids,” to be 
found in a collection of papers entitled Fluid Mechanics 
and Statistical Methods in Engineering (University of 
Pennsylvania Press, Philadelphia, 1941), page 25. 


and reflected pulses, i.e., 


f pr,dx = f paudx+ f ‘pat (2) 


where p and u are density and material velocity, 
respectively, and where the subscripts 7 and r 
refer to incident and reflected pulses. The space 
integration is extended over the whole lengths of 
the incident and reflected shocks. In Fig. 2 the 
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Fic. 2. Impulse given to reflector as function of y. 
Ordinate is actual impulse given to rigid wall divided by 
approximate acoustic impulse. Abscissa is constant (y) 
appearing in the equation of state. ¢ is shock intensity. 
¢=1.0 corresponds to limiting case of sound wave. 


impulse /o" pdt is normalized by dividing by the 
acoustic impulse, i.e., 2/p.dt, where p; is the 
value the over-pressure would have at the posi- 
tion of the wall if the incident wave moved 
forward without change of shape and with the 
velocity of sound corresponding to zero over- 
pressure. The intensity of the incident shock is 
also normalized in Fig. 2: it is specified by 
$=(p+7)/z. In Fig. 2 the normalized impulse is 
plotted as a function of y; the most interesting 
physical result of our calculation—the difference 
of behavior shown by compressible and incom- 
pressible media—is then brought out: in the case 
of gases (1 < y £1.67) the impulse is greater than 
the prediction of acoustic theory ; but in the case 
of liquids and solids the impulse is Jess than the 
acoustic approximation. Non-linear theory hap- 
pens to agree with the acoustic estimate in the 
neighborhood of y=2; y=2 also corresponds to 
the analogy between gravity waves and gas flow. 
The reason for this difference between gaseous 
and condensed matter is that the duration of the 
blow exceeds the acoustic approximation in the 
gaseous case, whereas the situation is reversed for 
solids and liquids, as one can see by the following 
qualitative consideration of the reflection process. 
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Fic. 3. Reflection in x space. 


After the instant of impact the head, H, of the 
shock wave reverses direction and begins to 
travel toward its tail, T (Fig. 3). As soon as the 
direction of H has been reversed by the wall, the 
total disturbance may be regarded approximately 
as a (non-linear) composition of two waves of 
finite amplitude. One of these moves toward the 
wall and is bounded by T and W. The other 
moves away from the wall and is bounded by W 
and H. These two waves penetrate each other 
until TZ reaches the wall and thereby marks the 
end of the reflection. It is clear that the tail T 
moves toward the wall at speed, a, the velocity of 
sound corresponding to p=0, until it meets H. It 
is not so easy to guess the speed of H with respect 
to the wall, but it also turns out to be nearly a for 
the shock strengths of interest here. The reason 
is, roughly speaking, that the high pressure 
behind Hf produces a supersonic velocity not with 
respect to the wall but with respect to fluid 
rapidly moving toward the wall. The net effect is 
that H and T approach each other approximately 
with velocity a. Hence if the initial thickness of 
the pulse is Z, the time until H and T meet is 
L/2a. After passing through the head the tail 
then moves with the velocity c—u, where c is the 
local velocity of sound and u is the material 
velocity of the fluid in the region M. Hence the 
total time of reflection is L/2a+L/2(c—u)» where 
(c—u)w is the average velocity of the tail in M. 

_On the other hand the duration according to the 
acoustic approximation is L/a. Hence if u/c is so 
large that (c—u)y <a, as in air, the pressure on 
the wall lasts for a longer time than acoustic 
theory indicates; but if u/c is so small that 
(c—1u)w ><a, as in water, then we have the oppo- 
site, situation. These qualitative considerations 
are supported by the detailed calculations to 
which we now turn. 


ANALYTICAL METHOD 


Consider a plane shock, S, impinging on a rigid 
wall, W, at normal incidence. Let the plane of the 
wall be x=0 and let the x axis be directed toward 
the oncoming shock. The time, ¢, is measured 
from the moment at which S strikes W. Since the 
problem is one-dimensional, there are only two 
independent variables, x and ¢, and it is con- 
venient to speak in terms of the x, ¢ plane shown 
in Fig. 4. Here D is the world line of the reflected 
shock front, and the world line of a typical 
element, initially at xo, is shown with a dis- 
continuity where it crosses D. My and M are the 
regions between D and the x axis and ¢ axis, re- 
spectively. The following boundary conditions 
are given on the axes and on D: 


1, Distribution of pressure, ~;, and velocity, u;, on the 
x axis. This is the distribution of pressure and velocity in 
the incident wave. 

2. u=0 on the ¢ axis since the wall is rigid. 

3. Shock equations across the discontinuity, D, ex- 
pressing the conservation of (a) mass, (b) momentum, 


(c) ‘energy. 


Subject to these boundary conditions three differ- 
ential equations which again express the con- 
servation of mass, momentum, and energy must 
be satisfied in Mo and M. In these differential 
equations appear three unknown functions, e.g., 
pressure, velocity, and-entropy, and the problem 
is to find these functions in the regions My and M 
and in addition to find the equation of the curve 
D. Employing the usual approximation, however, 
we assume that the flow in My and M is isentropic. 
In addition the stronger approximation is made 
that the entropy change across D can be neg- 
lected. The entropy change across D is actually 
only third order in the volume change; but a 
stronger argument is that results of the calcula- 
tion made under this assumption agree very well 
with the results of the rigorous numerical method 
in which the change of entropy is properly taken 
into account. Under these assumptions the 
entropy becomes a constant of the motion and 
may be ignored. There are then only two unknown 
functions, say, p and u, and the exact conserva- 
tion of energy condition is replaced by the 
approximate conservation of entropy condition. 

The three conservation conditions in M and Mo 











may then be written 


a dp 
—(pu)+—=0, (3) 
ax ot 
op ou = au 
Ox ot Ox 
p+x=kp’ (an approximation), (S) 


where p is density. There are two unknown func- 
tions and two differential equations. On the other 
hand the three conservation conditions across the 
shock front are 





ee 
somstule =|, 6) 
Vo—V 
u=uot[(p—po)(vo—v) }*, (7) 


(potm)vo’ =(p+7)v” (an approximation), (8) 


where s is the shock velocity (or slope of D) and 
v=1/p is the specific volume. The variables in 
these equations are of course to be evaluated on 
D; subscripts refer to the two sides (M and Mo) 
of D. 

By Riemann’s method’ Eqs. (3)-(5) are re- 
written 


oP oP 
—=—(c+a)—, (9) 
, Ox 
Spun, (10) 
ot Ox 


where P and Q are the Riemann functions which 
are defined in this paper as 


P=2c/(y—1)+u, (11) 


Q=2c/(y—1)—4, (12) 
where 
c=(dp/dp)*= velocity of sound. (13) 


It is convenient to rewrite the shock equations (6) 
and (7) with the aid of (8) and (13) in terms of 
these five variables: c, u, Co, %o, and s. They 
become to the first order in c/co—1 


S=uUpt+ect+eoco, (14) 


Q=Qo, (15) 


3See, for example, H. Lamb, Hydrodynamics (Cam- 
bridge University Press, 1932), sixth edition, page 481. 
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Fic. 4. Reflection in x, ¢ space. OD and EF are world 
lines of reflected shock and tail of incident shock, re- 
spectively, and ab is a characteristic on which Q is constant. 


where the following abbreviations have been 
introduced 
e=(y+1)/2(y—1), 


éo=1—e. 


It has been pointed out that Q is nearly constant! 
across a rather strong shock front in air, even 
without neglect of the entropy change. 


INITIAL CONDITIONS 


We come now to the formulation of the bound- 
ary conditions on the x axis or, in other words, to 
a description of the incident shock wave. For , 
simplicity it is assumed that at ¢=0 the incident 
wave is progressive and that the velocity distri- 
bution in it is linear; the pressure distribution is 
then fixed. Equations (9) and (10) show how the 
progressive condition may be expressed. Ac- 
cording to these equations a general one-dimen- 
sional hydrodynamical disturbance may be 
decomposed (non-linearly) into two waves, de- 
scribed by P and Q which travel in opposite 
directions with the velocities c+u and c—u, re- 
spectively. The condition for a progressive wave 
is that either P or Q be constant; and in this case 
it must be P since the pulse is traveling toward 
—«, The required boundary conditions on the 
x axis are then 


uo(x, 0) = —B+Bx 0<xM1, (16) 
uo(x, 0) =0 x>1, B>0O, 
Po(x, 0) =2/(y—1). (17) 


Here B is the peak velocity. The units of length 
and time are so chosen that the initial length of 
the incident pulse is unity and so that c=1 when 
p=0. Equation (17) together with (5) and (13) 
determine the pressure as the following function 





















46 





of the velocity 
ptr=a[1—(y—1)u/2 Vo, (18) 


The pressure computed from (18) and (16) is very 
nearly linear for shocks in which (p+7)/x<1.6. 
In other words a pulse of the type considered here 
is essentially linear in both pressure and velocity. 

A remark may be made about the progressive 


assumption at this point. Although a continuous: 


wave of finite amplitude can propagate without 


increase of length, the same is no longer true as _. 


soon as it has evolved into a shock. After it has 
become a shock, reflection of the continuous por- 
tion of the pulse from the discontinuity begins; as 
a result the disturbance propagates backwards as 
well as forwards, and hence the progressive as- 
sumption cannot be rigorously satisfied by a shock 
any time after it is formed. Nevertheless it is 
often a good approximation, and as we have just 
seen, the pulse considered here is essentially linear 
in both pressure and velocity. 

Although one could specify the pressure dis- 
tribution in the incident pulse instead of the 
velocity distribution in it, we have not done so 
because it is analytically more convenient to 
regard the velocity as given. Since, however, ex- 
perimental data generally concern pressure, it is 
necessary to have a way of determining the maxi- 
mum velocity B from the measured peak pres- 
sure. B can be calculated from the following 
formula which follows from (7) 


B=((¢—-1)(1-5-")/y}', (19) 


where {=(pn+7)/x and pm is measured peak 
pressure in the incident wave. This Eq. (19) 
is not consistent with (18) since (19) is a shock 
equation while (18) is based on the progressive 
assumption; and as we have just seen, they are 
not compatible. The procedure followed here is 
this: ¢ is regarded as computed from the meas- 
ured peak pressure. From (19) B is found. From 
(16) u is found. Equations (16) and (17) then 
completely determine the wave. The numerical 
discrepancy between (18) and (19) is actually 
insignificant for shocks of the intensity considered 
here ({ < 1.6). 


SOLUTION IN M, 


The differential equations can be solved readily 
in Mo if the disturbance in My remains pro- 
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gressive. That it does remain nearly progressive 


may be verified by comparing the result of the 
rigorous Riemann numerical integration with the 
analytic progressive solution. Therefore we 
assume 


P(x, t) =2/(y—1). (20) 


Let 
A=Co— Uo. | (21) 
Then by (20), (21), (10)—(12) it follows that 
0A dA 
—=A—. (22) 
ot ox 


The solution of (22) is 
t4+x=w(A), (23) 


where w is a function to be determined by the 
boundary conditions. When these conditions are 
(16) and (17) then 


A=[((y+1)(—Bx+B)+2)/[B(y+1i+2] 
x<1—t, (24) 
=1 x>1-2. 


The discontinuity in the derivatives of A along 
the line x+¢=1 (Fig. 4) has its origin in the 
discontinuity at the tail of the incident pulse 
(Eq. (16)). The straight line EF is the world line 
of the tail of the incident shock and F represents 
its intersection with the head of the reflected 
shock. Here we are interested only in the part 
OEF of the region M: and so only the first 
representation of A in (24) will be carried. The 
quantities %o, Co, and Q» are now readily found. 


uo=2B(x+t—1)/(2+B(y+1)t) x<1-t, (25) 
Co=1—(y—1)uo/2, (26) 
Qo = 2¢0/(y—1) — uo. (27) 


With these expressions the problem is solved 
in M 0- 


SOLUTION IN M 


To determine the solution in M one has to 
overcome a difficulty characteristic of shock 
problems; namely, although the boundary con- 
ditions on D are known, the equation of the line 
D is unknown and has to be determined simul- 
taneously with the solution of the differential 
equations in M. The shock equations (14) and 
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(15) contain five variables: uo, Co, u, c, and s, of 
which uo and co are known from (25) and (26). 
Still one more condition is necessary to fix u, c, 
and s, and this additional information is bound 
up in the solution of the differential equations. In 
order to separate the determination of the shock 
line from the solution of the differential equations 
we make the following approximation 


P(g(t), t) =P(0, 0), (28) 


where 


x = g(t) (29) 


is the equation of D. The best arguments for the 
condition (28) are that the decay of P along D, 
computed by the rigorous numerical method, is 
small and that the pressure-time curve, obtained 
on the basis of (28), agrees very closely with the 
corresponding curve computed according to the 
numerical method (as shown in Fig. 1). The 
function P actually decreases along D as ¢ in- 
creases and this decrease causes the shock to 
decay. Therefore the physical meaning of the 
approximation (28) is that decay of the reflected 
shock is negligible during the period of reflection. 
Aside from its numerical success this assumption 
has in its favor two facts which make it plausible: 
first, the time of reflection is short; and second, 
since the shock travels faster with respect to Mo 
than with respect to-the wall, the wall behaves as 
a sustaining piston behind the shock. 

The equations (14), (15), (25), (26), and (28) 
lead to the following differential equation for the 
shock line D 


dx f+hx+jt 

dt 1/B+(y+1)t/2’ 
f=1/B+(3y—5)/4, 
h=(3—y)/2, 
J=2+B(y+1)?/8. 


(30) 





The solution of this equation is 
%=[2/(y—1)B+(Sy—3)/4(y—-1)] 
X[1— (1+ (y+1)Bt/2)@-M/04 7] 
+(2/(y—1))[1+B(y+1)?/16]¢. (31) 


It is more convenient and numerical comparison 
with the complete equation again shows that it 


is little less accurate to use only the first two 
terms of the expansion of (31) in powers of ?, i.e., 


x=at+ fi’, (32) 
where 


a=(37—5)B/4+1, | (33) 
B=((Sy—3)B*/4+2B](3—y)/4. (34) 


Equation (32) fixes the boundary line D. The 
boundary conditions on it are 


P(x, t)=P(0, 0), (35) 
Q(x, t) =Qo(x, t). (36) 


Qo is given by (27) and (32). If the value of Qo so 
obtained is put in (36) one gets 


Q(x, t) =got+gqitt+qol*, (37) 


where 


qo=2(1/y—1+8B), (38) 


qi = —4B[1+(Sy—3)B/8], (39) 
g2= —B(3y—1)qi/4. (40) 


In (37) higher order terms were again found small 
enough to neglect. Finally 


P(x, t) =P(0, 0) =Q(0, 0) =go. (41) 


The boundary conditions are most convenient in 
the forms (41) and (37). 

The boundary D and the boundary conditions 
on it have now been found so that the differential 
equations in M can be solved. Numerical inte- 
gration shows that in M the characteristics of Q 
are almost straight. Now the actual slope of a Q 
characteristic is c—u. Let the slope be assumed 
constant and equal to —A which is defined as 
follows 


2A = (c—u)at(c—u)s, (42) 


where the points a and 6 are shown in Fig. 4. One 
then finds 


(c—u),=co=(y—1)Q./2 =(y—1)Q./2, 
and by (41) 
(c—u)e=(y—3)P./4+(¥+1)Q./4 
= (y—3)g0/4+(¥+1)Q./4. 


Hence 


A= (y¥—3)go/8+ (37 —1)Q./8. (43) 


Then the equation of the characteristic passing 
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through (0, 4) is 
x=X(h—2), (44) 


where \ is given by the preceding equation. 

To find Q, as a function of # it is convenient to 
regard ¢, as the independent variable. The rela- 
tion between Q and # is then given by the 
parametric equations: 


Qs= got gita+ gata’, (45) 
ty =tat+x(te)/A(ta), (46) 
according to (37) and (44), where 

x(ta) = ate+ ft’, 

Ata) = (y —3)go/8+ (37 —1)Qa(ta)/8, 
according to (32) and (43). The time ¢, is the 
retarded value of the time & in the sense that 
t—t, is the time required for a given value of Q 
to propagate itself from the shock front to the 


wall. If the shock and Q both traveled with the 
velocity of sound (unity), then one would have 


exactly #,=0.5t. This relation is nearly satisfied 


‘in any case. 
CALCULATION OF PRESSURE ON REFLECTOR 


In order to calculate (0, &) from # it is suffi- 
cient to find Q(0, &); because at the wall, where 
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u=0, 


c= (y—1)0/2, (47) 


and from c the pressure follows according to the 
equation 


pPtrareVo-d, (48) 
It is convenient to combine (47) and (48) 
pb+mr=a[(y—1)Q/2]}?/o-». (49) 


Since Q(0, t) is already known by (45) and (46), 
this equation completes the determination of 
p(0, &) as a function of 4. 

The complete method for computing the pres- 
sure-time curve at the wall then runs as follows, 
The incident shock is specified by its pressure- 
ratio, ¢. The maximum material velocity, B, is 
first found from Eq. (19). The following con- 
stants, go, 91, g2, a, and 8, which are functions of 
B and y only, are then found from (38)—(40), 
(33), and (34). These constants then determine 
the auxiliary functions x(t.), Q(ts), and A(t.) 
which are defined in Eqs. (32), (37), and (43). In 
terms of these functions the parametric equations 
of the pressure-time curve are (46) and (49). The 
duration of the pressure on the wall is determined 
by putting Q=2/(y—1) in (37) and then 
using (46). 
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The author recently proposed a generalized classical electrodynamics in which the usual 
assumption that the field of a moving electron is given by its retarded field is abandoned, and 
a combination of retarded and advanced fields is considered. This paper considers the applica- 
tion of the generalized equations of motion to investigate the motion of the electron in the 
hydrogen atom. It is shown that, unlike the Lorentz-Dirac equations, these generalized equa- 
tions permit the electron to spiral round and fall into the nucleus. 





I. INTRODUCTION 


N some recent papers I have considered the 

application of Dirac’s classical theory of 
radiating electrons' to various familiar problems; 
and the evidence gathered from these applica- 
tions is that the equations of motion do not 
always have physically understandable solu- 
tions. There are, for instance, the self-accel- 
erating motions that were considered by Dirac. 
Again, if an electron is moving along the x axis 
in an electric field of potential u|x|, that is, in 
the field due to a thin infinite charged plate in 
the yz plane, then the only motion allowed by 
the equations of motion is one in which the 
electron ultimately accelerates away from the 
plate with rapidly increasing speed, despite the 
attraction of the plate.? Further, when one 
investigates the motion of the electron in the 
hydrogen atom, it is seen that the electron does 
not spiral round and fall into the nucleus owing 
to loss of energy by radiation, as one would 
expect, but spirals outwards and escapes to 
infinity with increasing speed. The attraction 
of the nucleus has the effect of accelerating the 
electron away from the nucleus.* All these 
results point to the suggestion that the signs of 
certain terms in the Lorentz-Dirac equations 
do not fit in appropriately. 

With a view to bring about the necessary 
alterations in sign, I have developed a generalized 
form of classical electrodynamics,‘ by abandoning 
the usual assumption that the field of a moving 
electron is given by its retarded field, and con- 


. A. M. Dirac, Proc. Roy. Soc. (A) 167, 148 (1938). 

. J. Eliezer, Bull. Calcutta Math. Soc. 37, 125 (1945). 
. J. Eliezer, Proc. Camb. Phil. Soc. 39, 173 (1943). 

J. Eliezer, Nat. Acad. Sci., India (1946). 


sidering, instead, a combination of retarded and 
advanced fields. If F., denotes the field due to 
an electron, and Fy, and Fyay the retarded and 
advanced fields, respectively, then we take 


Foi = Free +h( Fret — Faav)s (1) 


where k is a constant, the value of which is left 
arbitrary for the present. The equations of 
motion corresponding to the assumption (1) are 
then obtained in the usual way from the prin- 
ciples of conservation of energy, momentum, and 
angular momentum. 

Suppose there are a number of electrons inter- 
acting with each other and with a given external 
field that is described by the field quantities 
Fext*”. If Zun(Sa), (u=0, 1, 2, 3), denote the time 
and space coordinates of the mth electron of 
charge e, and rest mass m,, and dots denote 
differentiation with respect to the proper time 
S,, and the metric tensor g*’ is such that g®=1, 
gi=g*%=g%*=—1 while all the other com- 
ponents vanish, then the equations of motion 
of this electron are 


qe women eH(Z) | 
n — Fen igen n 
" ds ds* ds ™ 


= €,Ven|_F*, ext ‘ (2) 
+ ) { Fam, rot + (Fm, rot — Fm, adv) } J; 


min 





Usa? = 1. J 


As mentioned earlier the constant & is arbi- 
trary. When k=0, Eqs. (2) reduce to the 
Lorentz-Dirac equations. When k= —}3, then 


Fei=3(Free+ Faav), (3) 
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and the equations of motion (2) become 
MV yn/ds = Cen F%y, ext 

$4 L (Fem rot Fym,adv)], ¢ (4) 
Une? = 1. 


Equations (4) do not contain the usual damping 
terms. If k< —4, then the sign of the radiation 
damping term is opposite to that of the cor- 
responding terms in the Lorentz-Dirac equations, 
and in the particular case when k= —1, the 


equations are the same, except for this dif- - 


ference of sign. 


Il. APPLICATION TO THE HYDROGEN ATOM 


If we take the nucleus of the atom to be fixed 
(by considering its mass to be infinite), then with 
reference to axes whose origin is at the position 
of the nucleus, the equations of motion of the 
electron are 


mo -4eQk+1) —+|(<) - ) E 


n dsr 


oa 


where the (?, x, y, z) notation is now used. When 
k>-—4, the behavior of the solution of Eqs. 
(5) will be similar to that of the solution of 
the Lorentz-Dirac equations, which have al- 
ready been considered and which do not allow 
the electron to fall into the nucleus. When k = —}, 
Egs. (5) reduce to the usual equations of motion 
without the radiation damping terms, and the 
solutions of these equations are well known. 

In this paper we shall consider the behavior of 
the solution of Eqs. (5) when k< —}. It will be 
seen that the equations then allow the electron 
to spiral round and approach the nucleus closer 
and closer and ultimately fall into it. 


Ill. RECTILINEAR MOTION 


We shall first consider the simpler case of recti- 
linear motion of an electron moving along the 
x axis, the proton being situated at the origin of 
coordinates. The equations of motion (5)gthen 
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give 


d’x 
a—— (2k+1) 
ds? 


«(i aaa) E+(Z) J] 


fe(Q)T oe 


where a=3m/2e?. It is unfortunately not pos. 
sible to obtain an explicit solution of this 
equation. But our purpose here is to know 
roughly the nature of the motion, and this can 
be done without solving the equation. We write 
dx/ds=u and transform the independent vari- 
able from s to x. Then d’x/ds?=uu', d*x/ds* 
=u(uu'’+u'*), where dashes denote differen- 
tiation with respect to x, and Eq. (6) becomes 
u’? ! all 3 G+" 





(ak-+1){w"+ 


To discuss the behavior of the solution of Eq. 
(7) it is convenient to transform from u to y, 
where 


=(1+4?)}, (8) 


and y is taken to have positive values only. The 
cases u positive and u negative lead to different 
differential equations. For u>0, we have 


(y?—1)4y””+Ay’+ux-*=0, (9) 


where A= —a/(2k+1)>0, p= —3/2(2k+1)>0; 
and for u<0, we have 


(y?—1)4y” —Ay’ — wx? =0, (10) 


where \ and yz are as before. 

Suppose initially the electron is projected 
towards the proton from the point x=x». Its 
motion immediately afterwards is given by Eq. 
(10). The electron may be brought to rest at 
some point x =8 before it reaches the proton, or 
it may collide with the proton. In the former 
case the solution will be defined only in 
B<&x<€xo, and in the latter case in OS x< x. 
Equation (10) is such that if y’=0 at any point 
where the solution is defined, then y’’>0; and 
therefore y has no maxima and therefore at most 
one minimum. Also, if y’’=0, then (y?—1)y’” 
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+2ux-*=0, and hence y’’’<0. Therefore y’’ =0 
at most once in 0<x<xo. We may therefore 
choose an interval O<x<y so that y’ and y” 
are of fixed signs in the portion of the interval 
in which the solution is defined. 

Suppose y’ <0. Then y increases as x decreases. 
Hence, as x0, y tends to a finite limit 1(>1) 
or to infinity. It is easy to see that y cannot tend 
to a finite value /, because if it does, the equation 
obtained by integrating Eq. (10) with respect to 
x gives 


dx 
‘—\ cosh y=p f ———+ constant 
. x2(y*—1)! 


(11) 
= a +&)(?—1)—'+constant, 
x 


where 0 as x—0. Integrating again with 
respect to x we see that the left-hand side 
remains finite as x0, but the right-hand side 
becomes infinite, being approximately of the 
order u(l?—1)—} log (1/x), for small x. We thus 
obtain a contradiction. Hence, if y’<0, the 
solution is such that y © as x—0. 

We have already seen that if initially y’ is 
negative, then as x decreases, y’ will continue to 
be negative, since y has no maximum. The con- 
dition that y’ is negative is equivalent to the 
condition that the acceleration is towards the 
proton, since d*x/ds*=uu'=yy’. We thus have 
the result that if the electron is projected with 
initial velocity and acceleration which are both 
towards the proton, then the electron continues 
to approach the proton with steadily increasing 
velocity and ultimately collides with the proton, 
with velocity tending to that of light. 

It is of interest to obtain an approximate 
solution valid for small values of x. To do this, 
we transform from x to r where x =e~". We then 
have 


ot 12 4 ety pret] =0. (12) 
dr dr 


To obtain a solution which is valid for large 
values of 7, we note that y+ and dy/dr>0 as 
t+, It is seen that dy/dr cannot tend to 
infinity, because if it did, then from Eq. (12), 


d*y/dr* must tend to minus infinity; hence dy/dr . 


must be decreasing with increasing +, which 


gives a contradiction. Hence dy/dr—k20 as 
7—o, and therefore d*y/dr?—— k. 

But according to a theorem on limits, if a 
function f(x) tends to a finite limit as x tends to 
infinity, and f’(x) is known to tend to a limit, 
then this limit must be zero. Hence k=0. Thus 
dy/dr—0 and d*y/dr*-0 as r> &. 

For large values of r, Eq. (12) is approximately 


(13) 


Writing dy/dr =p and changing the independent 
variable to » and the dependent variable to y 
we have 


iM 
P—+ p _ar 0, 
dy y 
which is easily seen to have approximate solution 
pb=n/y, and hence y?=2u7r-+ constant. Hence 


dx 1\3 
—a = (2 log -) , for small x, (14) 
ds x 


where p= —3/2(2k+1). 

If initially y’>0, that is, if the initial accel- 
eration is directed away from the proton, then 
in the subsequent motion, y’ may change sign 
soon afterwards or may continue to remain 
positive until the electron is brought to rest at 
x=B. If y’ changes sign and thus becomes 
negative, the subsequent motion will then be as 
discussed above, namely a motion in which the 
electron continues to approach the proton with 
steadily increasing velocity and ultimately col- 
lides with it. If, however, y’ remains positive 
until y reaches the value 1, then the electron is 
brought to rest at x=8 before it reaches the 
proton. The approximate solution near x=8 is 
of the form 


1+1( + (n+) (15 
j= x -B 3/3 B (x B ’ ( ) 


where y—>1 and y’—/(>0) as x8. If, however, 
y’—0 as x—>8 then the approximate solution near 
x= 6 is given by 


~1)ta- “wep. 16 
(y—1) rr) B) (16) 

















The electron being thus brought to rest before 
it reaches the proton, it would then turn back 
and move away from the proton. The subsequent 
motion is then described by Eq. (9). We do not 
discuss here the motion in detail, but merely 
give a summary of the results. It is easy to show 
that there is no solution with y>«, x0. The 
only possible types of solution are those in which 
either the electron is brought to rest while 
moving away from the proton, or in which the 
electron escapes to infinity while the velocity 
ultimately decreases and tends to a finite value. 
In the first case, the approximate solution near 
the point where the electron is brought to rest 
is similar in form to (15). The subsequent 
motion is such that the electron steadily ap- 
proaches the proton with increasing velocity. In 
the second case, the solution for large x is given 
approximately by 


pl 
y=C+--(1+8), (17) 
Ax 


where C2 1 and 6-0 as x &. 

Under given initial conditions the final motion 
consists of either:—(i) the electron moves 
towards the proton with increasing velocity and 
collides with it, the approximate solution being 
(14); or (ii) the electron recedes away from the 
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proton with steadily diminishing velocity, the 
approximate solution being of the form (17), An 
interesting feature in this result is that in both 
motions the final acceleration is directed towards 


the proton, which is in accord with the ideas of — 


Newtonian mechanics where the acceleration jg 
always in the direction of the force, which, in the 
problem under investigation, is towards the 
proton. 


IV. TWO-DIMENSIONAL MOTION 


Now we shall consider the more general 
motion of an electron which is spiralling round a 
fixed proton. The general motion will be three. 
dimensional, since the initial position, velocity 
and acceleration can all be arbitrary. Under 
certain initial conditions, however, the motion 
will be two-dimensional, and it is these motions 
that we shall consider here, as these are the only 
ones which are likely to be physically permissible, 
We shall chiefly be concerned here with the 
existence of solutions in which the electron 
spirals round and falls into the nucleus. Such 
motions were found to be not allowable by the 


Lorentz-Dirac equations® and it is of interest to , 


determine whether the modified Eqs. (2) lead to 


different results. Taking axes in the plane of . 


motion and the origin at the position of the 
proton, the equations of motion are 


dx . dx/@t dx? dy? 3 dt 
cecal are se ty. 2, 
ds? ds* ds ds? ds? ds? 2 ds 
> (18) 
d*y (2k-+1) ao dt d*x? —)|- 3 dt , 
a-— — —+—f —-—-— -—-— _ } } = —-- yr 
ds? ds? dst dst ds? 2ds° J 


where 


We transform to V and y, where 


Toa 


is the velocity, and r—y is the angle between the radius vector and the velocity, so that 


dx 
—=V cos (r+0-—y), 
ds 


We change the independent variable from s to r where r =e~’ and let dashes now denote differentiation 


7 dx? 
ds 











dy*\? 
ds 







‘dy 
—=Vsin (r+0-y). 
dt 










18) 


on 
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with respect to r. The equations of motion then give 


V"+V'+ cs Vy’ tany— V(ii+ V)(y o»| ot See. (19) 
” 7 la Um hd Pe ‘— tan -—— = - 
k+)| V(1+ V2) Vcosp 2 V?cosp 


d 3V’ 
cb+1)|—(W cosy —siny) + ( +—) (Y cosy—siny)| 
T 
ae~* 3(1+V%)! 
(y’ cosy —siny) =—- —————tany. (20) 
V cosy 2 


ys 





‘ We seek those solutions according to which the electron spirals round and ultimately falls into the 


nucleus. Assuming that this happens, let us consider the behavior of the solution as r—0, that is, as 
TIO. 

We follow the method used in an earlier paper and easily see that except possibly for the case 
when V-—> and ¥y—0 as r— ~~, there is no solution in which the electron falls into the nucleus. We 
shall therefore concentrate our attention on the possible solution in which V-+>” and y—0. For 
large r, Eq. (20) is approximately 

VW’ -V+(14+3V'/V)(W —y) =0, 


which has solution 


VY —y=Ce"/V*, (21) 
where C is an arbitrary constant of integration. Substituting in (19) and approximating we obtain 
(2k+1){ V"+V'+ V"/ V3} = —3/2V. (22) 


This equation can be shown to be equivalent to Eq. (10) which has already been considered above. 
Hence we have as an approximate solution of (22) V=(2ur)!. Substituting in (21), we obtain 


v —y=Ce~*/(2ur)!. 





Hence 
? e72 
ye"= C —dr (23) 
« (ur)! 
Cer 1 3 1 3-5 1 
= _ —S =. | (24) 
pe!2 | 273/292 5/2 28 TI? 
Hence for a first approximation we take 
y=Ae"r, (25) 


where A is a positive constant of integration. Therefore 


rd0/dr = —tan y= —Ae~*r-, for large 7; 
and hence 
d0/dr= Aer. 
Therefore 


9-4 fer tar=Aer| —— +--— — — -—- + 


Hence the equation of the curve is given approximately, for small r, by 
6=Ar (logi/r)-1. (26) 


Thus, the equations of motion (2) allow the electron to spiral round and fall into the nucleus, whereas 
the Lorentz-Dirac equations do not permit such a collision. 
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Covariant conditions for rigid-body motion are set up. They are equivalent to those pro- 
posed by Born and lead to the linear speed-distance law for the case of rotation about an axis. 
This result and also a discussion of the transformation equations in going over to a rotating 
frame of reference are used as arguments for the desirability of retaining the concept of rigid 
rotation with the linear speed-distance law, contrary to the opinion expressed by Hill. The 
rotation question is also considered from the standpoint of angular velocity, and one is lead 
rather naturally to two fluid velocity distributions, one of which was found by Hill. The 
expression for the spatial distance between two points on a rotating disk obtained by Berenda 
is derived without the use of the assumption introduced by the latter. 





HE question of whether the geometry on 
a rotating disk is Euclidean or not has 
interested many writers. Berenda! has subjected 
the problem to a careful analysis. However, 
there are still a number of questions connected 
with the problem of rotation in relativity theory 
that require further investigation, and some of 
them will be considered in the present work. 


1. MOTION OF A RIGID BODY 


In discussing the rotating disk, one frequently 
begins with an inertial frame of reference de- 
scribed by a cylindrical polar coordinate system 
with coordinates 70, 9, 20, to, and one then goes 
over to a rotating frame of reference attached to 
the rotating body, with coordinates r, 0, 2, ?, 
given by 


To=7, O=O0+t, 2=2, to=t, (1) 


where w, the angular velocity about the z axis, 
is constant. Starting with the usual expression 
for the line element in the non-rotating coordinate 
system, one gets in the rotating system: 


ds? = — (dr?+r'd#+dz?+ 2wr*d6dt) 
+(1—w*r?)d#, (2) 


with units in which the velocity of light in empty 
space is unity. 

Hill? pointed out (as von Laue* had done 
earlier) that a body rotating in this way would 
have a limit to the radius it could have, for the 

1C. W. Berenda, Phys. Rev. 62, 280 (1942). Other 
references will be found in this pew. 

*E. L. Hill, Phys. Rev. 69, 488 (1946). 


3 M. von Laue, Relativitatstheorie (Vieweg, Braunschweig, 
1921), Vol. 2, p. 24. 
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velocity would vary linearly with the radius and 


would exceed that of light for r>1/w. For this 
reason, Hill was led to conclude that the speed- 
distance law must be non-linear. 

While one can conceive of a fluid rotating with 
an arbitrary speed-distance law, the case of 
greatest interest is that of the speed-distance 
law for a rigid body. To investigate this, one 
must first have some criterion for rigid-body 
motion. Such a criterion was first given in rela- 
tivity theory by Born.‘ 

In classical physics one can characterize the 
motion of a rigid body by the fact that the rate 
of strain vanishes. In a Cartesian coordinate 
system, the velocity components satisfy the re- 
lation 


Ov; /dx,+00,/dx;=0 (4, k= 1, 2, 3). (3) 


In a relativistic treatment one looks for a co- 
variant equation which reduces to (3) in a 
Galilean system if the velocity is small (compared 
to that of light). The obvious generalization is to 
introduce, in an arbitrary coordinate system, the 
symmetrical tensor 


Pes = 4 (ty:»t+Uy:y)» (4) 
where the velocity 4-vector u* is given by 
uw =dx/ds (5) 


and a semi-colon denotes covariant differentia- 
tion, and to take as the condition for rigid-body 


4M. Born, Ann. d. Physik 30, 1 (1909). See also W. 


Pauli, Jr., Encyklopidie der math. Wissenschaften, Vol. 5, 
Part 2, pp. 689-691. 
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tion 
- P,,=0. (6) 


However, this condition represents too severe 
a restriction. In view of the fact that the vector 
u> satisfies the identity 


wu, =1, (7) 
it follows from (6) on multiplication by u’ that 
Uy;»t” = 0. (8) 
This expression is just the covariant form of the 
acceleration vector 
6u*/ds =du/ds+ { as} usu, (9) 
so that every particle of the body moves along a 
geodesic, or in gravitation-free space, along a 
straight line. 
It is therefore necessary to weaken the condi- 


tion imposed on the motion. For this purpose 
we replace py» by “ 


Pur =} (Upset Uo:p — Up; aU "Uy — Uy; qlt*U,), (10) 


since we then have the identity 


Pu? =0. (11) 

Let us now take as the condition for rigid-body 
motion 

Pur = 0. (12) 


In a Galilean coordinate system, at a point 
where the velocity 3-vector vanishes, this re- 
duces to the classical condition (3). Therefore 
it is equivalent to the condition proposed by 
Born.‘ 

Since Eq. (12) is a tensor equation it can be 
applied in any frame of reference. If we take an 
inertial system with Cartesian coordinates then 
we can describe rotation about an axis by 
setting 


uw=—ul=oy, U=—u?=—ox, 
; Us=—v=0, ug=ut=(1+0'r?)!, (13) 


with 
c=oa(r), P=x?+y?’. (14) 
From Eg. (12) one gets the single equation 
da/dr=a'r (15) 


which has for its solution 


g=w/(1—w*r’)!, (16) 
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where w is a constant. Going over to the three- 
dimensional velocity v by the relations 


u = (u;*+us*)! = or =v/(1—v*)!, (17) 
we get 
v=wr. (18) 


This result was obtained by Herglotz® by means 
of a different procedure. 

That it is possible to have a covariant law of 
rigid-body motion (even if it has only a small 
number of solutions) appears to be of importance 
in connection with the inner consistency of the 
theory of relativity. Since the fundamental con- 
cepts of the theory include those of the rigid 
frame of reference and the rigid measuring rod, 
it would be unsatisfactory if no rigid bodies or 
rigid-body motions existed in the theory. 

In 1911, von Laue® presented an argument of 
a general nature against the existence of a rigid 
body in relativity theory. It is based on the fact 
that a rigid body is-expected to have a finite 
number of degrees of freedom, while, on the 
other hand, one can set up N disturbances near N 
separated points of the body, and they will be 
non-interacting, i.e., independent, for a suffi- 
ciently short interval of time (because of the 
finite propagation time required by relativity 
theory) so that there are at least N degrees of 
freedom, where N can be increased indefinitely 
for a continuous medium. However, this argu- 
ment assumes the possibility of setting up arbi- 
trary deformations in the body initially. Since 
however the initial deformations can only be 
obtained by displacements with velocity com- 
ponents governed by Eq. (12) it follows that it 
is not possible to have arbitrary initial dis- 
turbances, but only those which can be arrived 
at by means of motions described by Eq. (12), 
i.e., displacements which can be obtained by 
integrating this equation, starting from an un- 
deformed state. 

It might be remarked here that the difficulty 
pointed out by Hill does not appear to provide 
a sufficiently strong argument for giving up the 
linear speed-distance law. It is true that no rigid 
body would have a radius equal to, or exceeding, 
1/w, but there seems to be no reason why an 


5G. Herglotz, Ann. d. Physik 31, 393 (1910). 
*M. von Laue, Relativitdtstheorie (1921), fourth edition, 
Vol. 1, pp. 203-204. 















idealized rigid body with a smaller radius could 
not rotate according to the linear law. That it is 
desirable to retain the concept of the rigid 
rotation can be seen from the following con- 
siderations: 

Let us go back to Eqs. (1) and (2), and let us 
look for the most general transformation that 
can be introduced in place of Eq. (1) which 
(a) will be linear in ¢ and @, (b) will have axial 
and cylindrical symmetry, (c) will be single- 
valued, and (d) will give a static line element. 
The conditions (a), (b), and (c) lead to the 
transformation equations 


ro=e(r)+f(r)t, G=0+g(r)t+h(r), 
Zo=2, to=p(r)itg(r). (19) 


Applying the condition (d), one finds that 
f=g'=p'=0, (20) 


where a prime denotes a derivative with respect 
to r. What is of interest to us here is that unless 
g’=0 (@—6» a linear function of ¢ with a constant 
coefficient), it will not be possible to have a 
static metric. This condition is, of course, just 
the linear relation between speed and distance 
from the axis. 


2. ANGULAR VELOCITY IN RELATIVITY THEORY 


Although it appears that the motion of a 
rigid body is best discussed by means of the 
criterion of the previous section, it is neverthe- 
less interesting to examine the rotation problem 
by a consideration of the angular velocity dis- 
tribution. This is the general method that was 
used by Hill.? 

In classical physics the angular velocity is 
defined as the vector 


o=4}VXv. (21) 


The obvious relativistic generalization consists 
in introducing an antisymmetric tensor 


Qu» = 3 (0u,,/Ox” — Ou, /dx") (22) 


in an arbitrary four-dimensional coordinate 
system. 

In an inertial system with Cartesian coordi- 
nates, if we have rotation with the components 
u* given by Eqs. (13) and (14), the following 
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components of Q,, are obtained : 


Q:2=0+40'r, (23) . 
with Geer, Gyo sy, (24) 
F=}(0?+00'r)/(1+0*")!, (25) 
If we now set 
Q32=w =constant, (26) 


then from (23) we find for the singularity-free 
solution’ 
c=. (27) 
This gives : 
u=wr, v=wr/(1+w*r")!. (28) 
The components wi, and wes are given by Eq, 
(24), with 
F=}w*/(1+w"r*)}, (29) 


so that, for small speeds, they are proportional 
to the centripetal acceleration components. 

Now, it is not satisfactory for the fluid motion 
to be determined by Eq. (26) since the latter 
is a condition on only one component and is 
not covariant. One can improve matters by 
introducing 


S= Pu, =v ;,u" = 5u/ds, (30) 
the acceleration vector. In the present case one 
finds 

Si=—}jo'%x, S=—jor'y, S=St=0, (31) 

so that the condition ~ 
S*;4=0 (32) 
would give the same solution as above. Another 


possibility, not as simple, is to take as the 
condition for the motion, 


oS=0. (33) 


If, instead of the above coordinate system, we 


follow Hill? in introducing a Galilean system — 


(coordinates x’, y’, 2’, t’) with a given point as 
origin and moving so that the fluid there is 
momentarily at rest with respect to it (u:’ =’ 
=u; =0, us’=1), we get for the angular velocity 
components 


Qi2’ = (0 +40'r+4o°r?)/(1+0°r")', (34) 
214’ = 4o"r, (35) 
7 Nathan Rosen, Phys. Rev. 70, 93 (1946). 
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provided we take the X’ axis in the radial direc- 
tion and the Z’ axis parallel to the Z axis. If 


one sets 
(36) 


242’ => 
and introduces a new variable 


f=0/r, 
o=t/(1—£7")!, 


t’r+2§=20(1—#r*), 
which is essentially the equation obtained by 
Hill. 

It is possible to put the condition expressed 
by Eq. (36) into covariant form. For this purpose 
let us introduce what might be called the “spatial 
angular velocity tensor” 


(37) 
that 
_ (38) 


one gets 
(39) 


Ops = Que _ Q, all“, + QU "Uy,, (40) 


which satisfies the identity 
(41) 


Wy» =0. 


It might be noted that the vanishing of Q,, leads 
to the absence of all acceleration; that is not 
the case with w,,. Corresponding to Eqs. (34) 
and (35) one finds as the only non-vanishing 
components 

(42) 


_ — , 
@12 = —wo1 =Qy2. 


Hence, an invariant condition corresponding to 
Eq. (36) can be written® 


Wy rw” = 2”, (43) 


We see that, depending on whether one works 
in the fixed coordinate system or the system in- 
stantaneously moving with the fluid, one arrives 
naturally at one or the other of two velocity 
distributions, differing from that for the rigid 
body as discussed in the previous section. Other 
distributions are also possible, of course. How- 
ever, at this time there does not appear to be 
any strong reason for giving up the linear speed- 
distance law in favor of any other distribution. 


*In reference 7 the author questioned Hill's procedure 
because it made use of a uniformly moving coordinate 
system in describing the accelerated motion of the fluid, 
so that the acceleration components of the fluid were not 
taken into account. However, we see now that this is not 
an objection if the description is given in terms of wy, 
since the acceleration components in this case are not 
involved. 
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3. SPATIAL DISTANCE IN A ROTATING SYSTEM 


Returning now to the rigid-body rotation as 
discussed in Section 1, let us consider a question 
treated by Berenda in his investigation of the 
nature of the geometry on the rotating disk. 
The question was this: given a frame of reference 
with the line element of Eq. (2), what is the 
correct expression for the distance dl between 
two neighboring points at rest in this reference 
frame? Berenda' obtained this distance as the 
length of a four-vector which had for its space 
components the three spatial coordinate differ- 
entials, but also had for its fourth component a 
value chosen to make this vector orthogonal to 
the world-lines of the two points. The identifica- 
tion of “rigid metric rod readings on the disk 
with spatial dimensions and geometry derived 
from such orthogonality’’ Berenda considered a 
matter of assumption. However since the ex- 
pression for di? derived from this assumption 
determines the spatial geometry of the disk, it 
seems desirable to derive it by another procedure 
free from such an assumption. 

It is well known that in the general relativity 
theory the coordinates used in locating events in 
space-time need have no simple relations to the 
readings of measuring rods and clocks. What 
requires stressing here is the fact that even the 
value of the interval or line element need have 
no simple relation to the readings of scales and 
clocks in a given frame of reference. This is 
associated with the fact that the readings of 
scales and clocks will be influenced by gravita- 
tional and inertial forces. Only in a Galilean 
frame of reference can intervals be obtained 
simply from clock and scale readings, according 
to the ideas of special relativity. Hence, con- 
sidering the two points of the body in the 
neighborhood of some particular values of r, 8, z, ¢, 
let us introduce a local Galilean coordinate 
system momentarily at rest relative to the 
rotating system by the relations 


dy’=Ad0, dz'=dz, 
dt' = Bdi+-Cdé, 


dx’ =dr, 
(44) 


where A, B, and C are to be determined so as 


to make 


ds* = —dx"*—dy" —dz"+dt”. (45) 
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One readily finds that 


A=r/(1—w*’)!, B=(1—w*r’)}, 
= —wr?/(1—w*r)!, (46) 


This is essentially the same system as the one 
used by Hill.2 From the form of Eq. (45) it 
follows that 


dP = dx"*+-dy"+dz”, (47) 


G. H. LIVENS 





or, on the basis of the transformation equations, 
dP =dP+r(1—wr)de+dz?, (4g) 


which is the result obtained by Berenda ang 
shows that the spatial geometry on the surface 
of a rotating disk ( = constant) is non-Euclidean, 
It is also interesting to note that the time 
interval between two events dt’ is given by 


dt’ = (1—w*r*) dt —wr?(1—w*r*)-1d0. (49) 
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Note on Magnetic Energy 


G. H. LIvens 
Department of Mathematics, University College of South Wales and Monmouthshire, Cathays Park, Cardiff, Wales 


1. 


N a recent paper! under this same title E. A. 
Guggenheim discussés certain results from a 
paper? published by me recently and claims that 
they are special cases of more general ones ob- 
tained by him on a previous occasion.’ As the 
results quoted from my paper are themselves 
special cases of quite general ones contained 
therein, and his note generally misrepresents the 
scope and purpose of my work, perhaps I may be 
permitted a reiteration of the main outlines of my 
arguments and the opportunity to discuss the 
relation they bear to the methods employed by 
him. 
The following symbols are used: 
B magnetic induction, 
H magnetic force intensity, 
I; induced magnetization intensity, 
I, permanent magnetization intensity, 
4, linear electric current, 
N, magnetic flux through #,, 
u permeability, 


£& Lagrangian function, 
3 Hamiltonian function. 


Guggenheim bases his discussion not on 


1E. A. Guggenheim, Phys. Rev. 68, 273 (1945). 
*G. H. Livens, Phil. Mag. 36, 1 (1945).-Cf. also Proc. 
Roy. Soc. A93, 200 (1916), and Phil. Trans. Roy. Soc. 
A220, 207 (1919). 
* E. A. Guggenheim, Proc. Roy. Soc. A155, 49 (1936). 
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Maxwell’s theory in its original form but on the 
modification given to it by Cohn.‘ The essence of 
this form of the theory, like that proposed by 
Hertz which it follows closely, is that it in- 
corporates the induced polarizations and the 
aether, whatever this may be, into a single 
transmitting medium whose elastic quality is 
summed up in the characteristic constant, the 
permeability y. This hypothesis of a single 
medium, excluding as it does the possibility of a 
displacement of the polarized medium from one 
position of the field to another, proves however to 
be a fatal handicap in a theory which has eventu- 
ally to be extended to cover electromagnetic 
phenomena in moving media. And it was pre- 
cisely for this reason that Larmor and Lorentz 
were forced back to the views held by Kelvin and 
Maxwell that the only really satisfactory treat- 
ment of these affairs interprets them in terms of a 
separate universal transmitting medium with its 
own stress on which is superposed the polarized 
media with their reacting mechanical forces. This 
implies that it is absolutely essential to distin- 
guish between the parts of the field vectors and 
energy which belong to the aether and remain 
with it and the parts which belong to the matter 


*E. Cohn, Das elektromagnetische Feld (1900). My 
knowledge of this book is derived from this earlier edition. 
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and are carried along with it in its motion through 
the field.® 

However, so far as the present discussion of the 
quasi-static energy relations of the field are con- 
cerned there is no essential difference between the 
different forms of the theory, except in the matter 
of interpretation. The differences, however, begin 
to appear when these energy formulae are applied 
to determine the forces on the magnetic media by 
applying the principle of virtual work, and it has 
long been known that the theory in the Helmholtz- 
Hertz-Cohn form leads, when properly inter- 
preted, to quite untenable results.*® 


2. 


My own discussions in this subject are based 
on the generalizations of Maxwell’s theory formu- 
lated by Larmor and Lorentz, and so far as the 
mechanical relations of the magnetic field are 
concerned they start from an assumption which 
is really the basis of all electrodynamic theories 
from Maxwell onwards. This, in its most general 
form, asserts that the magnetic part &,, of the 


Lagrangian function of a system consisting of a 


series of linear currents 7, enclosing fluxes NV, and 
a distribution of permanent polarity J, in the 
presence of magnetizable media distributed over 
the field in any manner is such that 


5 Ln = Ee Nidig + f (Fél,)dv 
: (1) 
~— f (B&H)dv+ f (Fél,)dv- +. 


This formula is interpreted for the moment in 
terms of a theory which regards F as the intensity 
of the magnetic force acting on the polar ele- 
ments. The usual assumption is that F=H, but 
there are grounds for thinking that F=B is a 
better choice. The integrands of all space inte- 
grals are understood to relate simply to the local 
conditions at a point, and their variation pro- 


5 This conclusion should be compared with Guggenheim's 
statement (in Proc. Roy. Soc. A155, 99 (1936)) that “as 
the ether and the matter occupy the same space, there is 
no means of distinguishing between the energy of the 
ether and the energy of the matter. We can therefore 
attach no physical meaning to this division of the energy 
into two parts.” 

*G. H. Livens, Phil. Mag. 32, 162 (1916). 
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duced in any process of building up the system 
from infinite dispersion of its elements. 

Formula (1) is of the utmost possible generality, 
especially in this virtual differential form, and it 
defines the complete electrodynamic relations of 
the quasi-stationary system, and this includes 


-not only the relations of the currents and 


permanent magnets, but also those of the induced 
magnets as well. There is no implied necessity for 
any relations between the vectors involved in it, 
anything in the nature of a law of induction, that 
would make it the complete differential of a 
unique function 2&,, although naturally the 
greatest’ mathematical interest attaches to the 
cases when such relations exist. 

Guggenheim starts not from the Lagrangian 
function, but from the associated Hamiltonian or 
available (conserved) energy function which he 
writes in the integrated form 


n= f " udN., 


and then, assuming with Maxwell, that the 
velocities in the electrical coordinates are i, he 
derives. the corresponding Lagrangian function 
£» in the form 


tis 
>> f Nei. 
N.=0 

In both of these formulae the summations are 
taken not only over the ordinary currents of the 
system, but also over a whole series of other 
currents introduced to mask the effects of the 
permanent magnetism in the early stages of the 
build-up process represented in the integrals. 


3. 


The choice by Guggenheim of the Hamiltonian 
function as the more fundamental is based on the 
possibility of deducing it directly from Maxwell's 
equations. The argument used by him for this 
purpose depends however essentially on the con- 
cept of ‘‘zero states,”’ and to define these. he uses 
the following postulate, I quote his own words :? 

“If, however, each permanent magnet is sur- 
rounded by suitable circuits it is theoretically 
possible for currents to flow in these circuits such 


7 E. A. Guggenheim, reference 3, page 55. 
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that magnetic induction due to them is every- 
where equal and opposite to the magnetic induc- 
tion due to the permanent magnetization.” 

This, of course, is true only when the perma- 
nent magnetism is distributed uniformly through- 
out the volume of each such magnet. If we try to 
solve Maxwell's equations for a steady magnetic 
system with currents on the surfaces of the 
permanent magnets and with B=H+4rI=0 
everywhere, we soon realize that, in the general 
case, it is necessary to supplement the surface 
currents by a distribution of volume currents 
throughout each magnet with a density at each 
place proportional to curl J, and these-volume 
currents vanish only when J is constant through- 
out each magnet. In the general case the two sets 
of currents, surface and volume, constitute a 
distribution which is then the exact opposite of 
that which is the effective equivalent of the 
electron current whirls ultimately constituting 
the magnetism. The two systems would thus 
cancel out completely in every detail except as 
regards their magnetic forces, the definition of 
which for a current distribution being somewhat 
different to that for the equivalent polarization.*® 
Properly interpreted therefore Guggenheim’s 
“zero states” are identical with the usual concept 
of the zero fields of the infinitely dispersed system 
in the more conventional building up procedure, 
except that there is in the interior of the magnets 
a residual magnetic force Ho’ resulting from the 
different interpretations of this vector for the two 
cancelling systems. 

The difficulty is, of course, the introduction of 
permanent magnetism and its associated energy 
into the simple form of the energy equation 
derived from Maxwell’s theory. But if we are 
careful to interpret the energy equation properly 
this difficulty is soon overcome. The equation 
itself determines a formula for W,,, the part of the 
field energy associated with the magnetism, which 


® Cf. Larmor, Aether and Matter, p. 106, “‘On the specifi- 
cation of a magnetic distribution in terms of a continuous 
distribution of currents.’’ The same idea probably underlies 
Guggenheim’s statement (Proc. Roy. Soc. A88, 100 (1936)) 
that “if we could imagine all matter removed leaving 
behind the molecular currents it is B, not H, which would 
remain unaltered.’’ Actually, of course, it is the change in 
the specification of the magnetization as a distribution of 
poles to a distribution of currents which necessitates the 
change in H, and the presence, or otherwise, of the mag- 
netically inert parts of the matter has nothing whatever 
to do with it. 


is such that 
1 
ibW,,=— f (H6éB)dv. 
4n 


But this formula represents the magnetic part of 
the corresponding Hamiltonian or conserved 
(available) energy only in so far as the magnetic 
connections of the system controlling the build. 
up of the polarizations and currents involve no 
dissipation other than that which occurs in the 
linear current circuits, and allowed for directly in 
the complete equation. But it is of the nature ofa 
permanent magnet that any energy expended by 
the field in separating its poles in the general 
build-up of the field is immediately locked up in 
the rigid connections holding the magnetism 
permanent, so that it can no longer be reckoned 
as part of the available energy. In this case, 
therefore, the Hamiltonian or conserved energy 
function is diminished from the value W,, by the 
amount expended by the field in building up the 
permanent magnetism. We have therefore 


1 
550, =— f (H6éB)dv— f (FéI,)dv, (2) 
4 ; 
and as we have assumed a kinetic origin for the 
linear currents only this corresponds to a 
Lagrangian function £,, which is such that 


1 
bon =— f (BsH)dv-+ f (FSI,)dv, 


the formula (1) above. 


4. 


Guggenheim also converts his general formulae 
for £, and 3,, into space integrals over the field 
obtaining results which are equivalent to 


1 H 
Ln =— fac (BdH’), 
4n Ho’ 


1 
Kn =— fof H'dB, 
4r 0 


Hy)! being the value of H’, more clearly defined 
below, in the zero field when H=B=0 every- 
where. These are apparently very different 
formulae from (1) and (2) above. But we must 
remember that in Guggenheim’s treatment there 
are other currents involved in the build-up, so the 


(3) 








3) 





magnetic field intensity which occurs in the 
integrands, and now called H ’, is a vector which 
is circuital with respect not only to the specified 
currents of the system, but also to all these 
additional currents as well. In other words it is a 
very different vector from the H employed in the 
usual treatment and used in our (1) and (2) above, 
although of course the upper limit of the two 
yectors, taken when the masking currents have 
disappeared, will be the same. The relation be- 
tween the two vectors at any stage of the build- 
up is easily obtained if Guggenheim’s method is 
properly interpreted to include the volume cur- 
rents as well as the surface ones, and it shows that 
if H is the usual vector which is circuital with 
respect only to the original currents of the system 
‘at this stage then the H’ of Guggenheim’s 
formulae is 
H’=H—4rTl, 


if the currents in this transition stage are just 
sufficient to cancel a distribution of magnetic 
polarity of intensity J.. The formula for £,, in (3) 
is therefore in the ordinary notation 


tan fof (Bax) — f dv F (BdI.), (4) 


where I.9 is the value of J, which completely 
cancels J, together with any induced mag- 
netization J,» in the zero state, 


I= pt lio. 


Substituting then J,=J,—I., dI,=—dI,, in 
the second integral, it reduces to 


1 pl Ip 
£.=— f dv f BdH+ f dv f (BdI,) 
4r 0 —Tio 


which is the equivalent of our (1) with F=B 
apart from a trivial constant which in fact 
measures the quasi-elastic energy associated with 
any induced magnetization which may be present 
in the zero field. 

We have also 


1 
-—faof (as 
Run ~ fe 4 ( d ) 


-— f (HB) —— fac J ; (BdH’) 
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~~ fof uaB)— f dof (Bat, 


which is the equivalent of formula (2). 
Guggenheim considers two special cases of 
these formulae corresponding to two possible 
laws of induction. As H is zero in the initial field 
we have 
Hy’ = —4rlo= —4r(I,+T io). 


If then we write J 4o=(u—1)Ho’/4r, we have 
Hy = —4rI,/pu, Tio=—([(u—1)/u lp; 
while if Io =(u’ —1)Bo/4r=0, we have 
Hy’ = —4rT;,. 


In this latter case the formulae (3) are actually 
identical with (1) and (2) interpreted of course 
with F=B, the lower limit of the second integral 
being zero. 


5. 


Guggenheim seems quite unaware of the differ- 
ence between the H’ of his formulae—he, of 
course, calls it H—and the more conventional H 
employed by Cohn and myself, for he criticizes 
Cohn’s choice of a zero lower limit for his H, 
when this is precisely the value which corre- 
sponds to his own choice of —42J, or —4rI,/u as 
the lower limit of H’. 

The real difference between Cohn’s formulae 
and those given by Guggenheim lies really in 
another direction altogether, and not, as surmised 
by Guggenheim, in the choice of a lower limit for 
the integrals involved. Cohn’s whole treatment is 
in fact interpreted not in terms of the usual 
induction vector B, as implied by Guggenheim’s 
quotations, but of another vector M, variously 
called by Hertz the polarization or magnetization 
of the field, and which is related to B by the 
equation 

M=B-—4r!I;,. 


This makes Cohn’s formula for %,,, viz., 


1 p™ 
Rn =— HdM 


4n v5 


identical with our (2) with F=H. Guggenheim, 
thinking in terms of his own treatment seems 
inclined to write off the second integrals in dJ, in 
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all our formulae as having zero values when the 
configuration is kept constant. As he says in a 
letter to me “I need hardly point out that at 
constant configuration as J, is constant 


f HdI,=0,” 


and this generally speaking seems to. form the 
basis of his claim to have generalized the results 
given by Cohn. But when properly interpreted by 
the conventional build-up process employed by 
Cohn and myself these integrals do not neces- 
sarily vanish. They will vanish in Guggenheim’s 
treatment but only because an equivalent non- 
vanishing integral 
f Hal. 


is included with them. 

It is this unwillingness to recognize the con- 
ventional build-up procedure used by all previous 
writers, as a legitimate method of evaluating, at 
constant configuration, integrals of the kind here 
under review which is the cause of Guggenheim’s 
difficulties with previous treatments. And as he 
interprets ‘‘constant configuration” to include 
also a constant specification of the distribution of 
permanent magnetism he is obliged to employ the 
device of covering the magnetism with a current 
shield or mask, which has then to be removed 
gradually during the build-up of the finite cur- 
rents which form the rest of the generating 
system. As we now see the two processes, re- 
moving the mask and building up the magnetism 
from zero are in effect almost the exact, equiva- 
lents of each other both mathematically and 
physically and they lead to identical formulae for 
the available energy in the field in every case. 
One method cannot therefore be described as any 
more general than the other. 


6. 


The rest of the argument of my paper was con- 
cerned with an attempt to derive the simplest 
and most natural form of the expression for the 
general Lagrangian function, or rather of its 
virtual variation, expressed as an integral over 
the field and in a form which exhibits its non- 
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integrability when it occurs, as arising simply and 


.solely in the hysteresis effects associated then 


with the induction of the polarization. 
In the usual case when F=H a familiar argu- 
ment soon shows that 


‘ . 
ion —5 f {2(BH) —H®}dv— f (H8I,)do: 
Tv 
while when F=B the same argument gives 
1 
i2_=—bf Bido— { (Baldo. 
8r 


These expressions are independent of any par- 
ticular law of induction, or even of whether a law 
exists at all, and they apply to every type of 
medium, isotropic or crystalline. In each of them 
the second integral measures the energy stored in 
the quasi-elastic connections holding the induced 
magnetism against the action of the field. The 
first integral, represented in both cases by a 
complete differential, thus represents the purely 
magnetic part. of the function, the part not 
definitely located in the latter. And the fact that 
in the second case this purely magnetic part of 
the energy assumes a very simple natural form in 
terms of the single vector B provides one of the 
arguments in favor of the choice of this vector as 
the fundamental force vector of the theory. 

This was the main conclusion of my paper, but 
being anxious naturally to show the relation be- 
tween these general formulae and those given in 
the more familiar treatments, which invariably 
use linear laws of induction, a good deal of 
attention was devoted to the special form of the 
results when these linear laws are assumed. It 
was then shown that the formula for £,, which, 
in the case when F=H, becomes 


£.=— f cars) -nrr\av 


is identical with that usually employed in the 
dynamical theory of currents, and that, while it 
apparently differs from the more familiar form 
when the field arises from permanent magnets it 
still, even in such cases, gives a total for the whole 
field which is identical with that given by the 
more usual formula. 
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When the law of induction is linear there is, of 
course, also an energy or Hamiltonian function, 
the special case of formula (2). This is shown to be 
represented in all cases by the familiar result 


1 
=— | »pH*dv 
8r 


so the theory provides a self-consistent set of 
simple formulae for determining all the details of 
the mechanical and electrical behavior of the 
system in all cases, whether the law of induction 
is linear or not. 

Under the impression, as he says, that I have 
confined myself to what I call linear laws of 
induction, it is these special results from my 
paper which Guggenheim claims as particular 
examples of his own general formulae. The dis- 
covery as we now see is not a very surprising one, 
seeing that one form of the general formulae of 
my paper, to which, however, he makes no ref- 
erence, is in fact identical with the general 
formulae given by him, when these are properly 
interpreted, however different they may be in 
appearance. 


7. 


There are other imperfections and misunder- 
standings involved in the rest of Guggenheim’s 
work which will be dealt with in detail in a more 
comprehensive survey of the whole subject to be 
published elsewhere. Enough has now been said 
to justify the conclusion that his claim to have 
produced more general formulae which include 
those obtained by Cohn and myself as special 
cases is entirely without justification. His for- 
mulae when properly interpreted prove to be 
identical with formulae given by both of us many 
years previously so that they are certainly no 
more general than the results given by Cohn and, 
in so far as his discussion is framed on the re- 
stricted basis of Cohn’s theory, his results have a 
far less degree of generality than the equivalent 
ones obtained by myself on the much wider basis 
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of the Kelvin-Maxwell theory. The only ad- 
vantage of Guggenheim’s formulae is in their 
mathematical compactness, but, as his subse- 
quent discussion shows, this may have its dangers 
when attempts are made to interpret the relevant 
parts of them for dynamical purposes. 

There are two other small points raised in 
Guggenheim’s note which can be disposed of 
quickly before closing this already over-long 
note. 

The question of the compatibility of the two 
interpretations of the linear law of induction 
given in my paper, and discussed at some length 
by him, is entirely irrelevant. There is really no 
conceivable reason why two such alternative 
physical laws should be compatible among them- 
selves, but it is hoped that eventually one of them 
will prove to be more consistent with experience 
than the other. 

Then at the end of the note Guggenheim 
complains that I use the term “fundamental 
(aethereal) force vector’’ without explaining 
what I mean. In following, however, the example 
of Cohn and ignoring the physical processes 
involved in the induction of the polarizations, he 
thereby overlooks one of the most important 
functions of the magnetic force vector. The 
forces on the currents are determined in terms of 
the induction vector B and, as the same vector 
can also, without ambiguity, be used to determine 
the forces on the permanent magnets, it is the 
only field vector which functions as a force 
intensity in this restricted theory. In the more 
detailed theory, however, a choice has to be made 
for the force vector which is effective in the 
polarization process and tradition, following the 
electrostatic analogy, has chosen H for this 
purpose. My contention is that, if B is chosen for 
this purpose as well, so that it becomes the one 
and only fundamental vector which functions as 
a force intensity then the whole of the energy and 
mechanical relations of the field, in the only 
really consistent form so far framed, assume their 
simplest and most natural form. 
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Rotational Spectra of Some Linear Molecules 
near l-cm Wave-Length 


C. H. Townes, A. N. HoLpEN, AND F. R. MERRITT 
Bell Telephone Laboratories Inc., Murray Hill, New Jersey 
November 20, 1946 


OTATIONAL spectra of several linear molecules have 

been found in the region near 1-cm wave-length. The 

lines are observed on an oscilloscope on which frequency is 

plotted against variation in absorption. Absorption takes 

place in a 4-meter length of wave guide containing a few 

tenths of one mm Hg pressure of the absorbing gas. The 
technique is similar to that already described.' 

In some cases there are a number of weaker lines ac- 
companying the main rotational line. The measured fre- 
quencies and relative intensities are given in Table I. 
For the stronger lines frequency measurements have a 
probable error of about two megacycles, although frequency 
differences are somewhat more accurate. Relative intensi- 
ties are not highly accurate, particularly for the weakest 
lines. The CICN lines at 23,885 and 23,389 mc appear to 
be caused by Cl*®CN and CI7CN molecules, respectively, 
which are in the ground vibrational state and undergo a 
transition from J=1 to J=2. Similarly the BrCN lines at 


TABLE I, Measured frequencies and intensities of lines in 
rotational spectra near one cm wave-length. 








Intensity relative to 


Frequency in 
strongest line 


megacycles 


BrCN 24,468 
24,506 
24,570 
24,576 
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TABLE II. Computed moments of inertia. 
































ee 
Moment of inerti 
Molecule (g cm? x10-#) 
Br?®CN 203.64 
BriCN 204.83 
CCN 140.47 
CR7CN 143.45 
Ocs 137.93 
TABLE III. Internuclear distances. 
Bond Bond distance f 
Bond distance electron diffraction! 
C-N 1.15A 1.13 +0.04A4 
Br-—C 1.79 1.79+0.02 
ci-—C 1.64 1.67 +0.02 














24,713 and 24,570 mc are attributed to molecules of Br™CN 
and Br®CN, respectively, in the ground vibrational state 
which make a transition from J=2 to J=3. The single 
OCS line is caused by the rotational transition J=1 to 
J=2. With these lines and assignments, and with no cor- 
rections for vibrations in the ground state, the moments of 
inertia are as in Table II. Using the isotopic frequency 
shift in CICN, the CN distance is calculated to be 1.15A, 
The isotopic shift in BrCN is also consistent with this value, 
but does not allow a very accurate determination of the 
CN distance. Assuming the CN distances in BrCN and 
CICN are the same, we obtain the values for the inter- 
nuclear distances given in Table III. 

The C—N distance determination is very sensitive to 
small changes in the isotopic shift. Until a satisfactory ex- 
planation of the subsidiary lines is obtained which will at 
the same time give assurance that no unknown effects are 
contributing to the isotopic shift, the internuclear distances 
given above must be regarded as tentative. 

The data can be used to distinguish between different 
possible structures for these molecules. The moments of 
inertia determined here and the electron diffraction data are 
consistent with the structures CI-C—N and Br—C-—N, 
but not with Cl—N—C and Br—N-—C which have some- 
times been proposed.? If one assumes an arrangement 
CI—N—C and takes 1.09A (the smallest reasonable value 
from electron diffraction data) as the N—C distance, the 
above value for the moment of inertia gives a total length 
for the molecule of 2.86A. This is in distinct disagreement 
with the electron diffraction value for the total length 
2.80+0.02A. From the above table it may be seen that the 
spectroscopic data give, assuming the Cl-C—N arrange- 
ment, a total molecular length 2.79A which is in good agree- 
ment with electron diffraction results. Similarly Br -C—N 
gives a total length 2.94A as compared with the electron 
diffraction value 2.92+0.02A, but if the Br—N-—C struc- 
ture is assumed the total length of the molecule cannot Le 
less than 3.01A. 

The frequency of the OCS line is in agreement with that 
recently reported by Dakin, Good, and Coles.‘ 

1C, H. Townes, Phys. Rev. 70, 665 (1946). 

2J. Y. Beach and A. Turkevich, J. Am. Chem. Soc. 61, 299 (1939). 

3 For a discussion see Si ick, Organic Chemistry of Nitrogen (Oxford 


University Press, New York, 1937), p. 326. 
( 4T. W. Dakin, W. E. Good, and D. K. Coles, Phys. Rev. 70, 560 
1946). 
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